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PREFACE 


This publication describes the IBM 
System/360 FORTRAN IV language for the 
IBM System/360 Operating System, the IBM 
System/360 Model 44 Programming System, 
and the IBM System/360 Disk Operating 
System. A reader should have some know- 
ledge of FORTRAN before using this pub- 
lication. A useful source for this 
information is the set of programmed 
instruction texts, FORTRAN IV for IBM 
System/360, Forms R29-0080 through 
R29-0087. 


The material in this publication is 
arranged to provide a quick definition 
and syntactical reference to the various 
elements of the language by means of a 
box format. In addition, sufficient 
text describing each element, with 
appropriate examples as to possible use, 
is given. 


Appendixes contain additional infor- 
mation useful in writing a FORTRAN IV 
program. This information consists of a 
table of source program characters, a 
list of other FORTRAN statements 
accepted by FORTRAN IV, a list of 
FORTRAN-supplied mathematical subpro- 
grams and service subprograms, lists of 
differences between FORTRAN IV and Basic 


Eighth Edition (October, 1968) 


FORTRAN IV and USA FORTRAN IV, and 
sample programs. Out-of-line mathemat- 
ical subprograms and service subprograms 
are described in the publication IBM 
System/360: FORTRAN IV Library Subpro- 
grams, Form C28-6596. Compiler restric- 
tions and programming considerations are 
contained in the programmer's guide for 
the respective system. The programmers' 
guides are as follows: 


Form C28-6639 


IBM System/360 Operating System: _ 
FORTRAN IV _(H) Programmer's Guide, 
Form C28-6602 


IBM System/360_ Model _ 44 Programming 
System: Guide to System Use for 
FORTRAN Programmers Form C28-6813 
IBM System/360 Disk Operating System, 


FORTRAN IV Programmer's Guide, Form 
C28-6397 


A comparison of FORTRAN IV compilers 


This is a major revision of, and makes obsolete, Form C28-6515-6. 
This edition clarifies text and corrects errors that appeared in the 
previous edition, and should be reviewed for revised, added, and 


deleted material. 


Changes to text, and small changes to illustrations, 


are indicated by a vertical line to the left of the change; changed or 
added illustrations are denoted by the symbol e to the left of the 


caption. 


Changes are periodically made to the specifications herein; 


any such 


changes will be reported in subsequent revisions or Technical 


Newsletters. 


Requests for copies of IBM publications should be made to your IBM 
representative or to the IBM branch office serving your locality. 


Address comments concerning the contents of this publication to IBM 
Corporation, Programming Publications, 1271 Avenue of the Americas, 


New York, N. Y. 
level being used. 


10020. 


Comments should mention the compiler and 
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INTRODUCTION 


IBM System/360 FORTRAN IV for the Operating System, the Model 44 Pro- 
gramming System, and the Disk Operating System consists of a language, a 
library of subprograms, and a compiler. 


The FORTRAN IV language is especially useful in writing programs for 
applications that involve mathematical computations and other manipula- 
tion of numerical data. The name FORTRAN is derived from FORmula 
TRANSlator. 


Source programs written in the FORTRAN IV language consist of a set 
of statements constructed by the programmer from the language elements 
described in this publication. 


In a process called compilation, a program called the FORTRAN com- 
piler analyzes the source program statements and translates them into a 
machine language program called the object program, which will be suit- 
able for execution on IBM System/360. In addition, when the FORTRAN 
compiler detects errors in the source program, it produces appropriate 
diagnostic error messages. The FORTRAN programmers’ guides listed in 
the Preface contain information about compiling and executing FORTRAN 
programs. 


The FORTRAN compiler operates under control of an operating system 
which provides the FORTRAN compiler with input/output and other ser- 
vices. Object programs generated by the FORTRAN compiler also operate 
under operating system control and depend on it for similar services. 


The IBM System/360 FORTRAN IV language is compatible with and encom- 
passes the United States of America (USA) FORTRAN, including its mathe- 
matical subroutine provisions. It also contains, as a proper subset, 
Basic FORTRAN IV. Appendixes F and G contain lists of differences 
between FORTRAN IV and BaSic FORTRAN IV and USA FORTRAN IV. 
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ELEMENTS OF THE LANGUAGE 


STATEMENTS 


Source programs consist of a set of statements from which the com- 


piler generates machine instructions, constants, and storage areas. A 
given FORTRAN statement effectively performs one of three functions: 


ie 


26 


Ss 


Causes certain operations to be performed (e.g., addition, multi- 
plication, branching) 


Specifies the nature of the data being handled 


Specifies the characteristics of the source program 


FORTRAN statements usually are composed of certain FORTRAN key words 


used in conjunction with the basic elements of the language: constants, 
variables, and expressions. The categories of FORTRAN statements are as 
follows: 


1. 


2e 


Arithmetic and Logical Assignment Statements: These statements 
cause calculations to be performed or conditions to be tested. The 
result replaces the current value of a designated variable or sub- 
scripted variable. 


Control Statements: These statements enable the user to govern the 
flow of and to terminate the execution of the object program. 


Input/Output Statements: These statements, in addition to control- 
ling input/output devices, enable the user to transfer data between 
internal storage and an input/output medium. 


FORMAT Statement: This statement is used in conjunction with cer- 
tain input/output statements to specify the form of a FORTRAN 
record. 


NAMELIST Statement: This statement is used in conjunction with 
certain input/output statements to specify the form of a special 
kind of record. 


DATA Initialization Statement: This statement is used to assign 


SSSR RRND SOD SEY TAGES SEY SEE EES TRIP RTP SNE SEE ND CARS SL GERD RE SONA SEEDED RE EE MID ESD METER eaemirD OEENS 


initial values to variables. 


Specification Statements: These statements are used to declare the 
properties of variables, arrays, and functions (such as type and 
amount of storage reserved) and, in addition, can be used to assign 
initial values to variables and arrays. 


Statement Function Definition Statement: This statement specifies 
Operations to be performed whenever the statement function name 
appears in the program. 


Subprogram Statements: These statements enable the user to name 
and define functions and subroutines, which can be compiled separ- 
ately or with the main program 


The basic elements of the language are discussed in this section 
The actual FORTRAN statements in which these elements are used are dis- 
cussed in following sections The term program unit refers to a main 
program or a subprogram; the term executable statements refers to those 
Statements in groups 1, 2, and 3. 


The order of a FORTRAN program unit (other than a BLOCK DATA subpro- 
gram) is as follows: 


1. Subprogram statement, if any. 
2.- IMPLICIT statement, if any. 


3. Other specification statements, if any. (Explicit specification 
statements that initialize variables or arrays must follow other 
specification statements that contain the same variable or array 
names. ) 


4. Statement function definitions, if any. 
5. Executable statements, at least one of which must be present. 
6. END statement. 


FORMAT, NAMELIST, and DATA statements may appear anywhere after the 
IMPLICIT statement, if present, and before the END statement. DATA 
statements, however, must follow any specification statements that con- 
tain the same variable or array names. A NAMELIST statement declaring a 
NAMELIST name must precede the use of that name in any input/output 
statement. 


The order of statements in BLOCK DATA subprograms is discussed in the 
section "BLOCK DATA Subprogram." 


CODING FORTRAN STATEMENTS 


The statements of a FORTRAN source program can be written on a stan- 
dard FORTRAN coding form, Form X28-7327. Each line on the coding form 
represents one 80-column card. FORTRAN statements are written one toa 
card within columns 7 through 72. If a statement is too long for one 
card, it may be continued on as many as 19 successive cards by placing 
any character, other than a blank or zero, in column 6 of each continua- 
tion card. For the first card of a statement, column 6 must be blank or 
ZCLO. 


As many blanks as desired may be written in a statement to improve 
its readability. They are ignored by the compiler. Blanks that are 
inserted in literal data are retained and treated as blanks within the 
data. 


Columns 1 through 5 of the first card of a statement may contain a 
statement number consisting of from 1 through 5 decimal digits. Blanks 
and leading zeros in a statement number are ignored. Statement numbers 
May appear anywhere in columns 1 through 5 and may be assigned in any 
order; the value of statement numbers does not affect the order in which 
the statements are executed in a FORTRAN program. : 


Columns 73 through 80 are not significant to the FORTRAN compiler and 


may, therefore, be used for program identification, sequencing, or any 
other purpose 
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Comments to explain the program may be written in columns 2 through 
80 of a card if the letter C is placed in column 1. Comments may appear 
between FORTRAN statements; a comments card may not immediately precede 
a continuation card. Comments are not processed by the FORTRAN com- 
piler, but are printed on the source program listing. Blanks may be 
inserted where desired to improve readability. 


CONSTANTS 


A constant is a fixed, unvarying quantity. There are four classes of 
constants -- those that specify numbers (numerical constants), those 
that specify truth values (logical constants), those that specify liter- 
al data (literal constants), and those that specify hexadecimal data 
(hexadecimal constants). 


Numerical constants may be integer, real, or complex numbers; logical 
constants may be «TRUE. or .FALSE.; literal constants may be a string 
of alphameric and/or special characters; and hexadecimal constants must 
be hexadecimal (base 16) numbers. 


INTEGER CONSTANTS 


] It occupies four locations of storage (i.e., four bytes). 
| 


Maximum Magnitude: 2147483647 (i.e., 231-1). 


An integer constant may be positive, zero, or negative; if unsigned, 
it is assumed to be positive. Its magnitude must not be greater than 
the maximum and it may not contain embedded commas. 


Examples: 


Valid Integer Constants: 


0 

941 

173 
—2147483647 


Invalid Integer Constants: 


27. (Contains a decimal point) 
3145903612 (Exceeds the allowable range) 
5, 396 (Contains an embedded comma) 
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REAL CONSTANTS 


Real Constant -- has one of three forms: a baSic real constant, a 


basic real constant followed by a decimal exponent, or an integer 
constant followed by a decimal exponent. 


mal point. If the string contains fewer than eight digits, the 


pies eight storage locations (bytes). 


The storage requirement (length) of a real constant can also be 
explicitly specified by appending an exponent to a basic real con- 


or the letter D followed by a Signed or unsigned 1- or 2-digit 
the letter D specifies a constant of length eight. 


Magnitude: (either four or eight locations) 0 or 167-®5 (approxi- 
: mately 10-78) through 16°63 (approximately 1075) 


Precision: (four locations) 6 hexadecimal digits 
(approximately 7.2 decimal digits) 
(eight locations) 14 hexadecimal digits 
(approximately 16.8 decimal digits) 


A real constant may be positive, zero, or negative (if unsigned, 
is assumed to be positive) and must be of the allowable magnitude. 
May not contain embedded commas. The decimal exponent permits the 


A basic real constant is a string Of decimal digits with a deci- 


basic real constant occupies four storage locations (bytes); if the 
String contains eight or more digits, the basic real constant occu- 


Stant or an integer constant. An exponent consists of the letter E 


integer constant. The letter E specifies a constant of length four; 


A A EE OT IE NS SD EY A AED EE A ND ED RRND NS SEE SEES GED SED AND ERED EN EA ETD EEE SAD TEED GED ERED GEE GAT EES CORED ED ERED CERED ERED SERED OUTED EEE GREP ED ERS EE CE GREED GEE ETD SES CUED cD GUE! GD GED QeEE cD ERD GED eee EEE coe aE a eee CED coe 


expression of a real constant as the product of a basic real constant or 


integer constant times 10 raised to a desired power. 


Examples: 


Valid Real Constants (four storage locations): 


+0. 

-~999.9999 

7. OE+0 (i1.e., 7.0 x 10° = 7.0) 

19761. 25E+1 (i.ee., 19761.25 x 101 = 197612. 5) 
TekS 

7.053 (1.e., 7.2.0 x 107 = 7000.0) 
7.O0F+03 

TE-03 {i.e., 7.0 x 107-3 = 0.007) 


Valid Real Constants (eight storage locations): 


1234567890123456.D-94 (Equivalent to .1234567890123456x10-78) 
21.98753829457168 

1.0000000 

7.9D03 

7.9D+03 (1.e€., 7.9 x 103 = 7900.0) 

7.9D+3 

7.9D0 (i.e., 7.9 x 10° = 7.9) 

7D03 (i.e., 7 0 x 103 = 7000. 0) 
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Invalid Real Constants: 


1 (Missing a decimal point or a decimal exponent) 
3,471.1 (Embedded comma) 
1.E (Missing a 1- or 2-digit integer constant fol- 


lowing the E. Note that it is not inter- 
preted as 10 x 10°) 


1. 2E+113 (E is followed by a 3-digit integer constant) 

23. 5E+97 (Magnitude outside the allowable range; that 
| is, 23.5 x 109751663) 

21. 3E-90 (Magnitude outside the allowable range; that 


is, 21.3 x 1079°<167-65) 


COMPLEX CONSTANTS 


Fr ge eS ge NG a a ey eG a eR wg Oe eee he ee ge 1 
] Definition | 
|--------------------------------------+-------------------------------- { 
|} Complex Constant - an ordered pair of signed or unsigned real con- | 
|] stants separated by a comma and enclosed in parentheses. The first | 
| real constant in a complex constant represents the real part of the | 
| complex number; the second represents the imaginary part of the com- | 
] plex number. Both parts must occupy the same number of storage | 
| locations (either four or eight). | 
peat gear Se SEN BN Qe ASU tm Ogee eB ga ep AN ye Nee PERO EE ee ET ee we J 


The real constants in a complex constant may be positive, zero, or 
negative (if unsigned, they are assumed to be positive), but they must 
be in the given range. 

Examples: 


Valid Complex Constants 


(3.2,-1.86) (Has the value 3.2 - 1.861) 

(-5. OE+03, .16E+02) (Has the value -5000. + 16.01) 
(4.7D+2,1.9736148) (Has the value 470. + 1.9736148i) 
(47D+2, 38D+3) (Has the value 4700. + 38000.i) 


Where i =\/=T 


(292704,1.697) (The real part is not a valid real constant) 
(.003E4,.005D6) (The parts differ in length) 
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LOGICAL CONSTANTS 


| Logical Constant - a constant that specifies a logical value There 
}] are two logical values: 

. TRUE. 
| 
| 


| 

| 

| 

- FALSE. | | 

Each occupies four storage locations. The words TRUE and FALSE must | 

| be preceded and followed by periods as shown above. | 
Ba a a a a a i ra J 


The logical constant .TRUE. or .FALSE. when assigned to a logical 
variable specifies that the value of the logical variable is true c~- 
false, respectively. (See the section “Logical Expressions.") 


LITERAL CONSTANTS 


a a ee 4 
] Definition 
SSS eS SRS as GC SNE I TE TO I RTC q 
] Literal Constant - a string of alphameric and/or special characters, | 
| delimited as follows: | 
| | 
| 1. The string can be enclosed in apostrophes. | 
| | 
| 2. The string can be preceded by wH where w is the number of | 
1 characters in the string. | 
L J 


Ae Ee cae en <> gh <n erme eae a come cern aD wenn ID ee tS gE A eS ET eng ee Ee OE ee ee eee a ae ee ee ee ee a ee ee ee ee ee ee ee ee ce ce ee ee ee ee ee ee ee ee ee ee eee ee ee ee 


The string may contain any characters (see Appendix A). Each 
character requires one byte of storage. The number of characters in the 
string, including blanks, may not be greater than 255. If apostrophes 
delimit the literal, a single apostrophe within the literal is repre- 
sented by two apostrophes. If wH precedes the literal, a single apos- 
trophe within the literal is represented as a single apostrophe. 


Literals can be used only in CALL statement or function reference 
argument lists, as data initialization values, or in FORMAT statements. 
The first form, a string enclosed in apostrophes, may be used in PAUSE 
statements. 


Examples: 


24H INPUT/OUTPUT AREA NO.2 

*DATA‘ 

* X-COORDINATE Y-~COORDINATE Z-COORDINATE® 
"3 14° 

A] DON® q Te 

5HDON'T 
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HEXADECIMAL CONSTANTS 


] Definition | 
ac ee ER ete a ae ee a ee Oe ee eee Oe Teen oN RE eh ete | 
| Hexadecimal Constant - the character Z followed by a hexadecimal | 
|} number formed from the set 0 through 9 and A through F. | 
a a as a re J 


Hexadecimal constants may be used only as data initialization values. 


One storage location (byte) contains two hexadecimal digits. Ifa 
constant is specified as an odd number of digits, a leading hexadecimal 
zero is added on the left to fill the storage location. The internal 
form of each hexadecimal digit is as follows: | 


0 —- 0000 4 —- 0100 8 - 1000 C - 1100 

1 - 0001 5 - 0101 9 - 1001 D- 1101 

2 -~ 9010 6 - 0110 A - 1010 E - 1110 

3 - 0011 7 - 0111 B - 1011 F - 1111 
Examples: 


Z1IC49A2F1 represents the bit string: 00011100010010011010001011110001 


ZBADFADE represents the bit string: 00001011101011011111101011011110 
where the first four zero bits are implied because an odd 
numker of hexadecimal digits is written. 


The maximum number of digits allowed in a hexadecimal constant 
depends upon the length specification of the variable being initialized 
(see "Variable Tyres and Lengths"). The following list shows the maxi- 
mum number of digits for each length specification: 


Length Specification , Maximum Number of 
of Variable Hexadecimal Digits 
16 32 
8 16 
4 8 
2 4 
1 2 


If the number of digits is greater than the maximum, the leftmost 
hexadecimal digits are truncated; if the number of digits is less than 
the maximum, hexadecimal zeros are supplied on the left. 


SYMBOLIC NAMES 


Fg er ne en er SE ee a eee gee ee 1 
} Definition | 


| Symbolic Name - from 1 through 6 alphameric (i.e., numeric, 0 

] through 9, or alphabetic, A through Z and $ ) characters, the first | 
] of which must be alphabetic. | 
L 


Symbolic names are used in a program unit (1.e€., a main program or a 
subprogram) to identify elements in the following classes. 


e An array and the elements of that array (see "Arrays") 
A variable (see "Variables"™) 


A statement function (see "Statement Functions") 


e An intrinsic function (see Appendix C) 


A FUNCTION subprogram (see “FUNCTION Subprograms") 


A SUBROUTINE subprogram (see “SUBROUTINE Subprograms") 


A block name (see “BLOCK DATA Subprogram") 


e An external procedure that cannot be classified as either a SUBROU- 
TINE or FUNCTION subprogram (see “EXTERNAL Statement") 


Symbolic names must be unique within a class in a program unit and 
can identify elements of only one class with the following exceptions. 


A block name can also be an array, variable, or statement function 
name in a program unit. 


A FUNCTION subprogram name must also be a variable name in the 
FUNCTION subprogram. 


Once a symbolic name is used as a FUNCTION subprogram name, a 
SUBROUTINE subprogram name, a block name, or an external procedure name 
in any unit of an executable program, no other program unit of that 
executable program can use that name to identify an entity of these 
classes in any other way. 


VARIABLES 


A FORTRAN variable is a symbolic representation of a quantity that 
occupies a storage area. The value specified by the name is always the 
current value stored in the area. 


For example, in the statement: 
A = 5.0+B 


both A and B are variables. The value of B is determined by some pre- 
vious statement and may change from time to time. The value of A is 
calculated whenever this statement is executed and changes as the value 
of B changes. 


VARIABLE NAMES 


The use of meaningful variable names can serve as an aid in document- 
ing a program. That is, someone other than the programmer may look at 
the program and understand its function For example, to compute the 
distance a car traveled in a certain amount of time at a given rate of 
speed, the following statement could have been written: | 
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X= Y * Z 


where * designates multiplication. However, it would be more meaningful 
to someone reading this statement if the programmer had written: 


DIST = RATE * TIME 


Examples: 


Valid Variable Names: 
B292S 
RATE 
SVAR 


Invalid Variakle Names: 


B292704 (Contains more than six characters) 
“ARRAY (First character is not alphabetic) 
SI.X (Contains a special character) 


VARIABLE TYPES AND LENGTHS 


The type of a variable corresponds to the type of data the variable 
represents. Thus, an integer variable represents integer data, a real 
variable represents real data, etc. There is no variable type asso- 
ciated with literal or hexadecimal data. These types of data are iden- 
tified by a name of one of the other types. 


For every type of variable, there is a corresponding standard and 
optional length specification which determines the number of storage 
locations (bytes) that are reserved for each variable. The following 
list shows each variable type with its associated standard and optional 
length: 


Variable Type Standard Optional 
Integer 4 2 
Real 4 8 
Complex 8 16 
Logical 4 1 


A programmer may declare the type of a variable by using the 
following: 


e Predefined specification contained in the FORTRAN language 

e Explicit specification statements 

e IMPLICIT statement 

An explicit specification statement overrides an IMPLICIT statement, 
which, in turn, overrides the predefined specification The optional 
length specification of a variable may be declared only by the IMPLICIT 
or explicit specification statements. If, in these statements, no 


length specification is stated, the standard length is assumed (see the 
section, "Type Statements"). 
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e Declaration by the Predefined Specification 


The predefined specification is a convention used to specify 
variables as integer or real as follows: 


1. If the first character of the variable name is I, J, K, L, M, or N, 
the variable is integer of a standard length 4. 


2. If the first character of the variable name is any other alphabetic 
character, the variable is real of a standard length 4. 


This convention is the traditional FORTRAN method of implicitly 
specifying the type of a variable as being either integer or real. In 
all examples that follow in this publication it is presumed that this 
specification applies unless otherwise noted. Variables defined with 
this convention are of standard length. 


Type Declaration by the IMPLICIT Statement 


The IMPLICIT statement allows a programmer to specify the type of 
variables in much the same way as was specified by the predefined con- 
vention. That is, in both the type is determined by the first character 
of the variable name. However, the programmer, uSing the IMPLICIT 
Statement, has the option of specifying which initial letters designate 
a particular variable type. The IMPLICIT statement can be used to 
specify all types of variables -- integer, real, complex, and logical -- 
and to indicate standard or optional length. 


The IMPLICIT statement overrides the variable type as determined by 
the predefined convention. For example, if the IMPLICIT statement 
specifies that variables beginning with the letters A through M are real 
variables and variables beginning with the letters N through Y are 
integer variakles, then the variable ITEM (which would be treated as an 
integer variable under the predefined convention) is now treated as a 
real variable. Note that variables beginning with the letters Z and $ 
are (by the predefined convention) treated as real variables. The 
IMPLICIT statement is presented in greater detail in the section "Speci- 
fication Statements." 


Type Declaration by Explicit Specification Statements 


ae wantin Salis mew ete ter SRN GRD ORES SEER Teme CNY UTED RD CONE EN SEEDED SRD GY NT AE RENSSS CN SY AES enane 


Explicit specification statements differ from the first two ways of 
specifying the type of a variable, in that an explicit specification 
statement declares the type of a particular variable by its name rather 
than as a group Of variables beginning with a particular character. 


For example, assume that an IMPLICIT statement overrode the prede- 
fined convention for variables beginning with the letter I by declaring 
them to be real and that a subsequent explicit specification statement _ 
declared that the variable named ITEM is complex. Then, the variable 
ITEM is complex and all other variables beginning with the character I 
are real. Note that variables beginning with the letters J through N 
are specified as integer by the predefined convention. 


The explicit specification statements are discussed in greater detail 
in the section "Specification Statements." 
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ARRAYS 


A FORTRAN array is a set of variables identified by a single variable 
name. A particular variable in the array may be referred to by its 
position in the array (e.g., first variable, third variable, seventh 
variable, etc.). Consider the array named NEXT which consists of five 
variables, each currently representing the following values: 273, 41, 
8976, 59, and 2. | 


NEXT(1) ais the location containing 273 
NEXT(2) is the location containing 41 
NEXT(3) is the location containing 8976 
NEXT(4) is the location containing 59 
NEXT(5) ais the location containing 2 


Each variable (element) in this array consists of the name of the 
array (1i.e., NEXT) immediately followed by a number enclosed in paren- 
theses, called a subscript quantity. The variables which the array con- 
prises are called subscripted variables. Therefore, the subscripted 
variable NEXT(1) has the value 273; the subscripted variable NEXT(2) has 
the value 41, etc. 


The subscripted variable NEXT(I) refers to the "Ith" subscripted _ 
variable in the array, where I is an integer variable that may assume a 
value of 1, 2, 3, 4, or 5. 


To refer to any element in an array, the array name must be sub- 
scripted. In particular, array name alone does not represent the first 
element. 


Consider the following array named LIST described by two subscript 
quantities, the first ranging from 1 through 5, the second from 1 
through 3: 


EY EEE CEES WORD UMN ua WHEE ETO 


Column 1 Column 2 Column 3 
7 


Row 1 82 u 

Row 2 12 13 14 
Row 3 91 1 31 
Row 4 24 16 10 
Row_5 2 8 2 


Suppose it is desired to refer to the number in row 2, column 3; this 
would be: 


LIST (2,3) 
Thus, LIST (2,3) has the value 14 and LIST (4,1) has the value 24. 
Ordinary mathematical notation might use LIST to represent any ele- 


ment of the array LIST. In FORTRAN, this is written as LIST(I,J) where 
I equals 1, 2, 3, 4, or 5 and J equals 1, 2, or 3. 


DECLARING THE SIZE AND TYPE OF AN ARRAY 


The size (number of elements) of an array is specified by the number 
of subscript quantities of the array and the maximum value of each sub- 
Script quantity. This information must be given for all arrays before 
using them in a FORTRAN program so that an appropriate amount of storage 
may be reserved Declaration of this information is made by a DIMENSION 
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Statement, a COMMON statement, or by one of the explicit specification 
statements; these statements are discussed in detail in the section 
"Specification Statements." The type of an array name is determined by 
the conventions for specifying the type of a variable name. Each ele- 
ment of an array is of the type specified for the array name. 


ARRANGEMENT OF ARRAYS IN STORAGE 


Arrays are stored in ascending storage locations, with the value of 
the eee, OF their subscript quantities increasing most rapidly and the 
value of the last increasing least rapidly. 


For example, the array LIST, whose values are given in the previous 
example, 1S arranged in storage as follows: 


82 12 91 24 24 13116 8 7 14 31 10 2 


The array named A, described by one subscript guantity which varies 
from 1 to 5, appears in storage as follows: 


AC1) AC2) AC3) AC4) AC5) 


The array named B, described by two subscript guantities with the 
first subscript quantity varying over the range from 1 to 5, and the 
second varying from 1 to 3, appears in ascending storage locations in 
the following order: 


B(1,1) B(2,1) B(3,1) B(4,1) B(5,1)--4 


a ee ee ome ee ee 2 eee ee oe ee ee ee ee ee ee ee ee ee ee ee ee ee es ee ee ee ee ee ee ee re ee eee re ee ee ee 


cr 
L->B(1,2) B(2,2) B(3,2) BC(4,2) B(5,2)—-, 


c 
L->B(1, 3) B(2,3) B(3,3) BC(4,3) B(5,3) 


Note that B(1,2) and B(1,3) follow in storage B(5,1) and B(5,2), 
respectively. 


The following list is the order of an array named C, described by 
three subscript quantities with the first varying from 1 to 3, the 
second varying from 1 to 2, and the third varying from 1 to 3: 


Cys CO 1s eG CG CO 2) 6,205 


OP AAD ES AS AGP A AE AE ENS ANS ED A ES A PS AP SS AS AE ES SE OS NS ND EE ED EY AY ES ET ES UE we EE ee ee ee A ED ee ee Se ee ee ee ee ee eee ae ee ee ee ee ee ee ee 


c 
4d 


OD Ee eee ee CE ED a ee ee ee ee a ED em ee a aD ee ee ee ee ee ee ee ee a es ee ee ee ee ee ee ee ee ee ee ee coe ee ee ee ee ee ee ee ee ee ee ee ee ee 


r 
L->c(1,1,3) C(2,1,3) C(3,1,3) C(1,2,3) C2,2,3) C(3,2,3) 


Note that C(1,1,2) and C(1,1,3) follow in storage C(3,2,1) and 
C(3,2,2), respectively. 


SUBSCRIPTS 


CGE CREO AED SY ENS CRORE NS AE 


A subscript is an integer subscript quantity or a set of integer sub- 
script guantities separated by commas, that is used to identify a parti- 
cular element of an array. The number of subscript quantities in any _ 
subscript must be the same as the number of dimensions of the array with 
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which the subscript is associated. 


A subscript is enclosed in paren- 
theses and is written immediately after the array name 


A maximum of 


seven subscript quantities can appear in a subscript. 


The following rules apply to the construction of subscript quanti- 


ties. (See the section "Expressions" 
the terms used below ) 


for additional information about 


1. Subscript quantities may contain arithmetic expressions that use 
any of the arithmetic operators: +t, -, *, 4, **. 

2. Subscript quantities may contain function references. 

3. Subscript quantities may contain subscripted names. 

4. Mixed mode expressions (integer and real only) within subscript 
quantities are evaluated according to normal FORTRAN rules. If the 
evaluated expression is real, it is converted to integer. 

5. The evaluated result of a subscript quantity should always be 
greater than zero and less than or equal to the size of the corres- 
ponding dimension. 

Examples: 
Valid Subscripted Variables: 
ARRAY (IHOLD) 
NEXT (19) 
MATRIX (I-5) 


BAK (I,J(K+1*L,.3*A(M,N))) 
ARRAY (I, d/4*K**2) 


ARRAY (-5) 
LOT (0) 


(A subscript quantity 

(A subscript quantity 
value of zero) 

(A subscript quantity 
false value) 

(A subscript quantity 
value) 


not be negative) 
never be nor assume a 


may 
may 
ALL(1.GE.TI) 


may not assume a true or 


NXT (1+(1.3,2.0)) may not assume a complex 


EXPRESSIONS 


FORTRAN IV provides two kinds of expressions: arithmetic and logi- 
cal. The value of an arithmetic expression is always a number whose type 
is integer, real, or complex. The value of a logical expression is 
always a truth value: .TRUE. or .FALSE.. Expressions may appear in 
assignment statements and in certain control statements. 


ARITHMETIC EXPRESSIONS 


The simplest arithmetic expression consists of a primary which may be 
a Single constant, variable, subscripted variable, function reference, 
or another expression enclosed in parentheses. The primary may be eith- 
er integer, real, or complex. | 
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In an expression consisting of a single primary, the type of the pri- 


Mary is the type of the expression 


Examples: 
Primary Type of Primary Type of Expression 
3 Integer constant Integer of length 4 
A Real variable Real of length 4 
3.14D3 Real constant Real of length 8 
(2.0,5.7) Complex constant Complex of length 8 
SIN(X) Real function reference Real of length 4 
(A*B+C) Parenthesized real Real of length 4 


expression 


More complicated arithmetic expressions containing two or more pri- 


maries may be formed by using arithmetic operators that express the 
computation(s) to be performed. 


Arithmetic Operators 


The arithmetic operators are as follows: 


Arithmetic Operator Definition 
* * Exponentiation 
* Multiplication 
/ Division 
+ Addition 
_ Subtraction 


RULES FOR CONSTRUCTING ARITHMETIC EXPRESSIONS: The following are the 
rules for constructing arithmetic expressions that contain arithmetic 
operators: 


1. 


All desired computations must be specified explicitly. That is, if 
more than one primary appears in an arithmetic expression, they 
must be separated from one another by an arithmetic operator. For 
example, the two variables A and B will not be multiplied if 
written: 

AB 


If multiplication is desired, the expression must be written as 
follows: 


A*B or B*A 


No two arithmetic operators may appear in sequence in the same 
expression. For example, the following expressions are invalid: 


A*/B and A*-B 
The expression A*-B could be written correctly as follows: 
A* (-B) 
In effect, -B will be evaluated first and then A will be multiplied 


with it. (For further uses of parentheses, see rule 3.) 
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Order of Computation: Computation is performed from left to right 
according to the hierarchy of operations shown in the following 
list 


Operation Hierarchy 
Evaluation of functions 1st 
Exponentiation (**) 2nd 
Multiplication and division (* and /) 3rd 
Addition and subtraction (+ and -) GY“th 


This hierarchy is used to determine which of two consecutive opera- 
tions is performed first. If the first operator is higher than or 
equal to the second, the first operation is performed. If not, the 
second operator is compared to the third, etc. When the end of the 
expression is encountered, all of the remaining operations are per- 
formed in reverse order. 


For example, in the expression A*B+C*D**I, the operations are per- 
formed in the following order: 


ae A*B Call the result X (multiplication) (X+C¥D*¥*T) 
be. D*¥*I Call the result Y (exponentiation) (X+C¥Y) 


Cc. CY Call the result Z (multiplication) (X+Z) 
de X+Z Final operation (addition) | 


If there are consecutive exponentiation operators, the evaluation 
is from right to left. Thus, the expression: 
A*¥*B*¥*C 
is evaluated as follows: 
ae B*¥*C Call the result Z 


be A**Z Final operation 


A unary plus or minus has the same hierarchy as a plus or minus in 
addition or subtraction. Thus, 


A=-B is treated as A=0-B 
A=-B*C is treated as A=-(B*C) 
=-B+C is treated as A=(-B)+C 

Parentheses may be used in arithmetic expressions, as in algebra, 
to specify the order in which the arithmetic operations are to be 
computed. Where parentheses are used, the expression within the 
parentheses is evaluated before the resuit is used. This is equi- 
valent to the definition above since a parenthesized expression is 
a primary. 
For example, the following expression: 


Bt ( (A+B) *C) +A*¥*2 


is effectively evaluated in the following order: 


ae (A+B) Call the result X B+ (X¥*C) +A*#*2 
be (X*C) Call the result Y BtYt+A**2 

Cc. Bty Call the result W WtA**2 

ad. A**2 Call the result Z W+Z 

e. W+tZ Final operation 


Table 1. Determining the Type and Length of the Result of +, -, *, / 
Operations 
i ara aa aa areas aaa ee aa oS eae oa a Ne eee ae ge 1 
| | | | | | 
JINTEGER | INTEGER ] REAL | REAL | COMPLEX] COMPLEX| 
j+ - * / } (2) | (4) ] (4) } (8) | (8) | (16) | 
----------- }-----~--4}---------}---------}---------4--------+-------- 
| INTEGER JInteger ] Integer | Real | Real {| Complex] Complex] 
| (2) } (2) | (4) ] (4) } (8) | (8) | (16) | 
eke ean : aaa: GLE eRe: Gan RG: (oh EA ese aa Gea ania 
{| INTEGER {Integer |] Integer | Real | Real {| Complex] Complex] 
] (4) } (4) | (4) } (4) } (8) | (8) | (16) | 
~----------- }--------4}---------4---------4---------4--------4-------- 
] REAL }Real | Real | Real { Real | Complex] Complex} 
| (4) } (4a) ] (4) | (4) | (8) | (8) | (16) | 
~----------- }--------4}---------4---------}---------4--------4-------- 
] REAL ]Real ] Real | Real |} Real | Complex] Complex] 
| (8) } (8) (8) (8) } (8) } (16) | (16) | 
~~---------- }-------- 4--------- 4--------- }---------4-------- 4-------- { 
] COMPLEX [Complex | Complex |] Complex | Complex | Complex] Complex} 
| (8) } (8) } (8) | (8) } (16) | (8) | (16) | 
{------------- }-------- --------- 4--------- }--------- -------- -------- 
|] COMPLEX {Complex | Complex | Complex | Complex |] Complex]| Complex} 
| (16) } (16) } (16) } (16) } (16) {| (16) | (16) | 
|—-————_~_.~—~— Detect a eg ige ee Eien peters te Pree os ane ene peepee een eae pe ame ener es 
|Note: When division is performed using two integers, the answer is 
Jtruncated and an integer answer is given. For example, 


]J=2, then the expression (I/J) would yield an integer answer of 4 


Jafter truncation. 


| 
if I=9 and | 
| 
| 


The type and length of the result of an operation depends upon the 
type and length of the two operands (primaries) involved in the | 
Operation. Table 1 shows the type and length of the result of the 
operations +, -, *, and /. 


Assume that the type of the following variables has been specified 
as follows: 


Variable Names Type Length Specification 
Cc Real variable 4 

Ey dy *K Integer variable ue, 2, 2 

D Complex variable 16 


Then the expression I*J/C**K+D is evaluated as follows: 


Subexpression Type _ and Length 
I*¥J (Call the result X) Integer of length 4 
C¥*K (Call the result Y) Real of length 4 
X/Y (Call the result Z) Real of length 4 
Z+D Complex of length 16 


Thus, the final type of the entire expression is complex of length 
16, but the type changed at different stages in the evaluation. 
Note that, depending on the values of the variables involved, the 
result of the expression I*J*C might be different from I*C*J. 
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5. The arithmetic operator denoting exponentiation (i e ,**) may only 
be used to combine the types of operands shown in Table 2. 


The type of the result depends upon the type of the two operands 


involved, as shown in Table 1. For example, if an integer is 
raised to a real power, the type of the result is real. 


Table 2. Valid Combinations with the Arithmetic Operator ** 


aa a aaa ae ama aa ae ae aa 1 
1 Base Exponent | 
}~------~------------~~--~-~-----------------------~--------------------- { 
] Integer (either length) Integer (either length) | 
| or + Or | 
1 Real (either length) Real (either length) | 
| 
| Complex (either length) ** Integer (either length) | 
Oe ah a a a ae ee ee a ee J 


LOGICAL EXPRESSIONS 


The simplest form of logical expression consists of a single logical 
primary, which can be a logical constant, logical variable, logical sub- 
Scripted variable, logical function reference, or logical expression 
enclosed in parentheses, which always has the value .TRUE. or .FALSE.. 


More complicated logical expressions may be formed by using logical 
and relational operators. These expressions may be in one of the fol- 


lowing forms: 


1. Relational operators combined with arithmetic expressions whose 
type is integer or real. 


2e Logical operators combined with logical primary. 


3. Logical operators combined with either or both forms of the logical 
expressions described in items 1 and 2. 


Item 1 is discussed in the following section, "Relational Operators;" 
items 2 and 3 are discussed in the section “Logical Operators." 


Relational Operators 


The six relational operators, each of which must be preceded and fol- 
lowed by a period, are as follows: | 


Relational Operator Definition 

-GT. Greater than (>) 

- GE. Greater than or equal to (2) 
- LT. Less than (<) 

- LE. Less than or equal to (S<) 

oe EQ. Equal to (=) 

» NE. Not equal to (#) 


The relational operators express an arithmetic condition which can be 
either true or false Only arithmetic expressions whose type is integer 
Or real may be corbined by relational operators. For example, assume 
that the type of the following variables has been specified as follows: 
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Variable Names Type 


ROOT, E Real variables 

A, I, F Integer variables 
L Logical variable 
Cc Complex variable 


Then the following examples illustrate valid and invalid logical 
expressions using the relational operators. 


Examples: 


Valid Logical Expressions Using Relational Operators: 
A .LT. I 

E**2.7 «EQ. (5*ROOT+4) 

~5 GE. .9*ROOCT 

E «EQ. 27.3D+05 


Invalid Logical Expressions Using Relational Operators: 


C .GE. (2.7,5.9E3) (Complex quantities may never appear in logical 


expressions) 
L .-EQ. (AtF) (Logical quantities may never be joined by 
relational operators) 
E**2 .EQ 97.1E9 (Missing period immediately after the relational 
operator) 
-GT. 9 (Missing arithmetic expression before the rela- 


tional operator) 


Logical Operators 


The three logical operators, each of which must be preceded and fol- 
lowed by a period, are as follows (where A and B represent logical con- 
stants or variables, or expressions containing relational operators): 


Logical Operator Use Meaning 
- NOT. -NOT.A If Ais .TRUE., then .NOT.A has the value 


-FALSE.; if A is .FALSE., then .NOT.A 
has the value .TRUE. 


e AND. A.AND.B If A and B are both .TRUE., then A.AND.B 
has the value .TRUE.; if either A or B 
or both are .FALSE., then A.AND.B has 
the value .FALSE. 


«OR. A.OR.B If either A or B or both are .TRUE., then 
A.OR.B has the value .TRUE.; if both A 
and B are .FALSE., then A.OR.B has the 
value .FALSE. 


Two logical operators may appear in sequence only if the second one 
is the logical operator .NOT.. 


Only those expressions which, when evaluated, have the value .TRUE. 
or «FALSE may be combined with the logical operators to form logical 
expressions. For example, assume that the type of the following vari- 
ables has been specified as follows: 
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Variable Names Type 


ROOT, E Real variables 

A, I, F Integer variables 
L, W Logical variables 
C Complex variable 


Then the following examples illustrate valid and invalid logical 
expressions using Loth logical and relational operators. 


Examples: 


Valid Logical Expressions: 


(ROOT*A .GT. A) «AND. W 
L e-AND. .NOT. (I .GT. F) 
(E+5.9D2 .GT. 2*E) .OR. L 
-NOT. W .AND. .NOT. L 
L «AND. «NOT. W.OR. I .GT. F 
(A**F .GT. RCCT) .-AND. .NOT. (1 .EQ. EB) 


Invalid Logical Expressions: 


A .AND. L (A is not a logical expression) 
-OR. W (.OR. must be preceded by a logical expression) 
NOT. (A .GT. F) (Missing period before the logical operator 
- NOT.) | 


(C .EQ. I) .AND. L (A complex quantity may never be an operand of 
a relational operator) 


L .AND. .OR. W (The logical operators .AND. and .OR. must 
always be separated by a logical expression) 
eAND. L. (.AND. must be preceded by a logical 
| expression) 


Order of Computations_in Logical Expressions: The order in which the 
operations are performed is: 


Operation Hierarchy 
Evaluation of functions 1st (highest) 
Exponentiation (¥**) 2nd | 
Multiplication and division (* and /) 3rd 

Addition and subtraction (+ and -) 4th 

UT eyo hej BOs ye NEsgs Gls 7 eGE. 5th 

~ NOT. 6th 

«AND. 7th 

OR. 8th 
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For example, the expression: 


A.GT.D**B. AND. .NOT.L.OR.N 


is effectively evaluated in the following order: 


1. D**B Call the result W (exponentiation) 

2. A.GT.W Call the result X (relational operator) 

3. .NOT.L Call the result Y (highest logical operator) 

4. X.AND.Y Call the result 2 (second highest logical operator) 
5. ZeOR.N Final operation 


Note: Logical expressions may not require that all parts be evaluated. 
Functions within logical expressions may or may not be called. For 
example, in the expression A.OR.LGF(.TRUE.), it should not be assumed 
that the LGF function is always invoked. 


Use of Parentheses in Logical Expressions: Parentheses may be used in 
logical expressions to specify the order in which the operations are to 
be performed. Where parentheses are used, the expression contained 
within the most deeply nested parentheses (that is, the innermost pair 
of parentheses) is effectively evaluated first. For example, the logi- 
cal expression: 


(I.GT. (B+C)) .AND.L 
is effectively evaluated in the following order: 


1.2 Btc Call the result X 
ac Le GlaX Call the result Y 
3. Y.AND.L Final operation 


The logical expression to which the logical operator .NOT. applies 
must be enclosed in parentheses if it contains two or more quantities. 
For example, assume that the values of the logical variables, A and B, 
are .FALSE. and .TRUE., respectively. Then the following two expres- 
sions are not equivalent: 


« NOT. (A. OR. B) 
-NOT.A.~OR.B 


In the first expression, A.OR.B, is evaluated first. The result is 
~TRUE.; but .~NOT.(.TRUE.) implies .FALSE.. Therefore, the value of the 
first expression is .FALSE. 


In the second expression, .NOT.A is evaluated first. The result is 


-TRUE.; but .~.TRUE..OR.~B implies .TRUE.. Therefore, the value of the 
second expression is .TRUE.. 
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ARITHMETIC AND LOGICAL ASSIGNMENT STATEMENT 


ee a a ae a ee 1 
] General Form | 
|---------------------------------------------------------------------- { 
| a=b | 
| ; | | 
] Where: a is a subscripted or nonsubscripted variable. | 
| 
| b is an arithmetic expression or logical expression. 
De aaa a ate Sn a ala ee te ek eta Ra J 


This FORTRAN statement closely resembles a conventional algebraic 
equation; however, the equal sign specifies replacement rather than 
equivalence. That is, the expression to the right of the equal sign is 
evaluated, and the resulting value replaces the current value of the 
variable to the left of the equal sign. 


If b is a logical expression, a must be a logical variable. If b is 
an arithmetic expression, a must be an integer, real, or complex vari- 
able. Table 3 gives the conversion rules used for placing the evaluated 
result of arithmetic expression b into variable a. 


Assume that the type of the following variables has been specified 


as: 
Variable Names ‘Type Length Specification 
I, J, W Integer variables 4,4,2 

Ay. By Cz. D Real variables 4,4,8,8 

E Complex variable 8 

G, H Logical variables 4,4 


Then the following examples illustrate valid arithmetic statements 
using constants, variables, and subscripted variables of different 
types: 


Statements Description 
A = The value of A is replaced by the current value of RB. 
W=B The value of B is truncated to an integer value, and 


this value replaces the value of W. 


A= I The value of I is converted to a real value, and this 
result replaces the value of A. 


I=.iI+i The value of I is replaced by the value of I + 1. 


ty 
I 


I**J+D I is raised to the power J and the result is con- 
verted to a real value to which the value of D is 
added. This result replaces the real part of the 
complex variable E. The imaginary part of the con- 
plex variable is set to zero. 


A = C*¥D The most significant part of the product of C and D 
replaces the value of A. 


A=E The real part of the complex variable E replaces the 
value of A. 
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Statements 
E= A 


G = .TRUE. 


H = .NOT.G 


G = 3..GT.I 


E= (1.0,2.0) 


Description | 
The value of A replaces the value of the real part of 


the complex variable E; the imaginary part is set 
equal to zero 


The value of G is replaced by the logical constant 
« TRUE. .« 


If G is .TRUE., the value of H is replaced by the 
logical constant .FALSE.. If G is .FALSE., the 
value of H is replaced by the logical constant 
- TRUE.. 


The value of I is converted to a real value; if the 
real constant 3. is greater than this result, the 
logical constant .TRUE. replaces the value of G. 
If 3. is not greater than I, the logical constant 
~FALSE. replaces the value of G. 


The value of the complex variable E is replaced by 
the complex constant (1.0,2.0). Note that the 
statement E = (A,B) where A and B are real vari- 
ables is invalid. 
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Table 3. 

Sr aoa 7 

| 

| | INTEGER*2 
| INTEGER*4 
| | 

Le £2, Ee ahaa EEE 
] INTEGER*2 Assign 

] INTEGER*4 

| 
}-_.______f#______ 








Ee EEE CE I OE SUNS EES SEY CY ER 


| 
| 
| 
4 

REAL*4 | Float and 
| assign 
| 
| 
4 


REAL*8 [DP Float 
jand assign 


4 
| COMPLEX*8 |Float and 


Jassign to 


Assign to 
real part; 


= 
{Real assign 
}real part; 


Conversion Rules for the Arithmetic Assignment Statement a = b 


a ere eee aaa a a aaa a | 
| | | | | 
| | | | | 
| REAL#*4 {| REAL*8 | COMPLEX*8 ]COMPLEX*16 | 
| | | | | 
}----------- 4—---------- 4----------- 1——--------- { 
| Fix and assign | Fix and assign real | 
| | part; imaginary part | 
| | not used. | 
+----------- ----------- 4----------- q----------- { 
| Assign }Real assign|Assign real|Real assign| 
| | |part; imag-]real part; | 
] | [inary part |imaginary | 
| | [not used. |part not | 
| | | used. | 

Sah tae ee Jeanne i fe 1] 
| Assign |Assign real part; imag-| 
| Jinary part not used. | 

~----~----}----------- }~----------7-----------}-----------7-----------4 


+ 
{| Assign }Real assign| 
| [real and 


| | 
] J}xreal part; Jimaginary |imaginary | [imaginary | 
| Jimaginary |part set [part set to| Jparts. 
Jpart set [to zero. |zero. | | | 
l jto zero. | | | | | 
pe ere 1 a aaa ae | teri a a Sait eR | 
] COMPLEX*16]DP float |] Assign to real part; | Assign | 
| Jand assign | imaginary part set to | | 
| jto real |} zero. | | 
Jpart; imag-| | | 
] Jinary part | | 
] JjJset to | | | 
| | ZEKO. | | | 
}--------—~— ele fee ea a 5 Seen ce eae een ey ey en Oe 8 Ie aye eg ae cr See reed ed Pek ere | 
| Notes: | 
]1. Assign means transmit the resulting value, without change. If | 
] the significant digits of the resulting value exceed the speci- | 
l fied length, results are unpredictable. | 
]2. Real Assign means transmit to a as much precision of the most | 
Significant part of the resulting value as REAL*4 data can | 
l contain. | 
]}3. Fix means transform the resulting value to the form of a basic | 
] real constant and truncate the fractional portion. | 
{4. Float means transform the resulting value to the form of a REAL*4 | 
| number, retaining in the process as much precision of the value | 
| aS a REAL*4 number can contain. | 
J5. DP Float means transform the resulting value to the form of a | 
| REAL*8 number. | 
]6. An expression of the form E=(A,B), where E is a complex variable | 
| and A and B are real variables, is invalid. The mathematical | 
j function subprogram CMPLX can be used for this purpose. See | 
Appendix C. 
Ma a ee a J 
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CONTROL STATEMENTS 


Normally, FORTRAN statements are executed sequentially. That is, 
after one statement has been executed, the statement immediately follow- 
ing it is executed. This section discusses the statements that may be 
used to alter and control the normal sequence of execution of statements 
in the program. 


GO TO STATEMENTS 


GO TO statements permit transfer Of control to an executable state- 
ment specified by number in the GO TO statement. Control may be trans- 
ferred either unconditionally or conditionally. The GO TO statements 
are: 


1. Unconditional GO TO statement 
2e Computed GO TO statement 


3. Assigned GO TO statement 


UNCONDITIONAL GO TO STATEMENT 


| GO TO xxxxx | 
| | 


} Where: xxxxx iS an executable statement number. | 


This GO TO statement causes control to be transferred to the state- 
ment specified by the statement number. Every subsequent execution of 
this GO TO statement results in a transfer to that same statement. Any 
executable statement immediately following this statement should have a 
statement number; otherwise it can never be referred to or executed. 


Example: 


GO TO 25 
10 A=B+#t#C 


Explanation: 


In this example, each time the GO TO statement is executed, control 
is transferred to statement 25. 
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COMPUTED GO TO STATEMENT 


<P eee en ee Ee ee eee eee ee ee ee ee ED Gee ee eee le eee ee ee op ee ee ee > ee ee ee ee ee oe ee ee ee ee ee ee ee oe ee ee ee ee oe eee ee ee ee ee ee es ee ee ee ee ee ee ee ee ce ee ee ee ee ee ee oe 


ee te: tree eee ee ee cy ee ce ee ee ee ee ee i ce ee Ee ee ee ee eee ee ee ee ee ce ee ee ge ee ee ee ee ee es ee ee ce ee ee es a es ee ee ee ee ss ee ee ee es ee ee ee ee ee ee ee ee ce ee ee es ee ee ee 


i is a nonsubscripted integer variable whose current value 


] Where: X4,X2a,e5eeXne are executable statement numbers. 
| is in the range: 1<i«<n 


This statement causes control to be transferred to the statement num- 
bered Xa, X2ay X3¢e0ee, OF Xn, Gepending on whether the current value of i 
is 1, 2, 3,.--, Or n, respectively. If the value of i is outside the 
allowable range, the next statement is executed. 


Example: 
GO TO (25, 10, 7), ITEM 


7 C = E*¥*2+A 
25 L=C 


10 B = 21.3E02 
Explanation: 


In this example, if the value of the integer variable ITEM is 1, 
Statement 25 will be executed next. If ITEM is equal to 2, statement 10 
is executed next, and so on. 


ASSIGN AND ASSIGNED GO TO STATEMENTS 


A ND CERT A EE EE EE ONS EE NEE EAD AT ED ND NT ED AE ED AD EE ED AOS EES ANS ND ND ERD AED ES ENS RED SEED SEED NN AEN IED COURS AR ENED ERED GANGY RED OEE SENET ENED GED REE EES GEL CHEE RTD GEE GUNES calnieD HUET GRE) GEER GENES CCONET GENIE CHEM ciate ANNE dinieee QUERY <i <clnD ED me 


GO TO my (X1—¢Xa—Xa¢ee%¢Xn) 


| 

| 

| 

| 

] Where: i is an executable statement number. It must be one of the 
| numberS X4¢Xa¢eX3ae0%%,Xne 
| 

| 

| 
L 


Xa—XaeXaeeee,Xn are executable statement numbers in the pro- 
gram unit containing the GO TO statement. 


m is a nonsubscripted integer variable of length 4 which is 
assigned one of the statement numberS: X4,Xa,sX3e9< «¢Xne 


The assigned GC TO statement causes control to be transferred to the 
Statement numbered X4,Xa,Xae «,OX Xn, Aepending on whether the current 
assignment of m iS X4,X2a—eX3ye00e,OX Xn, respectively. For example, in 
the following statement: 


GO TO N, (10, 25, 8) 


If the current assignment of the integer variable N is statement number 
8, then the statement numbered 8 is executed next. If the current 
assignment of N is statement number 10, the statement numbered 10 is 
executed next. If N is assigned statement number 25, statement 25 is 
executed next. 


At the time of execution of an assigned GO TO statement, the current 
value of m must have been defined to be one of the values x, Xa-.-Xn by 
the previous execution of an ASSIGN statement. The value of the integer 
variable m is not the integer statement number; ASSIGN 10 TO I is not 


the same as L120, 


Example 1: 


ASSIGN 50 TO NUMBER 
10 GO TO NUMBER, (35, 50, 25, 12, 18) 


Explanation: 


In example 1, statement 50 is executed immediately after statement 
10. 


Example 2: 


ASSIGN 10 TO ITEM 


13 GO TO ITEM, (8, 12, 25, 50, 10) 


10 B=C + D 
ASSIGN 25 TO ITEM 
GO TO 13 


25 C = E*¥*¥2 
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Explanation: 
In example 2, the first time statement 13 is executed, control is 


transferred to statement 10. On the second execution of statement 13, 
control is transferred to statement 25 | 


ADDITIONAL _CONTROL_ STATEMENTS 


ARITHMETIC IF STATEMENT 


A A SD EP EY AS EP ES AY a EN ee ES <a aN ENS AEP ES AE ED “ORE Eee cay GEES GSE <i Ge eae GE que <i <iiee carey <UnR aap SD nee ce GE ae Se ee ae le ee ee ee ae Es ee ee ee cee re ee ee i ee ee ce se ee ee cee ee ee oe ee 


| 

| 

] Where: ais any arithmetic expression except complex. 
! 

| 

L 


Xa~9XasX3 are executable statement numbers. 


AP ED A OED ND AE ED AD ED NS SD RE ERED RTD RY SD INS SEND ERED ES SEED NER GETS SEED ANE “EAD IRD AD ERISY RENTS EIS SET ENED INNS ERE EY ET NED SN NY “ED SUA INNS SNEED SEIN LIEN NDS ERR ONE “SEE AD GUE ERIS SESE ENED CED UB GANLEY SU ERED NED <ONIRD CORNY ENED <DUER ony SERED CUTE Ee 


The arithmetic IF statement causes control to be transferred to the 
statement numbered x4,X2, OXF X3 when the value of the arithmetic expres- 
Sion (a) is less than, equal to, or greater than zero, respectively. 

The first executable statement following the arithmetic IF statement 
should have a statement number; otherwise, it can never be referred to 
or executed. 


Example: 


10 E (F*B)/D+1 


Explanation: 


In this example, if the value of the expression (A(J,K)**3-B) is 
negative, the statement numbered 10 is executed next. If the value of 
the expression is zero, the statement numbered 4 is executed next. If 
the value of the expression is positive, the statement numbered 30 is 
executed next. 
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LOGICAL IF STATEMENT 


Te AEE et ND ED ERP ED ED eee CE ED EY ED ete ES ENS eR ES te RE EE ES ERED EES ETE EEE ETD ED ED CEE CEEED ES ee EE Gee ee aes AEE ee ae em eeeD Snee a ee Te ee Ee ee ee ee eee ee ee ee oe ee ee ee ee ee ee ee ee ee ee oe oe ee 


Ss is any executable statement except a DO statement or 


| 

) ; 

1 Where: ais any logical expression. 
| another logical IF statement. 


The logical IF statement is used to evaluate the logical expression 
(a) and to execute or skip statement s depending on whether the value of 
the expression is .TRUE. or .FALSE., respectively. 


Example 1: 


IF(A.LE.0.0) GO TO 25 


C=DteE 
IF(A.EQ.B) ANSWER = 2.0*A/C 
F = G/H 


Explanation: 


In the first statement, if the value of the expression is .TRUE. 
(i.e., A is less than or equal to 0.0), the statement GO TO 25 is 
executed next and control is passed to the statement numbered 25. If 
the value of the expression is .FALSE.(i.e., A is greater than 0.0), the 
statement GO TO 25 is ignored and control is passed to the second 
statement. 


In the third statement, if the value of the expression is .TRUE. 
(i1.e., A is equal to B), the value of ANSWER is replaced by the value of 
the expression (2.0*A/C) and then the fourth statement is executed. If 


the value of the expression is .FALSE. (i.e., A is not equal to B), the 


value of ANSWER remains unchanged and the fourth statement is executed 
next. 


Example 2: 


Assume that P and Q are logical variables. 


IF(P.OR.~.«NOT Q)A=B 
C = B¥F*2 
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Explanation: 


In the first statement, if the value of the expression is TRUE, the 
value of A is replaced by the value of B and the second statement is 
executed next. If the value of the expression is .FALSE., the statement 
A = B is skipped and the second statement is executed. 


DO STATEMENT 


{| General Form 
| a ene + 
End of DO Initial Test 
Range Variable Value Value Increment 
a Se ~~ —_—_— Se ee 
DO x 5. = Ma ¢ Ma2¢ Ms 


Where: x is an executable statement number appearing after the DO 
statement. 


i is a nonsubscripted integer variable. 


M4, Moa, and ms, are either unsigned integer constants great- 
er than zero or unsigned nonsubscripted integer variables 
whose value is greater than zero. Mz may not exceed 231-2 
in value. m3 is optional; if it is omitted, its value is 
assumed to be 1. In this case, the preceding comma must 
also be omitted. 


RY Re <a ED RN AEE EE aie ED eel EN EY ee ES ED NE emer cline deme imate Glline NEY am AENED anetD SERS Ee ED REET GES ED ERD GED ED ERY GUESSES ENE AO ED ED ED A EE RD aE GD GREE UT CED GEES cme ee RED Ge Gems CED CEP Ee eae ee Gee ae Ce eee em ame ome oe 


pe Sa eee te re ree ee ee oe aD Oe Gece GE eee eee Gees ose 


The DO statement is a command to execute, at least once, the state- 
ments that physically follow the DO statement, up to and including the 
statement numbered x. These statements are called the range of the Do. 
The first time the statements in the range of the DO are executed, i is 
initialized to the value mi; each succeeding time i is increased by the 
value m3- When, at the end of the iteration, i is equal to the highest 
value that does not exceed ma, control passes to the statement following 
the statement numbered x. Thus, the number of times the statements in 
the range of the DO are executed is given by the expression: 


where the brackets represent the largest integral value not exceeding 
the value of the expression within the brackets. If mz is less than m,, 
the statements in the range of the DO are executed once. Upon comple- 
tion of the DO, the DO variable is undefined and may not be used until 
assigned a value (e.g., in a READ list). 


There are several ways in which looping (repetitively executing the 
same statements) may be accomplished when using the FORTRAN language. 
For example, assume that a manufacturer carries 1000 different machine 
parts in stock. Periodically, he may find it necessary to compute the 
amount of each different part presently available. This amount may be 
calculated by subtracting the number of each item used, OUT(I), from the 
previous stock on hand, STOCK(I). 
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Example 1: 


T=0 

10 %I=1+1 
STOCK (I)=STOCK(1I)- OUT(I) 
IF(I-1000) 10,30, 30 

30 A=B+C 


The first, second, and fourth statements required to control the pre- 
viously shown loop could be replaced by a single DO statement as shown 
in example 2. 


Example 2: 


DO 25 I = 1,1000 
25 STOCK(I) = STOCK(I)-OUT(IL) 
A = B+tC 


In example 2, the DO variable, I, is set to the initial value of 1. 
Before the second execution of statement 25, I is increased by the 
increment, 1, and statement 25 is again executed. After 1000 executions 
of the DO loop, I equals 1000. Since I is now equal to the highest 
value that does not exceed the test value, 1000, control passes out of 


the DO loop and the third statement is executed next. Note that the DO 
variable I is now undefined; its value is not necessarily 1000 or 1001. 


Example 3: 
DO 25 I=1, 10, 2 
J = Itk 
25 ARRAY(J) = BRAY(J) 
A= B+C 


Explanation: 


In example 3, statement 25 is the end of the range of the DO loop. 
The DO variable, I, is set to the initial value of 1. Before the second 
execution of the DC loop, I is increased by the increment, 2, and the 
second and third statements are executed a second time. After the fifth 
execution of the DO loop, I equals 9 Since I is now equal to the high- 
est value that does not exceed the test value, 10, control passes out of 
the DO loop and the fourth statement is executed next. Note that the DO 


variable I is now undefined; its value is not necessarily 9 or 11. 
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PROGRAMMING CONSICDERATIONS IN USING A DO LOOP 


1. 
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The indexing parameters of a DO statement (i, M4, Ma, m3) should 
not be changed by a statement within the range of the DO loop. 


There may be other DO statements within the range of a DO state- 
ment. All statements in the range of the inner DO must be in the 
range of the outer DO. A set of DO statements satisfying this rule 
is called a nest of DO's. 


Example 1: 


po 50 I=1, 4 


A(I) = B(I)**2 


DO 50 J=1, 5 Range of 
| Range of Outer DO 
50 C(J+1) = A(T) Inner DO 

Example 2: 

DO 10 INDEX = L, M 

N = INDEX + K 

DO 15 Jd = 1, 100, 2 Range of 

Range of Outer DO 

15 TABLE(J) = SUM(J,N)-1 Inner DO 


10 B(N) = A(N) 


A transfer out of the range of any DO loop is permissible at any 
time. 


The extended range of a DO is defined as those statements in the 
program unit containing the DO statement that are executed between 
the transfer out of the innermost DO of a nest of DO‘s and the 
transfer back into the range of this innermost DO. The following 
restrictions apply: 


e Transfer into the range of a DO is permitted only if such a 
transfer is from the extended range of the DO. 


e The extended range of a DO statement must not contain another DO 
statement that has an extended range if the second DO is within 
the same program unit as the first. 


e The indexing parameters (i, M1, M2, M3) cannot be changed in the 
extended range of the DO. 


Note that a statement that is the end of the range of more than one 
DO statement is within the innermost DO. The statement label of 
such a terminal statement may not be used in any GO TO or arithme- 
tic IF statement that occurs anywhere but in the range of the most 
deeply contained DO with that terminal statement. 


Example: 


— . 
>: : 


In the preceding example, the transfers specified by the numbers 1, 
2, and 3 are permissible, whereas those specified by 4, 5, 6, and 7 
are not. 


5. The indexing parameters (1,M,,Ma,mM3) may be changed by statements 
outside the range of the DO statement only if no transfer is made 
back into.the range of the DO statement that uses those parameters. 


6. The last statement in the range of a DO loop (statement x) must be 
an executable statement. It cannot be a GO TO statement of any 
form, Or a PAUSE, STOP, RETURN, arithmetic IF statement, another DO 
statement, or a logical IF statement containing any of these forms. 


7. ‘The use of, and return from, a subprogram from within any DO loop 
in a nest of DO's is permitted. 


CONTINUE STATEMENT 


Fe eg ae ee ee pag eG ee I gr Get RE Op gO TR eS AR gg ae Te ee 7 
] General Form | 
mwa 4 
| CONTINUE | 
Bi a a aa ca ee J 


CONTINUE is a dummy statement that may be placed anywhere in the 
source program without affecting the sequence of execution. It may be 
used as the last statement in the range of a DO in order to avoid ending 
the DO loop with a GO TO, PAUSE, STOP, RETURN, arithmetic IF, another DO 
Statement, or a logical IF statement containing any of these forms 
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Example 1: 


DO 30 I = 1, 20 
7 IF (A(I)-B(I)) 5, 30,30 
5 ACI) =ACI) +1 0 

B(I) = B(I) -2.0 


GO TO 7 
30 CONTINUE 
C = A(3) + R(7) 


Explanation: 


In example 1, the CONTINUE statement is used as the last statement in 
the range of the DO in order to avoid ending the DO loop with the state- 
ment GO TO 7. 


DO 30 I=1, 20 
IF(A(I)-B(1))5, 40, 40 
5 ACI) = C(T) 
GO TO 30 
40 ACI) = 0.0 
30 CONTINUE 


Explanation: 


In example 2, the CONTINUE statement provides a branch point enabling 
the programmer to bypass the execution of statement 40. 


PAUSE STATEMENT 


A AD AD A A ED ANS DRY ENS ART EY OT AEE NEY EE SRT GUNNS <emmD <OURES <ERD CREED Ce ime SORE eeline ED eee CEES OUND Seen SEE GRD ERD ERS Ce SED eee Se ey Ene SET ata ew eee CEES EY En CREE ce Stn REED ERY GRIDS GEER camer <EYED <Emer lw Wey EID cine Come Gn «iNT qe ates ques eT ones 


EEE CNR CORD coweD AES: GRP SEER <iemte omer <imeP qnome cin EE coe em ape <i We come coe eRe EURO oe Se Gee ET eRe em ce ee eee EP ces epee ee a Re ED GERD cD ne cece ED wen ee CET? ame ABe a ee Ge GEES ee cree SEP Ee cS ee cn RD ED eee eS ee oe ee ee eee 


PAUSE n 
PAUSE ‘message' 


Where: n is a string of 1 through 5 decimal digits. 


‘message’ is a literal constant of one form only: specific- 
ally, a string of alphameric and/or special characters 
enclosed in apostrophes. 


ome ee eee ee cmp eee ee ee ee ee ee ee ee ee ee ee eee a ce ee eee ee ee ee ee eee ee ee ee ee ene ee eee ee ee ee ee ee ee em ee ee ee ee ee ee ee ae ee ee ce ee ce ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee eee 


PAUSE n, PAUSE message, or PAUSE 00000 is displayed, depending upon 
whether n, ‘message’ or no parameter was specified, and the program 
waits until operator intervention causes it to resume execution, start- 
ing with the next statement after the PAUSE statement. For further 
information, see the FORTRAN programmers’ guides listed in the Preface. 


STOP STATEMENT 


Re a a a aa a aaa 1 
] General Form | 

Sea ad aN Ne en I eae ny Rr OM EP re NON ne ee PCE ONE PN re f 
| STOP | 
| STOP n | 
| | 
} Where: nis a string of 1 through 5 decimal digits. | 
eae a a ad Ee SP etn eR ener are rare Sede Pee Pe ec Co eae a ee eee J 


The STOP statement terminates the execution of the object program and 
displays n if specified. For further information, see the FORTRAN pro- 
grammers’ guides listed in the Preface. 


END STATEMENT 


The END statement is a nonexecutable statément that defines the end 
of a source program or source subprogram for the compiler. Physically, 
it must be the last statement of each program or subprogram. It may not 
have a statement number, and it may not be continued. The END statement 
does not terminate program execution. To terminate execution, a STOP 
statement Or a RETURN statement in the main program is required. 
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INPUT/OUTPUT STATEMENTS 


Input/output statements are used to transfer and control the flow of 
data between internal storage and an input/output device, such as a card 
reader, printer, punch, magnetic tape unit, or disk storage unit. The 
data that is to be transferred belongs to a data set. Data sets are 
composed of one or more records. Typical records are punched cards, 
printed lines, or the images of either on magnetic tape or disk. 


Operation: In order for the input or output operation to take place, 
the programmer must specify the kind of operation he desires; READ, 
WRITE, Or BACKSPACE, for example. 


its data set reference number. (The FORTRAN programmers' guides, listed 
in the preface, explain how data set reference numbers are associated 
with data sets.) In the statement specifying the type of input/output 
operation, the programmer must give the data set reference number corre- 
sponding to the data set he wishes to operate on. 


I/Q_ List: Input/output statements in FORTRAN are primarily concerned 
with the transfer of data between storage locations defined in a FORTRAN 
program and records which are external to the program. On input, data 
is taken from a record and placed into storage locations that are not 
necessarily contiguous. On output, data is gathered from diverse 
Storage locations and placed into a record. An I/O list is used to spe- 
cify which storage locations are used. The I/O list can contain vari- 
able names, subscripted array names, unsubscripted array names, or array 
names accompanied by indexing specifications in a form called an implied 
DO. No function references or arithmetic expressions are permitted in 
an I/O list. 


If a variable name or subscripted array name appears in the I/O list, 
one item is transmitted between a storage location and a record. 


If an unsubscripted array name appears in the list, the entire array 
is transmitted in the order in which it is stored. (If the array has 
more than one dimension, it is stored in ascending storage locations, 
with the value of the first subscript quantity increasing most rapidly 
and the value of the last increasing least rapidly. An example is given 
in the section “Arrangement of Arrays in Storage.") 


If an implied DO appears in the I/O list, the elements of the 
array(s) specified by the implied DO are transmitted. The implied DO 
specification is enclosed in parentheses Within the parentheses are 
one or more subscripted array names, separated by commas with a comma 
following the last name, followed by indexing parameters i=m4, May, Ma. 
The indexing parameters are as defined for the DO statement. Their 
range is the list of the DO-implied list and, for input lists, i, m4, 
Mo, and m3 may appear within that range only in subscripts. 
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For example, assume that A is a variable and that B, C, and D are 
1-dimensional arrays each containing 20 elements. Then the statement: 


WRITE (6) A, B, (C(I), I=1,4), D4) 


writes the current value of variable A, the entire array B, the first 
four elements of the array C, and the fourth element of D. (The 6 fol- 
lowing the WRITE is the data set reference number.) 


Implied DO's can be nested if required. For example, to read an ele- 
ment into array B after values are read into each row of a 10 x 20 array 
A, the following would be written: 


READ (5) ((A(I,J), J=1,10), BCI), I=1,20) 


The order of the names in the list specifies the order in which the 
data is transferred between the record and the storage locations. 


A special kind of I/O list called a NAMELIST list is explained in the 
section “READ and WRITE Using NAMELIST." 


Formatted and Unformatted Records: Data can be transmitted either under 


EE RD ED CS ET UNOS ERE EAP UNE i GAT TY GED ce SUE ete uni wee Com Woe snes MARI ES ERED GEES GEEREGEGRED WED GE cD eons 


control of a FORMAT statement or without the use of a FORMAT statement. 


When data is transmitted with format control, the data in the record 
is coded ina form that can be read by the programmer or satisfies the 
needs of machine representation. The transformation for input takes the 
character codes and constructs a machine representation of an item. The 
output transformation takes the machine representation of an item and 
constructs character codes suitable for printing. Most transformations 
involve numeric representations that require base conversion. To obtain 
format control, the programmer must include a FORMAT statement in the 
program and must give the statement number of the FORMAT statement in 
the READ or WRITE statement specifying the input/output operation. 


When data is transmitted without format control, no FORMAT statement 
is used. In this case, there is a one-to-one correspondence between 
internal storage locations (bytes) and external record positions. A 
typical use of unformatted data is for information that is written out 
during a program, not examined by the programmer, and then read back in 
later in the program or in another program for additional processing. 


For unformatted data, the I/O list determines the length of the rec- 
ord. For example, an output record is complete when the current values 
of all the items in the I/O list have been placed in it, plus any con- 
trol words supplied by the input/output routines or Data Management. 
For further information, see the FORTRAN IV programmers’ guides listed 
in the Preface 


For formatted data, the I/O list and the FORMAT statement determine 
the form of the record. For further information see the section "FORMAT 
Statement™ and the FORTRAN IV programmers’ guides 
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There are two types of input/output statements: sequential and 
direct access Sequential input/output statements are used for storing 
and retrieving data sequentially. These statements are device indepen- 
dent and can be used for data sets on either sequential or direct access 
devices. 


The direct access input/output statements are used to store and 
retrieve data in an order specified by the user. These statements can 
be used only for a data set on a direct access storage device. 


SEQUENTIAL INPUT/CUTPUT STATEMENTS 


There are five sequential input/output statements: READ, WRITE, END 
FILE, REWIND, and BACKSPACE. The READ and WRITE statements cause 
transfer of records of sequential data sets. The END FILE statement 
defines the end of a data set; the REWIND and BACKSPACE statements con- 
trol the positioning of data sets. In addition to these five state- 
ments, the FORMAT and NAMELIST statements, although not input/ 
output statements, are used with certain forms of the READ and WRITE 
statements. 


READ STATEMENT 
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| PEAD(a,b, END=c, ERR=d) list 
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a is an unsigned integer constant or an integer variable 
that is of length 4 and represents a data set reference 
number. 


b is optional and is either the statement number or array 
name of the FORMAT statement describing the record(s) being 
read, or a NAMELIST name. 


END=c is optional and ¢c is the number of the statement to 
which transfer is made upon encountering the end of the data 
set. 


ERR=d is optional and d is the number of the statement to 


which transfer is made upon encountering an error condition 
in data transfer. 


list is optional and is an I/O list. 
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The READ statement may take many forms. The value of a must always 
be specified, kut under appropriate conditions b, c, d, and list can be 
omitted. The order of the parameters END=c and ERR=d can be reversed 
within the parentheses. 


ay 


Transfer is made to the statement specified by the END parameter when 
the end of the data set is encountered; i e., when a READ statement is 
executed after the last record on the data set has already been read. 
(No indication is given of the number of list items read into before the 
end of the data set was encountered.) If the END parameter is omitted, 
object program execution is terminated upon encountering the end of the 
data set. 


Transfer is made to the statement specified by the ERR parameter if 
an input/output device error occurs. No data is read into the list 
items and no indication is given of which record or records could not be 
read, only that an error occurred during transmission of data. If the 
ERR parameter is omitted, object program execution is terminated when an 
input/output device error occurs. 


The basic forms of the READ statements are: 








Form Purpose 

READ(a,b) list Formatted READ 
READ(a) list Unformatted READ 
READ (a, x) READ using NAMELIST 


The discussion of READ using NAMELIST is in the section "READ and 
WRITE Using NAMELIST." 


Formatted READ 


The form READ (a,b) list is used to read data from the data set asso- 
ciated with data set reference number a into the variables whose names 
are given in the list. The data is transmitted from the data set to 
storage according to the specifications in the FORMAT statement, which 
is statement number b. 


Example: 
READ (5,98) A,B, (C(1I,K),1I=1,10) 


Explanation: The above statement causes input data to be read from the 
data set associated with data set reference number 5 into the variables 
A, EB, Cl1,K), C(2,K),.<e<, C(10,K) in the format specified by the FORMAT 
statement whose statement number is 98. 


Unformatted READ 


The form READ(a) list is used to read a single record from the data 
set associated with data set reference number a into the variables whose 
names are given in the list. Since the data is unformatted, no FORMAT 
statement number is given. This statement is used to read unformatted 
data written by a WRITE(a) list statement. If the list is omitted, a 


record is passed over without being processed. 
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Example: 
READ (J) A,B,C 


Explanation: The above statement causes data to be read from the data 


set associated with data set reference number J into the variables A, B, 
and C. 


WRITE STATEMENT 


|General Form 

]WRITE(a,b) list 

jWhere: a is an unsigned integer constant or an integer variable 
that is of length 4 and represents a data set reference 
number. 


b is optional and is either the statement number or array 
name of the FORMAT statement describing the record(s) being 
written, or a NAMELIST name. 


list is optional and is an I/O list. 
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The WRITE statement may take many different forms. For example, the 
list or the parameter b may be omitted. 


The three basic forms of the WRITE statement are: 


Form Purpose 

WRITE(a,b) list Formatted WRITE 
WRITE(a) list Unformatted WRITE 
WRITE(a, x) WRITE using NAMELIST 


The discussion of WRITE using NAMELIST is in the section "READ and 
WRITE Using NAMELIST." 


Formatted WRITE 


The form WRITE(a,b) list is used to write data into the data set 
whose reference number is a from the variables whose names are given in 
the list. The data is transmitted from storage to the data set accord- 
ing to the specifications in the FORMAT statement, whose statement num- 


ber is b. 
Example: 
WRITE(7, 75)A, (B(I, 3), 1=1,10,2),C 
Explanation: The above statement causes data to be written from the 
variables A, B(1,3), B(3,3), B(5,3), B(7,3), B(9,3), C into the data set 


associated with data set reference number 7 in the format specified by 
the FORMAT statement whose statement number is 75. 
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Unformatted WRITE 


The form WRITE(a) list is used to write a single record from the 
variables whose names are given in the list into the data set whose data 
set reference number is a. This data can be read back into storage with 
the unformatted form of the READ statement, READ(a) list. The list can- 
not be omitted. 


Example: _ 


WRITE (L) ((A(I,J),1I=1,10,2), B(J,3), J=1,K) 


Explanation: The above statement causes data to be written from the 
variables A(1,1), A(3,1),..., A(9,1), B(1,3), A(1,2), A(3,2),..., 
A(9,2), B(2,3),..., B(K,3) into the data set associated with the data 
set reference number L. Since the record is unformatted, no FORMAT 
statement number is given. Therefore, no FORMAT statement number should 
be given in the READ statement used to read the data back into storage. 


READ AND WRITE USING NAMELIST 


The NAMELIST statement is used in conjunction with the READ(a,x) and 
WRITE(a,x) statements to provide for reading and writing data without 
including the list specification in the READ and WRITE statements. The 
NAMELIST statement declares a name x to refer to a particular list of 
variables or array names. Neither a dummy variable name nor a dummy 
array name may appear in the list. Thereafter, the forms READ(a,x) and 
WRITE(a,x) are used to transmit data between the data set associated 
with the reference number a and the variables specified by the NAMELIST 
name x. 


The format and rules for constructing and using the NAMELIST state- 
ments are described in the following text. 


IGeneral Form 
[NAMELIST/x/a,be.-C/y/dyepee-f/2/Gbyeeed 
hers: XeVYe and Z,--. are NAMELIST names. ! 
| ayC,yC,d,.e-. are variable or array names. | 
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The following rules apply to declaring and using a NAMELIST name: 
1. A NAMELIST name is a symbolic name. 
2. A NAMELIST name is enclosed in slashes. The list of variable or 


array names Ltelonging to a NAMELIST name ends with a new NAMELIST 
name enclosed in slashes or with the end of the NAMELIST statement. 
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3. A variable name or an array name may belong to one or more NAMELIST 
lists. 


4, A NAMELIST name must be declared in a NAMELIST statement before it 
is used in an input/output statement, and it may be declared only 
once. After it is declared, it may appear only in input/output 
statements. 


5. The rules for input/output conversion of NAMELIST data are the same 
as the rules for data conversion described in the section "FORMAT 
Statement." The NAMELIST data must be in a special form, described 
in the following sections. 


6. A NAMELIST name may not be used as an argument. 


NAMELIST Input _Data 


Input data must be in a special form in order to be read uSing a 
NAMELIST list. The first character in each record to be read must be 
blank. The second character in the first record of a group of data 
records must be an & immediately followed by the NAMELIST name. The 
NAMELIST name must be followed by a blank and must not contain any 
embedded blanks. This name is followed by data items separated by com- 
mas. (A comma after the last item is optional.) The end of a data 
group is signaled by &END. 


The form of the data items in an input record may be: 
e variable name = constant 


The variable name may be a subscripted array name or a single vari- 
able name. Subscripts must be integer constants. The constant may 
be integer, real, literal, complex, or logical. (If the constants 

are logical, they may be in the form T or .TRUE. and F or .FALSE.) 





® array name = set of constants (separated by commas) 


The array name is not subscripted. The set of constants consists of 
constants of the type integer, real, literal, complex, or logical. 
The number of constants must be less than or equal to the number of 
elements in the array. Successive occurrences of the same constant 
can be represented in the form k*constant. 


The variable names and array names specified in the input data set 
must appear in the NAMELIST list, but the order is not significant. A 
name that has been made equivalent to a name in the input data cannot be 
substituted for that name in the NAMELIST list. The list can contain 
names of items in COMMON but must not contain dummy argument names. 


Each data record must begin with a complete variable or array name or 
constant. Embedded blanks are not permitted in names or constants. 
Trailing blanks after integers and exponents are treated as zeros. 
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NAMELIST Output Data 


When output data is written using a NAMELIST list, it is written ina 
form that can be read using a NAMELIST list. All variable and array 
names specified in the NAMELIST list and their values are written out, 
each according to its type. The fields for the data are made large 
enough to contain all the significant digits. The values of a complete 
array are written out in columns. 


Example: Assume that A is a 3 by 1 array, I and L are 3 by 3 arrays, 
and that the following statements are given: 


NAMELIST /NAM1/A,B, 1,3, L/NAM2/C,J,1,L 
READ (5, NAM1) 
WRITE (6, NAM2) 


Explanation: The NAMELIST statement defines two NAMELIST lists, NAM1 
and NAM2. The READ statement causes input data to be read from the data 
set associated with data set reference number 5 into the variables and 


arrays specified by NAM1. Assume that the data cards have the form: 


Column 2 

| 

Vv 
First card ENAM1 1(2, 3) =5, J=4, B=3.2 
Last card A(3)=4.0,L=2, 3, 7*4, &END 


The first data card is read and examined to verify that its name is con- 
Sistent with the NAMELIST name in the READ statement. (If that NAMELIST 
name is not found, then it reads to the next NAMELIST group.) When the 
data is read, the integer constants 5 and 4 are placed in I(2,3) and Jd, 
respectively; and the real constants 3.2 and 4.0 are placed in B and 
A(3), respectively. Since L is an array name not followed by a sub- 
script, the entire array is filled with the succeeding constants. 
Therefore, the integer constants 2 and 3 are placed in L(1,1) and 
L(2,1), respectively, and the integer constant 4 is placed in L(3,1), 
Ll, 2) p0:80¢° C353) 


The WRITE statement causes data to be written from the variables and 
arrays specified by NAM2 into the data set associated with data set 
reference number 6. Assume that the values of J, L, and I(2,3) were not 
altered since the previous READ statement, that C was given the value 
428.0E+03, that I(1,3) was given the value 6, and that the rest of the 
elements of I were set to zero. Then, if the output is punched on 
cards, the form is: 


Column 2 
Vv 
First card §NAM2 
Second card C=428000.00, J=4, 1=0,0,0,0,0,0,6,5, 
Third card 0,L=2,3,4,4,4,4,4,4,4, 
Fourth card SEND 
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FORMAT STATEMENT 


Sh ees iN a a ea aa a ak a Na a 1 
}General Form | 
}----------------------------~----------------------------------------- { 
| XXXXX FORMAT (C4, CazeceygCn) | 
l | 
|Where: XXXXX is a statement number (1 through 5 digits). | 
| | 
] CaseCageeeeCyn are format codes. 
| 
|The format codes are: | 
| 
] alw (Descrikes integer data fields.) 
| paDw.d (Describes real data fields.) | 
| paEw.d (Describes real data fields.) 
| paFw.d (Describes real data fields.) | 
| aZw (Describes hexadecimal data fields.) | 
|] paGw.s (Describes integer, real, complex, Or logical data fields.) | 
| aLw (Describes logical data fields.) 
] aAw (Descrikes alphameric data fields.) | 
|*Literal® (Transmits literal data.) | 
] wH (Transmits literal data.) | 
] wX (Indicates that a field is to be skipped on input or filled | 
{ with blanks on output.) | 
| Tx (Indicates the position in a FORTRAN record where transfer | 
] of data is to start.) | 
Jat...) (Indicates a group format specification.) | 
| 
|Where a is optional and is an unsigned integer constant used to | 
] denote the number of times the format code is to be used. | 
| If a is omitted, the code is used only once. | 
| | 
| w is an unsigned nonzero integer constant that specifies the| 
] number of characters in the field. | 
| | | 
| d is an unsigned integer constant specifying the number of | 
| decimal places to the right of the decimal point; i.e., the | 
] fractional portion. | 
| | 
| s is an unsigned integer constant specifying the number of | 
| significant digits. | 
| | 
| ris an unsigned integer constant designating a character | 
] position in a record. | 
| | 
| p is optional and represents a scale factor designator of | 
| the form nP where n is an unsigned or negatively signed | 
| integer constant. | 
| | | | 
| (...-) is a group format specification. Within the paren- | 
] theses are format codes separated by commas or slashes. | 
| Group format specifications can be nested to a level of two. ] 
] The a preceding this form is called a group repeat count. | 
}--------------=------------------------------------------------------- { 
| Notes: | 
{ 1. Complex data fields in records require two successive D, hy 2; | 
| G, or A format codes. These codes may be grouped within | 
| parentheses | 
| 
} 2 Roth commas and slashes can be used as separators between format | 
| codes (see the section “Various Forms of a FORMAT Statement"). | 
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The FORMAT statement is used in conjunction with the I/O list in the 
READ and WRITE statements to specify the structure of FORTRAN records 
and the form of the data fields within the records. In the FORMAT 
statement, the data fields are described with format codes, and the 
order in which these format codes are specified gives the structure of 
the FORTRAN records. The I/O list gives the names of the data items to 
make up the record. The length of the list in conjunction with the FOR- 
MAT statement specifies the length of the record (see the section 
"Various Forms of a FORMAT Statement"). Throughout this section, the 
examples show punched card input and printed line output. The concepts 
apply to all input/output media. In the examples, the character b 
represents a blank. 


The following list gives general rules for using the FORMAT 
Statement: 


1. FORMAT statements are not executed; their function is to supply 
information to the object program. They may be placed anywhere in 
the source program. 


2. When defining a FORTRAN record by a FORMAT statement, it is impor- 
tant to consider the maximum size record allowed on the input/ 
output medium. For example, if a FORTRAN record is to be punched 
for output, the record should not be longer than 80 characters. If 
it is to be printed, it should not be longer than the printer's 
line length. For input, the FORMAT statement should not define a 
FORTRAN record longer than the record referred to in the data set. 


3. When formatted records are prepared for printing, the first 
character of the record is not printed. It is treated as a car- 
riage control character. It can be specified in a FORMAT statement 
with either of two forms of literal data: either ‘'x‘' or 1Hx, where 
x is one of the following: 


x Meaning 

blank Advance one line before printing 

0 Advance two lines before printing 
1 Advance to first line of next page 
+ No advance 


For media other than the printer, the first character of the record 
is treated as data. 


4. If the I/O list is omitted from the READ or WRITE statement, a 
record is skipped on input, or a blank record is inserted on out- 
put, unless the record was transmitted between the data set and the 
FORMAT statement (see "H Format Code and Literal Data"). 


Various Forms of a FORMAT Statement 


All of the format codes in a FORMAT statement are enclosed in a pair 
of parentheses. Within these parentheses, the format codes are delin- 
ited by the separators: comma and slash. 


Execution of a formatted READ or formatted WRITE statement initiates 
format control Fach action of format control depends on information _ 
provided jointly Ly the I/O list, if one exists, and the format specifi- 
cation. There is no I/O list item corresponding to the format descrip- 
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tors X, H, and literals enclosed in apostrophes. These communicate 
information directly with the record. 


Whenever an I, D, Ey, Fy, G, A, Ly, Or Z code is encountered, format 
control determines whether there is a corresponding element in the I/O 
list. If there is such an element, appropriately converted information 
is transmitted. If there is no corresponding element, the format con- 
trol terminates. 


If, however, format control reaches the last outer right parenthesis 
of the format specification and another element is specified in the I/O 
list, control is transferred to the group repeat count of the group for- 
mat specification terminated by the last right parenthesis that precedes 
the right parenthesis ending the FORMAT statement. 


The question of whether there are further elements in the I/O list is 
asked only when an I, D, E, F, G, A, L, or Z code or the final right 
parenthesis of the format specification is encountered. Before this is 
done, T, X, and H codes, literals enclosed in apostrophes, and slashes 
are processed. If there are fewer elements in the I/O list than there 
are format codes, the remaining format codes are ignored. 


Comma: The simplest form of a FORMAT statement is the one shown in the 
box at the beginning of this section with the format codes, separated by 
commas, enclosed in a pair of parentheses. One FORTRAN record is 
defined by the beginning of the FORMAT statement (left parenthesis) to 
the end of the FORMAT statement (right parenthesis). For an example, 
see the section “Examples of Numeric Format Codes." 


Slash: A slash is used to indicate the end of a FORTRAN record format. 
For example, the statement: 


25 FORMAT (13, F6.2/D10. 3, F6. 2) 


describes two FORTRAN record formats. The first, third, etc., records 
are transmitted according to the format I3, F6.2 and the second, fourth, 
etc., records are transmitted according to the format D10.3, F6.2. 


Consecutive slashes can be used to introduce blank output records or 
to skip input records. If there are n consecutive slashes at the begin- 
ning or end of a FORMAT statement, n input records are skipped or n 
blank records are inserted between output records. If n consecutive 
Slashes appear anywhere else in a FORMAT statement, the number of 
records skipped or blank records inserted is n-1. For example, the 
statement: 


25 FORMAT (1X,10I5//1X, 8E14.5) 


describes three FORTRAN record formats. On output, it causes double 
Spacing between the line written with format 1X,101I5 and the line writ- 
ten with the format 1X,8E14.5. 


I Format Code 


The I format code is used in transmitting integer data. For example, 
if a READ statement refers to a FORMAT statement containing I format 
codes, the input data is stored in internal storage in integer format. 
The magnitude of the data to be transmitted must not exceed the maximum 
magnitude of an integer constant. 
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Input: Leading, embedded, and trailing blanks in a field of the input 
card are interpreted as zeros. 


Output: If the number of significant digits and sign required to repre- 
sent the quantity in the storage location is less than w, the leftmost 
print positions are filled with blanks. If it is greater than w, 
asterisks are printed instead of the number. 


D, E, and F_ Format Codes 


The D, E, and F format codes are used in transmitting real or double 
precision data. The data must not exceed the maximum magnitude for a 
real or double precision constant. 


Input: Input must be a real or double precision number which, optional- 
ly, may have a D or E exponent. The decimal point may be omitted. If 
it is present, its position overrides the position indicated by the d 
portion of the format field descriptor, and the number of positions 
specified by w must include a place for it. If the data has a D or E 
exponent and the format field descriptor includes a P scale factor, the 
scale factor has no effect. Each data item must be right justified in 
its field, since leading, trailing, and embedded blanks are treated as 
zeros. These three format codes are interchangeable for input. It 
makes no difference, for example, whether D, E, or F is used to describe 
a field containing 12.42E+08. 


Output: For data written under a D or E format code, unless a P scale 
factor is specified, output consists of an optional sign (required for 
negative values), a decimal point, the number of significant digits 
specified by d, and a D or E exponent requiring four positions. The w 
specification must provide for all these positions, including the one 
for a sign when the output value is negative. If additional space is 


available, a leading zero may be written before the decimal point. 


For data written under an F format code, w must provide sufficient 
spaces for an integer segment if it is other than zero, a fractional 
segment containing d digits, a decimal point, and, if the output value 
is negative, a sign. If insufficient positions are provided for the 
integer portion, including the decimal point and sign (if any), 
asterisks are written instead of data. If excess positions are pro- 
vided, the number is preceded by blanks. 


For D, E, and F, fractional digits in excess of the number specified 
by d are dropped after rounding. 


Z Format Code 


The Z format code is used in transmitting hexadecimal data. 


Input: Leading, embedded, and trailing blanks in an input field are 
treated as zeros. One storage location (byte) in internal storage con- 
tains two hexadecimal digits; thus, if an input field contains an odd 
number of digits, the number will be padded on the left with a hexa- 
decimal zero when it is stored. 


Output: If the number of characters in the storage location is less 
than w, the leftmost print positions are filled with blanks If the 
number of characters in the storage location is greater than w, the 
leftmost digits are truncated and the rest of the number is printed. 
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G Format Code 


The G format code is a generalized code used to transmit integer, 
real, complex, or logical data according to the type ere ne orn of 
the corresponding variable in the I/O list. 


Input: The rules for input for G format code depend upon the type of 
the variable in the I/O list and the form of the number punched on the 
card. For example, if the variable is real and the number punched in 
the card has an E decimal exponent, the rules are the same as for the E 
format code. If the variable in the I/O list is integer or logical, the 
Ss portion of the format code, specifying the number of significant 
digits, can be omitted; if it is given, it is ignored. For complex and 
real data, the s portion gives the location of the implied decimal point 
for input -- just like the d specification for D, E, and F format codes. 


Output: The s portion of the format code can be omitted for integer and 
logical data and the numbers are printed according to the rules for I 
and L format codes. For complex and real data, the s is used to deter- 
mine the number of digits to be printed and whether the number should be 
printed with or without a decimal exponent. If the number, say n, is in 
the range 0.1< n < 10**s, the number is printed without a decimal 
exponent. Otherwise, it is printed with an E or D decimal exponent 
depending on the length specification of the variable in the I/O list. 
The w specification for complex and real data must include a position 
for a decimal point, four positions for a decimal exponent, and, if the 
value is negative, a position for a minus sign. All other rules for 
output are the same as those for the individual format codes. 


Examples of Numeric Format Codes 


The following examples illustrate the use of the format codes I, F, 
D, E, Z, and G. 


Example 1: 
75 FORMAT (13,F5.2,E10.3,G10. 3) 


READ (5,75) N,A,B,C 


FR re 


1. Four input fields are described in the FORMAT statement and four 
variables are in the I/O list. Therefore, each time the READ 
statement is executed, one input card is read from the data set 
associated with data set reference number 5. 


2. When an input card is read, the number in the first field of the 
card (three columns) is stored in integer format in location N, 
The number in the second field of the input card (five columns) is 
stored in real format, with no decimal exponent, in location A, 
etc. 


3. If there were one more variable in the I/O list, say M, another 
card would be read and the information in the first three columns 
in that card would be stored in integer format in location M. The 
rest of the data on the card would be ignored. 


4. If there were one fewer variable in the list (say C is omitted), no 
number would be stored according to the format G10.3. 
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5 This format statement defines only one record format. The section 
“Various Forms Of a FORMAT Statement" explains how to define more 
than one record format in a FORMAT statement 

Example 2: Assume that the following statements are given: 
75 FORMAT (Z4,D10.3, 2G10. 3) 
READ (5,75) A,B,C,D 


where A, C, and D are REAL*¥4 and B is REAL*8 and that the following 
input cards are read: 


Column 1 5 15 25 35 


] | | | 
Vv V V Vv ov 
b3F1156432D+02276. 38E+15bbbbbbbbbb 

Input 

ZAF 3155381+4+02b382506E+28276,. 38E+15 
Cards ] 

3ACbH346.18D-03485. 322836276. 38E+15 
Format Z4 D10.3 G10.3 G10.3 


Then the variakles A, B, C, and D receive values as if the following 
had been punched: 


A B ic D 

O3F1 156.432D+02 276. 38E+15 000000. 000 
2AF3 155. 381D+020 382. 506E+28 276. 38E+15 
3Ac0 346.18D-03 485. 322836 276. 38E+15 


Explanation: 


1. Leading, trailing, and embedded blanks in an input field are 
treated as zeros. Therefore, since the value for B on the second 
input card was not right justified in the field, the exponent is 20 
not 2. 


2. Values read into the variables C and D with a G format code are 


converted according to the type of the corresponding variable in 
the I/O list. 


Example 3: Assume that the following statements are given: 
76 FORMAT ('0°,F6.2,812.3,G614.6,15) 
WRITE (6, 76)A,B,C,N 


and that the variakles A, B, C, and N have the following values: 


A B c N 
034.40 123. 380E+02 123. 380E+02 031 
031.1 1156.1E+02 123456789. 130 
-354 32 834 621E-03 1234 56789 428 
01.132 83,121E+06 123380 D+02 000 
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Then, the following lines are printed: 


Column 2 3 
1 9 1 5 
34.40 0 123E 05 12338.0 31 


31.10 0.116EF 06 O 123457E 09 130 
KK 0.835E 00 1234.57 428 


1.13 0.831E 08 0.123380E 08 0 


Explanation: 


1. The integer portion of the third value of A exceeds the format spe- 
cification, so asterisks are printed instead of a value. The frac- 
tional portion of the fourth value of A exceeds the format specifi- 
cation, so the fractional portion is rounded. 


2. Note that for the variable B the decimal point is printed to the 
left of the first significant digit and that only three significant 
digits are printed because of the format specification E12. 3. 
Excess digits are rounded off from the right. 


3. The values of the variable C are printed according to the format 
specification G14.6. The s specification, which in this case is 6, 
determines the number of digits to be printed and whether the 
number should be printed with a decimal exponent. Values greater 
than or equal to 0.1 and less than 1,000,000 are printed without a 
decimal exponent in this example. Thus, the first and third values 
have no exponent. The second and fourth values are greater than 
1,000,000, so they are printed with an exponent. 


Scale Factor -_P 


The P scale factor may be specified as the first part of aD, EH, F, 
or G field descriptor to change the location of the decimal point in 
real numbers. The effect of the scale factor is: 


scale factor 
external numker = internal number x 10 


Input: A scale factor may be specified for any real data, but it is 
ignored for any data item that contains an exponent in the external 
field. For example, if the input data is in the form xx.xxxx and is to 
be used internally in the form .~.xxxxxx, then the format code used to 
effect this change is 2PF7.4. Or, if the input data is in the form xx. 
xXxXxXx and is to be used internally in the form xxxx.xx, then the format 
code used to effect this change is -2PF7.4. 


Output: A scale factor can be specified for real numbers with or 
without E or D decimal exponents. For numbers without an E or D decimal 
exponent, the effect is the same as for input data except that the 
decimal point is moved in the opposite direction For example, if the 
number has the internal form xx.xxxx and is to be written out in the 


form xxxx.xx, the format code used to effect this change is 2PF7.4. 


56 


For numbers with an E or D decimal exponent, when the decimal point 
is moved, the exponent is adjusted to account for it, i.ee., the value is 
not changed. For example, if the internal number 238. were printed 
according to the format E10.3, it would appear as 0 238Eb03. If it were 
printed according to the format 1PE10.3, it would appear as 2.380Eb02. 


A repetition code can precede the D, E, or F format code. For 
example, 2P3F7.4 is valid. 


Warning: Once a scale factor has been established, it applies to all 
subsequently interpreted D, E, F, and G codes in the same FORMAT state- 
ment until another scale factor is encountered. The new scale factor is 
then established. A factor of 0 may be used to discontinue the effect 


of a previous scale factor. 


L_Format Code 


The L format code is used in transmitting logical variables. 


Input: The first T or F encountered in the w characters of the input 
field causes a value of .TRUE. or .FALSE., respectively, to be assigned 
to the corresponding logical variable in the I/O list. If the field w 
consists entirely of blanks, a value of .FALSE. is assumed. 


Qutput: AT or F is inserted in the output record depending upon wheth- 
er the value of the logical variable in the I/O list was .TRUE. or 
-FALSE., respectively. The single character is right justified in the 
output field and preceded by w-1 blanks. 


A Format Code 


The A format code is used in transmitting data that is stored 
internally in character format. The number of characters transmitted 
under A format code depends on the length of the corresponding variable 
in the I/O list. Each alphabetic or special character is given a unique 
internal code. Numeric data is converted digit by digit into internal 
format, rather than the entire numeric field being converted into a 
Single binary number. Thus, the A format code can be used for numeric 
fields, but not for numeric fields requiring arithmetic. 


Input: The maximum number of characters stored in internal storage 
depends on the length of the variable in the I/O list. If w is greater 
than the variakle length, say v, then the leftmost w-v characters in the 
field of the input card are skipped and the remaining v characters are 
read and stored in the variable. If wis less than v, then w characters 
from the field in the input card are read and the remaining rightmost 
characters in the variable are filled with blanks. 


Cutput: If w is greater than the length of the variable in the I/O 
list, say the length is v, then the printed field will contain v charac- 
ters right-justified in the field, preceded by leading blanks. If w is 
less than v, the leftmost w characters from the variable will be printed 
and the rest of the data will be truncated. | 


Example 1: Assume that B has been specified as real of length 8, that N 


and M are integers of standard length 4, and that the following state- 
rents are given. 
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25 FORMAT (3A7) 
READ (5,25) B, N, M 


When the READ statement is executed, one input card is read from the 
data set associated with data set reference number 5 into the variables 
B, N, and M in the format specified by FORMAT statement number 25. The 
following list shows the values stored for the given input cards (b 
represents a blank). 


Input_Card B N M 
ABCDEFGHU6bLATH11234567 ABCDEFGhb ATb1 4567 


HIJKLMN7 6543213334445 HIJKLMNb 4321 GUUS 


Example 2: Assume that A and B are real variables of length 4, that C 
is a real variable of length 8, and that the following statements are 
given: | | 


26 FORMAT (A6,A5,A6) 
WRITE (6,26) A,B,C 


When the WRITE statement is executed, one line is written on the data 
set associated with data set reference number 6 from the variables A, B, 
and C in the format specified by FORMAT statement 26. The following 
list shows the printed output for values of A, B, and C (b represents a 
blank). | 


A RB C Printed Line 


ane ~ome Et AR TRS GENE ED EY VOROP ite EE wes SEED enn 


A1B2 C3D4 E5F6G7H8 bbA1B2bC 3D4E5F6G7 
BH Format Code ang titeret Data 


Literal data can appear in a FORMAT statement in one of two ways: it 
can be enclosed in apostrophes or it can follow the H format code. For 
example, the following FORMAT statements are equivalent: 


25 FORMAT (* 1968 INVENTORY REPORT') 
25 FORMAT (22H 1968 INVENTORY REPORT) 


No item in the I/C list corresponds to the literal data. The data is 
read or written directly into or from the FORMAT statement. (The FORMAT 
Statement can contain other types of format codes with corresponding 
variables in the I/O list.) | 


Input: Information is read from the input card and replaces the literal 
data in the FORMAT statement. (If the H format code is used, w charac- 
ters are read. If apostrophes are used, as many characters as there are 
spaces between the apostrophes are read.) For example, the following 
statements: | 


8 FORMAT (' HEADINGS‘) 
READ (5,8) 
cause the first nine characters of the next record to be read from the 
data set associated with data set reference number 5 into the FORMAT 


Statement 8, replacing the blank and the eight characters H, E, A, D, I, 
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Qutput: The literal data from the FORMAT statement is written on the 
output data set (If the H format code is used, the w characters fol- 
lowing the H are written. If apostrophes are used, the characters 
enclosed in apostrophes are written.) For example, the following 
statements: 


8 FORMAT (14HOMEAN AVERAGE:, F8.4) 
WRITE (6,8) AVRGE 


would cause the following record to be written if the value of AVRGE 
were 12.3456: 


MEAN AVERAGE: 12.3456 
Note: If the literal data is enclosed in apostrophes, an apostrophe 


character in the data is represented by two successive apostrophes. For 
example, DON'T is represented as DON'T. 


X_Format Code 


The X format code specifies a field of w characters to be skipped on 
input or filled with blanks on output. For example, the following 
statements: 

5 FORMAT (110,10X, 4110) 
READ (5,5) I,J,K,L,M 
cause the first ten characters of the input card to be read into vari- 
able I, the next ten characters to be skipped over without transmission, 


and the next four fields of ten characters each to be read into the 
variables J, K, L, and M. 


T Format Code 


The T format code specifies the position in the FORTRAN record where 
the transfer of data is to begin. (Note that for printed output, the 
first character of the output data record is used for carriage control 
and is not printed. Thus, for example, if T50, ‘Z* is specified ina 
FORMAT statement, a Z will be the 50th character of the output record, 
but it will appear in the 49th print position.) 

The following illustrates the use of the T code. 

5 FORMAT (T40,°1968 STATISTICAL REPORT!,T80, 
X DECEMBER',T1,"OPART NO. 10095") 
WRITE (6,5) 


cause the following line to be printed: 


Print Print Print 
Position 1 Position 39 Position 79 
| | | 

Vv Vv Vv 

PART NO. 10095 1968 STATISTICAL REPORT DECEMBER 


The T format code can be used in a FORMAT statement with any type of 
format code, as, for example, with FORMAT ('0',T40,15). 
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Group Format Specification 


The group format specification is used to repeat a set of format 
codes and to control the order in which the format codes are used. 


The group repeat count a is the same as the repeat indicator a which 


can be placed in front of other format codes. For example, the follow- 
ing statements are equivalent: 


10 FORMAT (13,2(14,15),16) 


10 FORMAT (13, (14,15,14,15),16) 


Group repeat specifications control the order in which format codes 
are used since control returns to the last group repeat specification 
when there are more items in the I/O list than there are format codes in 
the FORMAT statement (see "Various Forms of a FORMAT Statement"). Thus, 
in the previous example, if there were more than six items in the I/O 
list, control would return to the group repeat count 2 which precedes 
the specification {(1I4,15). 


If the group repeat count is omitted, a count of 1 is assumed. For 
example, the statements: 


15 FORMAT (13, (F6.2,D10.3)) 
READ (5,15) N,A,B,C,D,E 


cause values to be read from the first record for N, A, and B, according 
to the format codes 13,F6.2, and D10.3, respectively. Then, because the 
I/O list is not exhausted, control returns to the last group repeat spe- 
cification, the next record is read, and values are transmitted to C and 
D according to the format codes F6.2 and D10.3, respectively. Since the 
I/O list is still not exhausted, another record is read and a value is 
transmitted to E according to the format code F6.2 -- the format code 
bDi0.3 is not used. 


The format codes within the group repeat specification can be 
separated by commas and slashes. For example, the following statement 
is valid: 


The first record is transmitted according to the specification 213, the 
second, fourth, etc., records are transmitted according to the specifi- 
cation 3F6.2,F6.3, and the third, fifth, etc., records are transmitted 
according to the specification D10.3,3D10.2, until the I/O list is 
exhausted. 


Reading Format Specifications at _Cbject Time 


FORTRAN provides for variable format statements by allowing a format 
specification to ke read into an array in storage. The data in the 
array may then be used as the format specification for subsequent input/ 
output operations. The format specification may also be placed into the 
array by a DATA statement or an explicit specification statement in the 
source program. 
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1. The name of the array containing the variable FORMAT specification 
must appear in a DIMENSION, COMMON, or explicit specification 
statement, even if the array size is only 1. 

2e The format codes entered into the array must have the same form as 
a source program FORMAT statement, except that the word FORMAT and 
the statement number are omitted. 

3. If a format code read in at object time contains double apostrophes 
within a literal field that is defined by apostrophes, it should be 
used for output only. If an object time format code is to be used 
for input, and if it must contain a literal field with an internal 
apostrophe, the H format code must be used for the literal field 
definition. 

Example: Assume that the following statements are given: 
DIMENSION FMT (18) 
1 FORMAT (18A4) 
READ (5,1) EMT 
READ (5, FMT) A,B, (C(I),I=1,5) 


and that the first input card associated with data set reference number 
5 contains (2F10.3, 5F10.8). 


The data on the rest of the input cards is read, converted, and stored 
in A, B, and the array C, according to the format codes 2E10.3, 5F10.8. 


END FILE STATEMENT 


| | 
| | | 
| Where: ais an unsigned integer constant or integer variable that 
] is of length 4 and represents a data set reference number. ] 


The END FILE statement defines the end of the data set associated 
with a. 


REWIND STATEMENT 


] Where: ais an unsigned integer constant or integer variable that 
| is of length 4 and represents a data set reference number. 


The REWIND statement causes a subsequent READ or WRITE statement ref- 
erring to a tc read data from or write data into the first record of the 
data set associated with a. 
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BACKSPACE STATEMENT 


Se ce ist cues a Ss ps Sk ih i ni in gp Si ee iin <p“ ssn a i ss a ew sl ee Se ea Sc ei hie ca i es a es aa es ge a ei its ss ae ce 4 
| BACKSPACE a | 
| 
] Where: ais an unsigned integer constant or integer variable that | 
] is cf length 4 and represents a data set reference number. | 
becenes eee cae eee Re NS ero ne ee ele gee eee RT Se NSO SEE See ee One ane ent ee Paneer J 


The BACKSPACE statement causes the data set associated with a to Lkack+ 
Space one record. If the data set associated with a is already at its 
beginning, execution of this statement has no effect. For further 
information, see the FORTRAN IV programmers’ guides listed in the 
Preface. 


DIRECT ACCESS INPUT/OUTPUT STATEMENTS 


LOO ALT COTOS ES SLT TET AY AD LAOREES-CGREST LEUAEED SUROE ORNS Net UAE OEP CLEED “WRENS NAD GLAD Wek OWED GEMS rete WORT EE WAT SNES EEE RSE AUNEP CUED EE CUED RE GY me 


There are four direct access input/output statements: READ, WRITE, 
DEFINE FILE, and FIND. The READ and WRITE statements cause transfer of 
data into or out of internal storage. These statements allow the user 
to specify the location within a data set from which data is to be read 
or into which data is to be written. 


The DEFINE FILE statement specifies the characteristics of the data 
set(s) to be used during a direct access operation. The FIND statement 
overlaps record retrieval from a direct access device with computation 
in the program. In addition to these four statements, the FORMAT state- 
ment (described previously) specifies the form in which data is to be 
transmitted. The direct access READ and WRITE statements and the FIND 
statement are the only input/output statements that may refer to a data 
set reference number defined by a DEFINE FILE statement. 


Each record in a direct access data set has a unigue record number 
associated with it. The programmer must specify in the READ, WRITE, and 
FIND statements not only the data set reference number, as for sequen- 
tial input/output statements, but also the number of the record to be 
read, written, or found. Specifying the record number permits opera- 
tions to be performed on selected records of the data set, instead of on 
records in their sequential order. 


The number of the record physically following the one just processed 
is made available to the program in an integer variable known as the 
associated variable. Thus, if the associated variable is used in a READ 
Or WRITE statement to specify the record number, sequential processing 
is automatically secured. The associated variable is specified in the 
DEFINE FILE statement, which also gives the number, size, and type of 
the records in the direct access data set. 


DEFINE FILE STATEMENT 


wl 
vat 


The DEFINE FILE statement describes the characteristics of any data 
set to be used during a direct access input/output operation. To use 
the direct access READ, WRITE, and FIND statements in a program, the 
data set{s) must ke described with a DEFINE FILE statement. Each data 
set must be described once, and this description may appear once in each 
program Or subprogram. Subsequent descriptions have no effect. 
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The DEFINE FILE statement must logically precede (i.e., must be 
"executed" prior to) any input/output statement referring to the data 
set described in the DEFINE FILE statement. 


‘ 
|] General Form | 


Where: a represents an integer constant that is the data set 
reference number. 


m represents an integer constant that specifies the number 
of records in the data set associated with a. 


r represents an integer constant that specifies the maximum 
Size of each record associated with a. The record size is 
measured in characters (bytes), storage locations (bytes), 
Or storage units (words). (A storage unit is the number of 
storage locations divided by four and rounded to the next 
highest integer.) The method used to measure the record 
Size depends upon the specification for f£. 


| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
! 
£ specifies that the data set is to be read or written eith- | 
er with or without format control; f may be one of the fol- | 
lowing letters: | 
! 

either with or without format control. The maximum record | 
Size is measured in number of storage locations (bytes). | 
| 

| 

| 

| 

| 

| 

| 

| 

| 

| 

| 

| 

| 

| 

| 

| 


E indicates that the data set is to be read or written 
under format control (as specified by a format statement). 
The maximum record size is measured in number of charac- 
ters (bytes). 


U indicates that the data set is to be read or written 
without format control. The maximum record size is mea- 
sured in number of storage units (words). 


v represents a nonsubscripted integer variable called an 
associated variable. At the conclusion of each read or 
write operation, v is set to a value that points to the rec- 
ord that immediately follows the last record transmitted. 
At the conclusion of a find operation, v is set to a value 


| 
| 
| 
| 
| 
| 
{ 
| 
| 
] 
l 
| 
| 
| 
| L indicates that the data set is to be read or written 
| 
| 
] 
| 
| 
| 
| 
] 
l 
| 
| 
| that points to the record found. 
L 


The associated variable cannot appear in the I/O list of a READ or 
WRITE statement for a data set associated with the DEFINE FILE 
statement. 


Example: 
DEFINE FILE 8(50,100,L,12),9(100,50,L,d3) 


This DEFINE FILE statement describes two data sets, referred to by 
data set reference numbers 8 and 9. The data in the first data set con- 
Sists of 50 records, each with a maximum length of 100 storage loca- 
tions. The L specifies that the data is to be transmitted either with 
or without format control. 12 is the associated variable that serves as 
a pointer to the next record. 
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The data in the second data set consists of 100 records, each with a 
maximum length of 50 storage locations. The L specifies that the data 
is to be transmitted either with or without format control. J3 is the 
associated variable that serves as a pointer to the next record. 


If an E is substituted for the L in the preceding DEFINE FILE state- 
ment, a FORMAT statement is required and the data is transmitted under 
format control. If the data is to be transmitted without format con- 
trol, the DEFINE FILE statement can be written as: 


DEFINE FILE 8(50, 25,U,12),9(100,13,U,J3) 


NIRFCT ACCESS FROGRAMMING CONSIDERATIONS 


When programming for direct access input/output operations, the user 
must establish a correspondence between FORTRAN records and the records 
described by the CLEFINE FILE statement. All conventions of FORMAT con- 
trol discussed in the section “FORMAT Statement” are applicable. 


For example, tc process the data set described by the statement: 
DEFINE FILE 8(10, 48,L,K8) 


the FORMAT statement used to control the reading Or writing could not 
specify a record lcnger than 48 characters. The statements: 


FORMAT (4F12.1) or 
FORMAT (112, 9F4. 2) 


define a FORTRAN record that corresponds to those records described by 
the DEFINE FILE statement. The records can also be transmitted under 
FORMAT control by substituting an E for the L and rewriting the DEFINE 
FILE statement as: 


DEFINE FILE 8(10, 48,E, K8) 


To process a direct access data set without format control, the number 
of storage locations specified for each record must be greater than 
or equal to the maximum number of storage locations in a record to be 
written by any WRITE statement referencing the data set. For example, 
if the I/O list of the WRITE statement specifies transmission of the 
contents of 100 storage locations, the DEFINE FILE statement can be 
either: 


DEFINE FILE 8(50,100,L,K8) or 
DEFINE FILE 8(50, 25,U, K8) 


Programs may share an associated variable as a COMMON variable. The 
following example shows how this can be accomplished. | 


COMMON IUAR SUBROUTINE SUBI(A, B) 


DEFINE FILE 8(100,10,1, IUAR) COMMON IUAR 
ITEMP=IUAR | 


CALL SUBI (ANS, ARG) 


6 


In this example, the program and the subprogram share the associated 
variable IUAR. An input/output operation that references data set 8 and 
is performed in the subroutine causes the value of the associated vari- 
able to be changed. The associated variable is then tested in the main 
program in statement 8. | 


READ STATEMENT 


The READ staterent causes data to be transferred from a direct access 
device into internal storage. The data set being read must be defined 
With a DEFINE FILE statement. 


] General Form 
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READ (a'r, b, ERR=d) list 


Where: ais an integer constant or unsigned integer variable that 
is of length 4 and represents a data set reference number; a 
must be followed by an apostrophe ('). 


ris an integer expression that represents the relative 
position of a record within the data set associated with a. 


b is optional and, if given, is either the statement number 
of the FORMAT statement that describes the data being read 
or the name of an array that contains an object time format. 


ERR=d is optional and d is the statement number to which 
control is given when a device error condition is encoun- 
tered during data transfer from device to storage. 


list is optional and is an I/O list. 


a a es iin a a cms heen mali tum “in el i cance ie mmc’ 2s Bn “tien el 


The I/O list must not contain the associated variable defined in the 
DEFINE FILE staterent for data set a. 


Example: 


DEFINE FILE 8(500,100,L,1D1),9(100, 28, L, 1D2) 
DIMENSION M(10) 


10 FORMAT (5120) 
9 READ (8'16,10) (M(K),K=1,10) 
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READ statement 9 transmits data from the data set associated with 
data set reference number 8, under control of FORMAT statement 10; 
transmission begins with record 16. Ten data items of 20 characters 
each are read as specified by the I/O list and FORMAT statement 10. Two 
records are read to satisfy the I/O list, because each record contains 
only five data items (100 characters). The associated variable ID1 is 
set to a value of 18 at the conclusion of the operation. | 


READ statement 13 transmits data from the data set associated with 
data set reference number 9, without format control; transmission begins 
with record 26. Data is read until the I/O list for statement 13 is 
satisfied. Because the DEFINE FILE statement for data set 9 specified 
the record length as 28 storage locations, the I/O list of statement 13 
calls for the same amount of data (the seven variables are type real and 
each occupies four storage locations). The associated variable ID2 is 
set to a value of 27 at the conclusion of the operation. If the value 
of ID2 is unchanged, the next execution of statement 13 reads record 32. 


The DEFINE FILE statement in the previous example can also be written 
as: 


DEFINE FILE 8(500,100,E,1D1),9(100,7,U, ID2) 


The FORMAT statement may also control the point at which reading 
Starts. For example, if statement 10 in the example is 


10 FORMAT (//5120) 


records 16 and 17 are skipped, record 18 is read, records 19 and 20 are 
skipped, record 21 is read, and ID1 is set to a value of 22 at the con- 
clusion of the READ operation in statement 9. 


WRITE STATEMENT 


The WRITE statement causes data to be transferred from internal 
storage to a direct access device. The data set being written must Le 
defined with a DEFINE FILE statement. 
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Where: a is an integer constant or unsigned integer variable that 
is of length 4 and represents a data set reference number; 
a must ke followed by an apostrophe ('). 


r is an integer expression that represents the relative 
position of a record within the data set associated with a. 


b is optional and, if given, is either the statement number 
of the FORMAT statement that describes the data being writ- 
ten or the name of an array that contains an object time 
format. | 


list is optional and is an I/O list. 


Example: 


DEFINE FILE 8(500,100,L,ID1),9(100, 28,L, 1D2) 
DIMENSION M(10) 


ID2=21 


10 FORMAT (5120) 
8 WRITF (8'16,10) (M{K),K=1,10) 


11 WRITE (9'ID2+5) A,B,C,D,E,F,G 


WRITE statement 8 transmits data into the data set associated with 
the data set reference number 8, under control of FORMAT statement 10; 
transmission begins with record 16. Ten data items of 20 characters 
each are written as specified by the I/O list and FORMAT statement 10. 
Two records are written to satisfy the IYO list because each record con- 
tains 5 data items (100 characters). The associated variable ID1 is set 
to a value of 18 at the conclusion of the operation. 


WRITE statement 11 transmits data into the data set associated with 
data set reference number 9, without format control; transmission begins 
with record 26. The contents of 28 storage locations are written as 
specified by the IVC list for statement 11. The associated variable ID2 
is set to a value of 27 at the conclusion of the operation. Note the 
correspondence between the records described (28 storage locations per 
record) and the number of items called for by the I/VO list (7 variables, 
type real, each occupying four storage locations). 


The DEFINE FILE statement in the example can also be written as: 
DEFINE FILE 8(500,100,E,ID1), 9(100,7,U,ID2) 


As with the READ statement, a FORMAT statement may also be used to 
control the pcint at which writing begins. 


FIND STATEMENT 


The FIND statement causes the next input record to be found while the 
present record is being processed, thereby increasing the execution 
speed of the object program. The program has no access to the record 
that was found until a READ statement for that record is executed. 
(There is no advantage to having a FIND statement precede a WRITE 
statement. ) 


] General Form 
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is of length 4 and represents a data set reference number; a 
must be followed by an apostrophe ('). 


rx is an integer expression that represents the relative 
position of a record within the data set associated with a. 


ne EE ce EY CE EE EE ee eee aan Elites eS ED wale apt AED sale: eee gems <i citi wens CNT ile ey <P <iinte COE ime che Re ibe ie ae ee eS ER cue See ate Cee eee ape See ee Gee ae ee EE ae ee ee ee ee ee ee ee ee ee ee ee eee oe ee eee ee ee ee ee 


; 
| 
| 

Where: ais an integer constant or unsigned integer variable that 
| 
| 
| 
| 
J 
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The data set on which the record is being found must be defined with a 
DEFINE FILE statement 


Example: 


10 FIND (8'50) 


15 READ (8°50) A,B 


While the statements between statements 10 and 15 are executed, rec- 
ord 50, in the data set associated with data set reference number 8, is 
found. 


General Example -- Direct Access Operations 


DEFINE FILF 8(1000,72,L,ID8) 
DIMENSION A(100),B(100) ,C (100) ,D(100), £(100), F(100) 


15 FORMAT (6F12.4) 
FIND (8'5) 


ID8=1 

DO 100 I=1,100 

READ (8*ID84+4,15)A(I),B(I),C(1), D(1), E(1), F(L) 
100 CONTINUE 


DO 200 I=1,100 
WRITE (8*ID84+4,15)A(I),B(I), C(I), D(I), E(I), F(I) 
200 CONTINUE 


END 


The general example illustrates the ability of direct access state- 
ments to gather and disperse data in an order designated by the user. 
The first DO loop in the example fills arrays A through F with data from 
the 5th, 10th, 15th,...e, and 500th record associated with data set 
reference number 8. Array A receives the first value in every fifth rec- 
ord, B the second value and so on, as specified by FORMAT statement 15 
and the I/O list cf the READ statement. At the end of the READ opera- 
tion, each record has been dispersed into arrays A through F. At the 
conclusion of the first DO loop, ID8 has a value of 501. 


The second DO loop in the example groups the data items from each 
array, as specified by the I/O list of the WRITE statement and FORMAT 
statement 15. Each group of data items is placed in the data set asso- 
ciated with data set reference number 8. Writing begins at the 505th 
record and continues at intervals of five, until record 1000 is written, 
if ID8 is not changed between the last READ and the first WRITE. 
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DATA INITIALIZATION STATEMENT 
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DATA k1/d1/,k2/d2/,... ,kn/dn/ 


Where: Each k is a list containing variables, subscripted variables 
(in which case the subscripts must be integer constants), or 
array names. Dummy arguments may not appear in the list. 


decimal, logical, or literal), any of which may be preceded 
by i*. 


Fach i is an unsigned integer constant. When the form i* 
appears before a constant, it indicates that the constant is 


| 
| 
| 
| 
| Each @ is a list of constants (integer, real, complex, hexa- 
| 
l 
to be specified i times. 


bes mmm ees cee cess SE ees ets me comet Seen Seem emoees dene eoclens smn anol 


A DATA initialization statement is used to define initial values of 
variables, array elements, and arrays. There must be a one-to-one 
correspondence between the total number of elements specified or implied 
by the list k and the total number of constants specified by the corres- 
ponding list d after application of any replication factors, i. 


This statement cannot precede any other specification statement that 
refers to the same variables or arrays. It also cannct precede an 
IMPLICIT statement. Otherwise, a DATA statement can appear anywhere in 
the program. 


Example 1: 


DIMENSION D(5,10) 
DATA A, B, C/5.0,6.1,7.3/,D,E/25*1. 0, 25*2.0,5.1/ 


The DATA staterent indicates that the variables A, B, and C are to be 
initialized to the values 5.0, 6.1, and 7.3 respectively. In addition, 
the statement specifies that the first 25 variables in the array D are ° 
to be initialized to the value 1.0, the remaining 25 variables in D to 
the value 2.0, and the variable E to the value 5.1. 
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Example 2: 


DIMENSION A(5), B(3,3), L(4) 
DATA A/5*1.0/, B/9*2 O/, L/4*.TRUE./, C/'FOUR!/ 


Explanation: 


The DATA statement specifies that all the variables in the arrays A 
and B are to be initialized to the values 1.0 and 2.0, respectively. 
All the logical variables in the array L are initialized to the value 
«TRUE... The letters T and F may be used as an abbreviation for .TRUE. 
and .FALSE., respectively. In addition, the variable C is initialized 
with the literal data constant FOUR. 


An initially defined variable, or variable of an array, may not be in 


blank common. In a labeled common block, they may be initially defined 
only in a BLOCK DATA subprogram. (See the section "Subprograms.") 
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SPECIFICATION STATEMENTS 


The specification statements provide the compiler with information 
about the nature of the data used in the source program. In addition, 
they supply the information required to allocate locations in storage 
for this data. 


Specification statements must precede statement function definitions, 
which must precede the program part containing at least one executable 
statement. Within the specification statements, any statement describ- 
ing data must precede references to that data. In particular, the 
IMPLICIT statement, if used, must be the first specification statement. 


The specification statement EXTERNAL is described in the section 
"“Subprograms."“ 


DIMENSION STATEMENT 


ge ae rn eee ee ee SD RY PEN eT 7 
| General Form | 
Scat ae Ae se a a aa ns are ee ees | 
| DIMENSION a4 (K1),aalKa), az lK3) 4-020, an (Kn) | 
| | 
| | 
| | 
| Where: a4, Any Aaygeeey An AYE array names. | 
| 
| Kay Kay Kayeee,Kn are each composed of one through seven | 
| unsigned integer constants, separated by commas, represent- | 
| ing the maximum value of each subscript in the array. ky, | 
| through kn may be integer variables of length 4 only when | 
| the CIMENSION statement in which they appear is ina | 
] Subprogram. | 
Dc ea a ae ee 4 


The information necessary to allocate storage for arrays used in the 
source program may Le provided by the DIMENSION statement. The follow- 
ing examples illustrate how this information may be declared. 


Examples: 


DIMENSION A (10), ARRAY (5,5,5), LIST (10,100) 
DIMENSION B(25, 50), TABLE(5, 8, 4) 


TYPE STATEMENTS 


There are two kinds of type statements: the IMPLICIT specification 
Statement, and the explicit specification statements (INTEGER, REAL, 
CCMPLEX, and LCGICAL). 

The IMPLICIT statement enables the user to: 

e Specify the type of a group of variables or arrays according to the 

initial character of their names 
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Specify the amount of storage to be allocated for each variable 
according to the associated type 


The explicit specification statements enable the user to: 


Specify the type of a variable or array according to its particular 
name 


e Specify the amount of storage to be allocated for each variable 
according tc the associated type 


e Specify the dimensions of an array 
e Assign initial data values for variables and arrays 


The explicit specification statement overrides the IMPLICIT state- 
ment, which, in turn, overrides the predefined convention for specifying 


type. 


IMPLICIT STATEMENT 


] IMPLICIT typetslay,aaeece) pees - typer*s (aq, anyone) 


Where: type is one of the following: INTEGER, REAL, COMPLEX, or 
LOGICAL. 


*S 1S Optional and represents one of the permissible length 
specifications for its associated type. 


A1¢ Aayeee are Single alphabetic characters separated by 
commas, Or a range of characters drawn from the set A, B, 
eee LZ, $, in that order. The range is denoted by the first 
and last characters of the range separated by a minus Sign 
(e.g., (A-D)). 


amet eee Se Ke Sh ee Smee SS ee ee SR OED 


| 
| 
| 
| 
| 
| 
| 


The IMPLICIT specification statement must be the first statement ina 
main program and the second statement ina subprogram. There can be 
only one IMPLICIT statement per program or subprogram. The IMPLICIT 
specification statement enables the user to declare the type of the 
variables appearing in his program (1.e., integer, real, complex, or 
logical) by specifying that variables beginning with certain designated 
letters are of a certain type. Furthermore, the IMPLICIT statement 
allows the programmer to declare the number of locations (bytes) to be 
allocated for each in the group of specified variables. The types that 
a variable may assume, along with the permissible length specifications, 
are as follows: 


Type Length Specification 

INTEGER 2or 4 (standard length is 4) 
REAL 4 or 8 (standard length is 4) 
COMPLEX @ or 16 (standard iength is ®) 
LOGICAL 1 or 4 (standard length is 4) 


For each type there is a corresponding standard length specification. 
If this standard length specification (for its associated type) is 
desired, the *s may be omitted in the IMPLICIT statement. That is, the 
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variables will assume the standard length specification. For each type 
there is also a corresponding optional length specification. If this 
optional length specification is desired, then the *s must be included 
Within the IMPLICIT statement. 


Example 1: 


IMPLICIT REAL (A-H, O-S), INTEGER (I-N) 


All variables beginning with the characters I through N are declared 
aS INTEGER. Since no length specification was explicitly given (i.e., 
the *sS was omitted), four storage locations (the standard length for 
INTEGER) are allocated for each variable. 


All other variables (those beginning with the characters A through H, 
C through Z, and $) are declared as REAL with four storage locations 
allocated for each. 


Note that the statement in example 1 performs the same function of 
typing variables as the predefined convention (see “Type Declaration by 
the Predefined Specification"). 


Example 2: 


IMPLICIT INTEGER*2(A-H), REAL*8(I-K), LOGICAL(L,M,N) 


Explanation: 


All variables keginning with the characters A through H are declared 
as integer with two storage locations allocated for each. All variables 
beginning with the characters I through K are declared as real with 
eight storage locations allocated for each. All variables beginning 
With the characters L, M, and N are declared as logical with four loca- 
tions allocated for each. 


Since the remaining letters of the alphabet, namely, O through Z and 
S, are left undefined by the IMPLICIT statement, the predefined conven- 
tion will take effect. Thus, all variables beginning with the charac- 
ters O through Z and $§ are declared as real, each with a standard length 
of four locations. 


Example 3: 


IMPLICIT COMPLEX*16(C-F) 
Explanation: 


All variables Leginning with the characters C through F are declared 
as complex, each with eight storage locations reserved for the real part 
of the complex data and eight storage locations reserved for the 
imaginary part. The types of the variables beginning with the charac- 
ters A, B, G through Z, and $ are determined by the predefined 
convention ! 
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EXPLICIT SPECIFICATION STATEMENTS 


|] General Form 
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Where: Type is INTEGER, REAL, LOGICAL, or COMPLEX. 


¥S,_)¥*¥S4- *Sa,eee,*Syn are Optional. Each s represents one of 
the permissible length specifications for its associated 
type. 


F 

| 

| 

| 

| 

| 

| 

| 
| Ay Dygece,Z are variable, array, or function names (see the 
| section "“Subprograms"™) | 

| 

| 

| 

| 

| 

| 

| 

L 


(ka), (Ka) ,¢ee<, (Kn) are optional and give dimension informa- 
tion for arrays. Each k is composed of one through seven 
unsigned integer constants, separated by commas, represent- 
ing the maximum value of each subscript in the array. Each 
k may be an unsigned integer variable of length 4 only when 
the type statement in which it appears is in a subprogram. 


1X41 q1B2lyeeesg/¥yn/ are Optional and represent initial data 
values. 


a ee Ee Sn | 


The explicit specification statements declare the type (INTEGER, 
REAL, COMPLEX, or LOGICAL) of a particular variable or array by its 
name, rather than by its initial character. This differs from the other 
ways of specifying the type of a variable or array (1i.e., predefined 
convention and the IMPLICIT statement). In addition, the information 
necessary to allocate storage for arrays (dimension information) may be 
included within the statement. 


Initial data values may be assigned to variables or arrays by use of 
/X¥n/ where Xn is a constant or list of constants separated by commas. 
Lists of constants are used only to assign initial values to array ele- 
ments. x successive occurrences of the same constant can be represented 
by the form r*constant. If initial data values are assigned to an array 
in an explicit specification statement, the dimension information for 
the array must be in the explicit specification statement or in a pre- 
ceding DIMENSION or COMMON statement. An initial data value may not be 
assigned to a function name. 


Initial data values cannot be assigned to variables or arrays in 
blank common. The BLOCK DATA subprogram must be used to assign initial 
values to variables and arrays in labeled common. 


In the same manner in which the IMPLICIT statement overrides the pre- 
defined convention, the explicit specification statements override the 
IMPLICIT statement and predefined convention. If the length specifica- 
tion is omitted (i.e.,*#s), the standard length per type is assumed. 


Example 1: 
INTEGER*2 ITEM/76/, VALUE 


Explanation: 


This statement declares that the variables ITEM and VALUE are of type 
integer, each with two storage locations reserved. In addition, the 
variable ITEM is initialized to the value 76. 
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Example 2: 
COMPLEX C,D/(2.1,4.7)/, E*¥16 


This statement declares that the variables C, D, and E are of type 
complexe Since no length specification was explicitly given, the stan- 
dard length is assumed. Thus, C and D each have eight storage locations 
reserved (four for the real part, four for the imaginary part) and D is 
initialized to the value (2.1,4.7). In addition, 16 storage locations 
are reserved for the variable E. Thus, if a length specification is 
explicitly written, it overrides the assumed standard length. 


Example 3: 
REAL*8 BAKER, HOLD, VALUE*4, ITEM(5, 5) 


Explanation: 


This statement declares that the variables BAKER, HOLD, VALUE, and 
the array named ITEM are of type real. In addition, it declares the 
Size of the array ITEM. The variables BAKER and HOLD have eight storage 
locations reserved for each; the variable VALUE has four storage loca- 
tions reserved; and the array named ITEM has 200 storage locations 
reserved (eight for each variable in the array). Note that when the 
length is associated with the type (e.g., REAL*8), the length applies to 
each variable in the statement unless explicitly overridden (as in the 
case Of VALUE*4). 


Example 4: 
REAL A(5, 5)/20*6. 9E2, 5*1.0/, B(100)/100*0. 0/7, TEST*8 (5) /5*0,. OY 


This statement declares the size of each array, A and B, and their 
type (real). The array A has 100 storage locations reserved (four for 
each variable in the array) and the array B has 400 storage locations 
reserved (four for each variable). In addition, the first 20 variables 
in the array A are initialized to the value 6.9E2 and the last five 
variables are initialized to the value 1.0. All 100 variables in the 
array B are initialized to the value 0.0. The array TEST has 40 storage 
locations reserved (eight for each variable). In addition, each vari- 
able is initialized to the value 0.0. 


DOUBLE PRECISION STATEMENT 
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Where: a,b,ece,zZ represent variable, array, or function names (see 
the section “Subprograms"). 


(ka), (Ka), «e, (Kn) are optional. Each k is composed of one 
through seven unsigned integer constants, separated by com- 
mas, that represent the maximum value of each subscript in 

the array. 
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The DOUBLE PRECISION statement explicitly specifies that the vari- 
ables a,b,c,..-. are of type double precision. This statement overrides 
any specification of a variable made by either the predefined convention 
or the IMPLICIT statement. This specificaticn is identical to that of 
type REAL*8. This statement cannot be used to define initial data 
values 


In addition, FUNCTION subprograms may be typed double precision as 
follows: 


DOUBLE PRECISION FUNCTION name (a1,a2y83y++ey@n) 


COMMON STATEMENT 


ED A RD ER ENE RE END AR EE AND ED EY AEE AAS EN EE EN NE EES EE A A A ES A ED AS NS ES ED NS SE NS SATS EES GENES EID UES UTS SR SCE “OY EENS EEE GENES <emem SEED GREE GENE QUES GEE CURED GREED GENTE REE GEN GUE CUE GES GR GREEN GEEND ce ce SERED eS CNIS ES mT 


COMMON /r/a (Ky), blka),00./4¥/C(K3) ,d(k )yeee 


Where: a,DygeeeeCydeee are variable names or array names that can- 
not ke dummy arguments. 


KapKayeeeKa-k «+e are optional and are each composed of one 
through seven unsigned integer constants, separated by conm- 
mas, representing the maximum value of each subscript in the 
array. 


one through six alphameric characters, the first of which is 
alphabetic. These names must always be embedded in slashes. 


| 
] J/v/eee represent optional common block names consisting of 


The COMMON statement is used to define a storage area that can be 
referred to by a calling program and one or more supoprograms, and to 
specify the names of variables and arrays to be placed in this area. 
Therefore, variables or arrays that appear in a calling program or sub- 
program can be made to share the same storage locations with variables 
Or arrays in other subprograms. Also, a common area can be used to 
implicitly transfer arguments between a calling program and a subpro- 
gram. Arguments passed in common are subject to the same rules with. 
regard to type, length, etc., as arguments passed in an argument list 
(see the section “Subprograms"). 


If more than one COMMON statement appears in a calling program or 
subprogram, the entries in the statements are cumulative. Redundant 
entries are not permitted. 


Although the entries in a COMMON statement can contain dimension 
information, okject-time dimensions may never be used. 


The length of a common area can be increased by using an EQUIVALENCE 
statement. 


Since the entries in a common area share storage locations, the order 


in which they are entered is significant. Consider the following 
example: 
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Example: 
Calling Program 


COMMON A, B, C, R(100) 
REAL A,B,C 
INTEGER R 


CALL MAPMY (...) 


Subprogram 


SUBROUTINE MAPMY (...) 


COMMON X, Y, Z, S100) 
REAL X,Y,Z 
INTEGER S 


Explanation: 


In the calling program, the statement COMMON A,B,C,R(100) would cause 
412 storage locations (four locations per variable) to be reserved in 
the following order: 


Beginning ' A B C 1 Layout of 
of common | 4 lJccations 4 locations 4 locations ] storage 
area | 

}-------------------------------------------------- : 

] R(1) se 8 R(100) { 

] 4 locations 4 locations | 

ee eee ee ee eee ee oe a eee ee J 


The statement COMMON X, Y, Z, S(100) would then cause the variables 
X, Y, Z and S(1)...S(100) to share the same storage space as A, B, C, 
and R(1)...R{(100), respectively. Note that values for X, Y, Z, and S(1) 
eee5(100), because they occupy the same storage locations as A, B, C, 
and R(1)...R(100), do not have to be transmitted in the argument list of 
a CALL statement. 


PLANK AND LAPELED COMMON 


In the preceding example, the common storage area (common block) is 
called a blank common area. That is, no particular name was given to 
that area of storage. The variables that appeared in the COMMON state- 
ments were assigned locations relative to the beginning of this blank 
common area. However, variables and arrays may be placed in separate 
common areas. Each of these separate areas (or blocks) is given a name 
consisting of one through six alphameric characters (the first of which 
is alphabetic); those blocks which have the same name occupy the same 
storage space. This permits a calling program to share one common block 
with one subprogram and another common block with another subprogram and 
also facilitates program documentation. 


Those variables that are to be placed in labeled (or named) common 
are preceded by a common block name enclosed in slashes. For example, 
the variables A,B, and C will be placed in the labeled common area, 
HOLD, by the follcwing statement: 


COMMON/HOLD/A, B,C 


In a COMMON statement, blank common may be distinguished from labeled 
common by preceding the variables in blank common by two consecutive 
Slashes or, if the variables appear at the beginning of the COMMON 
Statement, by omitting any block name. For example, in the following 
statement: 
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COMMON A, B, C /ITEMS/ X, Y¥, Z2//D, E, F 


the variables A, B, C, D, E, and F will be placed in blank common in 
that order; the variables X, Y, and Z will be placed in the common area 
labeled ITEMS. , 


Blank and labeled common entries appearing in COMMON statements are 
cumulative throughout the program. For example, consider the following 
two COMMON statements: 


COMMON A, B, C /R/ D, E /S/ F 
COMMON G, H /S/ I, J /R/P//W 


These two statements have the same effect as the single statement: 
COMMON A, B, C, G, H, W /R/ D, E, P /S/ F, I, J 
Example: 


Assume that A, B, C, K, X, and Y each occupy four locations of 
storage, H and G each occupy eight locations, and D and E each occupy 
two locations. 


Calling Program Subprogram 


° SUBROUTINE MAPMY(...) 


COMMON H, A /R/ X, D // B Z 
P COMMON G, Y, C /R/ K, E 
CALL MAPMY(...) | : 


Explanation: 


In the calling program, the statement COMMON H, A /R/ xX, D //B causes 
16 locations (four locations each for A and B, and eight for H) to be 
reserved in blank common in the following order: 


Ct Oct ee ae ee ge Oe ee ep ae gi OR ye Pe ew 1 
Beginning |] H A B | 
of blank | 8 locations 4 locations 4 locations | 
common i | 

}~--~--~----------------------------------------------------- 


and also causes six locations (four for X, and two for D) to be reserved 
in the labeled common area R in the following order: 


aaa a aa aa aaa aaa a a a 7 
Beginning | X D | 
of labeled] 4 locations 2 locations | 
common R_ | | 

}~-~---------------------------------------------------------- { 
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The statement COMMON G,Y,C/R/K,E appearing in the subprogram MAPMY 
would then cause the variables G, Y, and C to share the same storage 
Space (in blank common) as H, A, and B, respectively. It would also 
cause the variables K and E to share the same storage space (in labeled 
common area R) as X and D, respectively. 


ARRANGEMENT OF VARIABLES IN COMMON 


Variables in a common block need not be aligned properly. However, 
considerable object-time efficiency is lost unless the programmer 
ensures that all of the variables have proper boundary alignment. 


Proper alignment is achieved either by arranging the variables ina 
fixed descending order according to length, or by constructing the block 
so that dummy variables force proper alignment. If the fixed order is 
used, the variables must appear in the following order: 


length of 16 (complex) 

length of 8 (complex or real) 

length of 4 (real or integer or logical) 
length of 2 (integer) 

length of 1 (logical) 


If the fixed order is not used, proper alignment can be ensured by 
constructing the block so that the displacement of each variable can be 
evenly divided by the reference number associated with the variable. 
(Displacement is the number of storage locations (bytes) from the begin- 
ning of the block to the first storage location of the variable.) The 
following list shows the reference number for each type of variable: 


Type of Length Reference 
Variable Specification _Number _ 
Logical a 1 
4 Ly 
Integer 2 Z 
uy | ny 
Real y- 4 
8 8 
Complex 8 8 
16 8 


The first variable in every common block is positioned as though its 
length specification were eight. Therefore, a variable of any length 
may be the first assigned within a block. To obtain the proper align- 
ment for other variables in the same block, it may be necessary to add a 
dummy variable to the block. For example, the variables A, I, and CMPLX 
are REAL*4, INTEGER*4, and COMPLEX*8, respectively, and form a COMMON 
block that is defined as: 


COMMON A, I, CMPLX 


Then, the displacement of these variables within the block is illus- 
trated as follows: 
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}< eee A--------- >|< -——— I---—~—-—--~—-—— >|< S eabenhentaeented CMPLX-——--~——--—---— > 
| 4 storage l 4 storage | 8 storage 

] locations | locations | locations 

| l 

Vv Vv Vv 

displacement displacement displacement 

0 storage | 4 storage 8 storage 

locations locations locations 


The displacements of I and CMPLX are evenly divisible by their reference 
numbers. However, if I were an integer with a length specification of 
2, then CMPLX is not properly aligned (its displacement of 6 is not 
evenly divisible by its reference number of 8). In this case, proper 
alignment is ensured by inserting a dummy variable with .a length speci- 
fication of 2 either between A and I, or between I and CMPLX. 


EQUIVALENCE STATEMENT 
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| 

i that may be subscripted. The subscripts may have two forms: 
| If the variable is singly subscripted, it refers to the 

| position of the variable in the array (i.e., first variable, 
] 25th variable, etc.). If the variable is multi-subscripted, 
J it refers to the position in the array in the same fashion 

l 


‘ 
| 

4 

| 

| 

Where: a, b, cy, d, e, £,--- are variables (not dummy arguments) | 
| 

| 

| 

| 

| 

as the position is referred to in an arithmetic statement. | 

J 


The EQUIVALENCE statement provides the option for controlling the 
allocation of data storage within a single program unit. In particular, 
when the logic of the program permits it, the number of storage loca- 
tions used can be reduced by causing locations to be shared by two or 
more variables of the same or different types. Equivalence between 
variables implies storage sharing only, not mathematical equivalence. 


Since arrays are stored in a predetermined order (see “Arrangement of 
Arrays in Storage"), equivalencing two elements of two different arrays 
may implicitly equivalence other elements of the two arrays. The 
EQUIVALENCE statement must not contradict itself or any previously esta- 
blished equivalences. 


Two variables in one common block or in two different common blocks 
cannot be made equivalent. However, a variable in a program or a sub- 
program can be made equivalent to a variable in a common block. If the 
variable that is equivalenced to a variable in the common block is an 
element of an array, the implicit equivalencing of the rest of the ele- 
ments of the array may extend the size of the common block (see example 
2). The size of the common block must not be extended so that elements 
are added before the beginning of the established common block. 


Example 1: 


Assume that in the initial part of a program, an array C of size 
100x100 is needed; in the final stages of the program C is no longer 
used, but arrays A and B of sizes 50x50 and 100, respectively, are used. 
The elements of all three arrays are of the type REAL*4. Storage space 
can then be saved by using the statements: 
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DIMENSION C(100,100), A(50,50), B(100) 
EQUIVALENCE (C(1), A(1)), (C(2501), B(1)) 


The array A, which has 2500 elements, can occupy the same storage as 
the first 2500 elements of array C since the arrays are not both needed 
at the same time. Similarly, the array B can be made to share storage 
with elements 2501 to 2600 of array C. 


Example 2: 


DIMENSION B(5), C(10, 10), D5, 10, 15) 
EQUIVALENCE (A, B(1), C(5,3)), (D(5,10,2), E) 


This EQUIVALENCE statement specifies that the variables A, B(1), and 
C(5,3) are assigned the same storage locations and that variables D(5, 
10,2) and E are assigned the same storage locations. It also implies 
that B(2) and C(6,3), etc., are assigned the same storage locations. 
Note that further equivalence specification of B(2) with any element of 
array C other than C(6,3) is invalid. 


The designations C(5,3) and D(5,10,2) could have been replaced with 
the designations C(25) and D(100) and the effect would have been the 
same. 


Example 3: 


COMMON A, B, C 
DIMENSION D(3) 
EQUIVALENCE (B,D(1)) 


Explanation: 


This would cause a common area to be established containing the vari- 
ables A, B, and C. The EQUIVALENCE statement would then cause the vari- 
able D(1) to share the same storage location as B, D(2) to share the 
Same storage location as C, and D(3) would extend the size of the common 
area, in the following manner: 


A (lowest location of the common area) 
B, D(1) 
C, D2) 

D(3) (highest location of the common area) 


The following EQUIVALENCE statement is invalid: 


COMMON A, B, C 
DIMENSION D(3) 
EQUIVALENCE (B, D(3)) 


because it would force D(1) to precede A, as follows: 


D(1) 
A, D2) (lowest location of the common area) 
B, D3) 
C (highest location of the common area) 
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STORAGE ARRANGEMENT OF VARIABLES IN EQUIVALENCE GROUPS 


Variables in an equivalence group may be in any order in main 
storage. However, considerable object-time efficiency is lost unless 
the programmer ensures that all of the variables have proper boundary 
alignment. 


Proper alignment is achieved either by arranging the variables ina 
fixed, descending order according to length, or by constructing the 
group so that dummy variables force proper alignment. If the fixed 
order is used, the variables must appear in the following order: 


6 (complex) 


length of 1 

length of 8 (complex or real) 

length of 4 (real or integer or logical) 
length of 2 (integer) 

length of 1 (logical) 


If the fixed order is not used, proper alignment can be ensured by 
constructing the group so that the displacement of each variable in the 
group can be evenly divided by the reference number associated with the 
variable. (Displacement is the number of storage locations (bytes) from 
the beginning of the group to the first storage location of the vari- 
able.) The reference numbers for each type of variable are given in the 
section "COMMON Statement." The first variable in each group is posi- 
tioned as if its length specification were eight. 


For example, the variables A, I, and CMPLX are REAL*4, INTEGER*4, and 
COMPLEX*8, respectively, and are defined as: 


DIMENSION A(10), 1(16), CMPLX(5) | 
EQUIVALENCE (A(1), I(7), CMPLX(1)) 


Then, the displacement of these variables within the group is illus- 
trated as follows: 


] 1I(1)<---------- q------------------ 64 storage locations--—--~-~--- >I (16) 
| | 

| | AC1)<---~--------- 4Q storage locations-——-—-—---- >A (10) 
| 

] | CMPLX(1)<--------40 storage locations-—-—-—-- >CMPLX(5) 
] | 

Vv Vv 

displacement displacement 

0 storage 24 storage 

locations locations 


The displacements of A and CMPLX are evenly divisible by their reference 
numbers. However, if the EQUIVALENCE statement were written as 


EQUIVALENCE (A(1), I(6), CMPLX(1)) 


then CMPLX is not properly aligned (its displacement of 20 is not evenly 
divisible by its reference number of 8). 


Note that this discussion applies solely to the manner in which the 
equivalence group is arranged in storage This arrangement is not 
affected by the order in which the variable and array names are listed 
in the EQUIVALENCE statement For example, the statement EQUIVALENCE 
(A(1),1(7),CMPLX(1)) has exactly the same effect as EQUIVALENCE 
(CMPLX(1),A(1),1(7)). 
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SUBPROGRAMS 


It is sometimes desirable to write a program which, at various 
points, requires the same computation to be performed with different 
data for each calculation. It would simplify the writing of that pro- 
gram if the statements required to perform the desired computation could 
be written only once and then could be referred to freely, with each 
subsequent reference having the same effect as though these instructions 
were written at the point in the program where the reference was made. 


For example, to take the cube root of a number, a program must be 
written with this object in mind. If a general program were written to 
take the cube root of any number, it would be desirable to be able to 
combine that program (or subprogram) with other programs where cube root 
calculations are required. 


The FORTRAN language provides for the above situation through the use 
of subprograms. There are two classes of subprograms: FUNCTION subpro- 
grams and SUBROUTINE subprograms. In addition, there is a group of 
FORTRAN supplied subprograms (see Appendix C). Functions differ from 
SUBROUTINE subprograms in that they return at least one value to the 
calling program, whereas SUBROUTINE subprograms need not return any. 


Statement functions are also discussed in this section since they are 
Similar to FUNCTION subprograms. The difference is that subprograms are 
not in the same program unit as the program unit referring to them, 
while statement function definitions and references are in the same pro- 
gram unit. 


NAMING SUBPROGRAMS 


A subprogram name consists of from one through six alphameric charac- 
ters, the first of which must be alphabetic. A subprogram name may not 
contain special characters (see Appendix A). The type of a function 
determines the type of the result that can be returned from it. 


e Type Declaration of a Statement Function: Such declaration may he 
accomplished in one of three ways: by the predefined convention, by 
the IMPLICIT statement, or by the explicit specification statements. 
Thus, the rules for declaring the type of variables apply to state- 
ment functions. 


° Type Declaration of FUNCTION Subprograms: The declaration may be 
made by the predefined convention, by the IMPLICIT statement, by an 


explicit specification in the FUNCTION statement, or by an explicit 
specification statement within the FUNCTION subprogram. 


The type of a SUBROUTINE subprogram cannot be defined because the 
results that are returned to the calling program are dependent only on 
the type of the variable names appearing in the argument list of the 
calling program and/or the implicit arguments in common. 
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FUNCTIONS 


A function is a statement of the relationship between a number of 
variables To use a function in FORTRAN, it is necessary to: 


1. Define the function (i.e , specify which calculations are to be 
performed) 


2. Refer to the function by name where required in the program 
Function Definition 
There are three steps in the definition of a function in FORTRAN: 


1. The function must be assigned a unique name by which it may be 
called (see the section “Naming Subprograms") 


2. The dummy arguments of the function must be stated 
3. The procedure for evaluating the function must be stated 


Items 2 and 3 are discussed in detail in the sections dealing with 
the specific subprogram (e.g., “Statement Functions," “FUNCTION Subpro- 
grams," etc.). 


Function Reference 


When the name of a function, followed by a list of its arguments, 
appears in any FORTRAN expression, it references the function and causes 
the computations to be performed as indicated by the function defini- 
tion. The resulting quantity replaces the function reference in the 
expression, and assumes the type of the function. The type of the name 
used for the reference must agree with the type of the name used in the 
definition. 


STATEMENT FUNCTIONS 


A statement function definition specifies operations to be performed 
whenever that statement function name appears as a function reference in 
another statement in the same program unit. 


|name(az, aay Aayeeee an) = Expression 


= 
D> 
1) 
KX 
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: name is the statement function name (see the section "Naming 
Subprograms"). 


Age AaeAzyeeeygan are dummy arguments. They must be unique 
(within the statement) nonsubscripted variables. 


expression is any arithmetic or logical expression that does 
not contain subscripted variables Any statement function 
appearing in this expression must have been defined 
previously. 
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The expression to the right of the equal sign defines the operations 
to be performed when a reference to this function appears in an assign- 
ment statement. The expression defining the function must not contain a 
reference to the function. | 


The dummy arguments enclosed in parentheses following the function 
name are dummy variables for which the arguments given in the function 
reference are substituted when the function reference is encountered. 
The same dummy arguments may be used in more than one statement function 
definition, and may be used as variables outside the statement function 
definitions. An actual argument in a statement function reference may 
be any expression of the same type as the corresponding dummy argument. 


All statement function definitions to be used in a program must pre- 
cede the first executable statement of the program. 
Example: The statement: 
FUNC(A,B) = 3.*AtB¥*¥2.+X+Y+Z 
defines the statement function FUNC, where FUNC is the function name and 
A and B are the dummy arguments. The expression to the right of the 
equal sign defines the operations to be performed when the function 
reference appears in an arithmetic statement. 
The function reference might appear in a statement as follows: 
C = FUNC(D,E) 
This is equivalent to: 
C = 3.*D+E*¥*¥2.+X+Y+Z 
Note the correspondence between the dummy arguments A and B in the func- 


tion definition and the actual arguments D and E in the function 
reference. 


Examples: 
Valid statement function definitions and statement function 
references: 
Definition Reference 
SUM(A, B,C, D) = A+B+C+D NET = GROS-SUM (TAX, FICA, HOSP, STOCK) 
FUNC(Z) = A+X*Y*¥Z ANS = FUNC (RESULT) 


VALID(A,B) = .NOT. A .OR. B. VAL = TEST .OR. VALID(D,E) 
BIGSUM = SUM(A,B,SUM(C,D,E,F) ,G) 


Invalid statement function definitions: 


SUBPRG(3,0,K) =3*I+0**3 (Arguments must be variables) 
SOMEF(A(1I), B)=A(1I)/Bt3. (Arguments must be nonsubscripted) 
SUBPROGRAM (A, B)=A*¥*24+B**2 (Function name exceeds limit of Six 
| characters) 
3FUNC(D)=3 14*E (Function name must begin with an 
alphabetic character) 
ASF (A)=A+E(T) ; (Subscripted variable in the expres- 
| Sion) 
BAD(A, B)=A+B+EAL(C, D) (Recursive definition not permitted) 
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Invalid statement function references (the functions are defined 
as above): 


WRONG = SUM(TAX, FICA) (Number of arguments does not agree 
| } with above definition). 
. MIX = FUNC(I) (Type of argument does not agree with 


above definition) 


FUNCTION SUBPROGRAMS 


The FUNCTION subprogram is a FORTRAN subprogram consisting of a 
FUNCTION statement followed by other statements including at least one 
RETURN statement. It is an independently written program that is 
executed wherever its name is referenced in another program. 
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Where: Type is INTEGER, REAL, DOUBLE PRECISION, COMPLEX, or LOSIC- 
AL. Its inclusion is optional. 


mame is the name of the FUNCTION. 


*s represents one of the permissible length specifications 
for its associated type. It may be included optionally only 
when Type is specified. It must not be used when DOUBLE 
PRECISION is specified. 


Age Gag Aay,eeegan are dummy arguments. They must be nonsub- 
Scripted variable, array, or dummy names of SUBROUTINE or 
Other FUNCTION subprograms. (There must be at least one 
argument in the argument list.) 
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A type declaration for a function name may be made by the predefined 
convention, by an IMPLICIT statement, by an explicit specification in 
the FUNCTION statement, or by an explicit specification statement within 
the FUNCTION subprogram. 


A type declaration in the FUNCTION statement Overrides any type 
declaration by an explicit specification statement within the FUNCTION 
Subprogram. The function must also be typed in the calling program, if 
the predefined convention is not used. 


Since the FUNCTION is a separate subprogram, the variables and state- 
ment numbers within it do not relate to any other program. 


The FUNCTION statement must be the first statement in the subprogram. 
The FUNCTION subprogram may contain any FORTRAN statement except a _ 
SUBROUTINE statement, another FUNCTION statement, or a BLOCK DATA state- 
ment. If an IMPLICIT statement is used in a FUNCTION subprogram, it 
must immediately follow the FUNCTION statement. 


The name of the function must be assigned a value at least once in 
the subprogram -- either as the variable name to the left of the equal 
Sign in an assignment statement, as an argument of a CALL statement, or 
in an input list (READ statement) within the subprogram. 
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The dummy arguments of the FUNCTION subprogram 
(isCeyg Ary Any Azy «¢,/An) may be considered to be dummy variable names. 
These are replaced at the time of execution by the actual arguments sup- 
plied in the function reference in the calling program Additional 
information about arguments is in the section “Arguments in a FUNCTION 
Or SUBROUTINE Subprogran. “ 


The relationship between variable names used as arguments in the cal- 
ling program and the dummy variables used as arguments in the FUNCTION 
subprogram is illustrated in the following example: 


Example 1: 
Calling Program FUNCTION Subprogram 
FUNCTION CALC (A,B,J) 
- I = J*¥2 
ANS = ROOT1*CALC(X, Y, 1) ‘ 
- CALC = A**T/B 
RETURN 
END 


Explanation: 


In this example, the values of X, Y, and I are used in the FUNCTION 
subprogram as the values of A, B, and J, respectively. The value of 
CALC is computed, and this value is returned to the calling program 
where the value of ANS is computed. The variable I in the argument list 
of CALC in the calling program is not the same as the variable I 
appearing in the subprogran. 


Example 2: 
Calling Program FUNCTION Subprogram 
INTEGER*2 CALC INTEGER FUNCTION CALC*¥2(I,J,K) 
ANS=ROOT1*CALC(N,M, P) CALC = ItJd+K*¥¥*2 


RETURN 
END 


Explanation: 


The FUNCTION subprogram CALC is declared as type INTEGER of length 2. 
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RETURN and END Statements in a FUNCTION Subprogram 


All FUNCTION subprograms must contain an END statement and at least 
one RETURN statement. The END statement specifies, for the compiler, 
the end of the subprogram; the RETURN statement signifies a logical con- 
clusion of the computation and returns the computed value and control to 
the calling program. There may be more than one RETURN statement in a 
FORTRAN subprogran. 


Example: 


FUNCTION DAV (D,E,F) 
10 A = Dt+2.0*E 


20 DAV = A+B*¥*2 


RETURN 
30 DAV = B*¥*2 


RETURN 
END 


SUBROUTINE SUBPROGRAMS 


The SUBROUTINE subprogram is similar to the FUNCTION subprogram in 
many respects. The rules for naming FUNCTION and SUBROUTINE subprograms 
are Similar. They both require an END statement, and they both contain 
the same sort of dummy arguments. Like the FUNCTION subprogram, the 
SUBROUTINE subprogram is a set of commonly used computations, but it 
need not return any results to the calling program, as does the FUNCTION 
subprogram. 


The SUBROUTINE subprogram is referenced by the CALL statement, which 
consists of the word CALL followed by the name of the subprogram and its 
parenthesized arguments. | 
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RETURN 


END 


Where: name is the SUBROUTINE name (see the section "Naming 
Subprograms"). 


A4¢a24,aA3eeee,an Are dummy arguments. (There need not be 
any.) Each argument used must be a nonsubscripted variable 
Or array name, the dummy name of another SUBROUTINE or 
FUNCTION subprogram, or of the form * where the character 
“*" denotes a return point specified by a statement number 
in the calling program. 


Since the SUBROUTINE is a separate program, the variables and state- 
ment numbers within it do not relate to any other program. 


The SUBROUTINE statement must be the first statement in the subpro- 
gram. The SUBROUTINE subprogram may contain any FORTRAN statement 
except a FUNCTION statement, another SUBROUTINE statement, or a BLOCK 
DATA statement. If an IMPLICIT statement is used in a SUBROUTINE sub- 
program, it must immediately follow the SUBROUTINE statement. 


The SUBROUTINE subprogram may use one or more of its arguments to 
return values to the calling program. Any arguments so used must appear 
to the left of an arithmetic statement in an input list within the sub- 
program, as arguments of a CALL statement, or as arguments in a function 
reference. The SUBROUTINE name must not appear in any other statement 
in the SUBROUTINE subprogram. 


The dummy arguments (a4, aay Azyeee,an) May be considered dummy vari- 
able names that are replaced at the time of execution by the actual 
arguments supplied in the CALL statement. Additional information about 
dummy arguments is in the section “Arguments in a FUNCTION or SUBROUTINE 
Subprogram." 


Example: The relationship between variable names used as arguments in 
the calling program and the dummy variable used as arguments in the SUB- 
ROUTINE subprogram is illustrated in the following example. The object 
of the subprogram is to "copy" one array directly into another. 


Calling Program SUBROUTINE Subprogram 
DIMENSION X(100),Y(100) 

‘ SUBROUTINE COPY(A,B,N) 

é DIMENSION A (100),B(100) 
K = 100 bdDO 10 I= 1, N 
CALL COPY (xX, Y,K) 10 BCI) =A (I) 

RETURN 
. END 
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CALL Statement 


The CALL statement is used to call a SUBROUTINE subprogram. 


Ger as aes ae fe a, se a ethane Fe Oe en gem ene ey ER ee ee ee ee 1 
| General Form | 
fare a a ee a er ee ee | 
] CALL name (a1,a2;a3¢+++,san) | 
| 
| | | 
|] Where: name is the name of a SUBROUTINE subprogram. | 
| 
| . | 
] A4¢42a~¢aA3,ee0e,an are the actual arguments that are being sup- | 
] plied to the SUBROUTINE subprogram. Each may be of the form | 
] én where n is a statement number (see “RETURN Statements in | 
a SUBROUTINE Subprogram"). | 
Da a a Nr ai a ak a et ce J 
Examples: 

CALL OUT 

CALL MATMPY (X,5,40,Y,7, 2) 

CALL QDRTIC (X,Y, Z,ROOT1, RCOT2) 

CALL SUB1(X+Y*5,ABDF, SINE) 

The CALL statement transfers control to the SUBROUTINE subprogran, 
and replaces the dummy variables with the value of the actual arguments 
that appear in the CALL statement. 

RETURN Statements in a SUBROUTINE Subprogram 
ae ee aN grrr Be yrs ee ny cape pe TG PTY ay Tog eg ee a ee ee 1 
| General Form | 
Sg ei A net Sa I ee ce ee he aide tk a ie i eee 4 
| RETURN | 
| | 
| RETURN i | 
| | 
] Where: i is an integer constant or variable of length 4 whose | 
| value, say n, denotes the nth statement number in the argu- | 
| ment list of a SUBROUTINE statement; i may be specified only | 
] in a SUBROUTINE subprogram. | 
a ee ee eee Be eee J 


The normal sequence of execution following the RETURN statement of a 
SUBROUTINE subprogram is to the next statement following the CALL in the 
calling program. It is also possible to return to any numbered state- 
ment in the calling program by using a return of the type RETURN i. 
Returns of the type RETURN may be made in either a SUBROUTINE or 
FUNCTION subprogram (see "RETURN and END Statements in a FUNCTION Sub- 
program"). Returns of the type RETURN i may only be made ina 
SUBROUTINE subprogram. In a main program, a RETURN statement performs 
the same function as a STOP statement 
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Example: 


Calling Program Subprogram 
‘ SUBROUTINE SUB (X,Y, Z, *, *) 


10 CALL SUB (A,B,C, &30, &40) 


20 Y=A+B 100 IF (M) 200,300,400 
° 200 RETURN 
e 300 RETURN 1 
. 400 RETURN 2 

30 Y=A+C END 


Explanation: 


In the preceding example, execution of statement 10 in the calling 
program causes entry into subprogram SUB. When statement 100 is 
executed, the return to the calling program will be to statement 20, 30, 
or 40, if Mis less than, equal to, or greater than zero, respectively. 


A CALL statement that uses a RETURN i form may be best understood by 
comparing it to a CALL and computed GO TO statement in sequence. For 
example, the following CALL statement: 


CALL SUB (P, &20,0, &35,R, &22) 
is equivalent to: 


CALL SUB (P,Q,R,I) 
GO TO (20,35,22),1 


where the index I is assigned a value of 1, 2, or 3 in the called 
subprogram. 


ARGUMENTS IN A FUNCTION OR SUBROUTINE SUBPROGRAM 


The dummy arguments of a subprogram appear after the FUNCTION or 
SUBROUTINE name and are enclosed in parentheses. They are replaced at 
the time of execution by the actual arguments supplied in the CALL state- 
ment oor function reference in the calling program. The dummy argu- 
ments must correspond in number, order, type, and length to the actual 
arguments. For example, if an actual argument is an integer constant, 
then the corresponding dummy argument must be an integer of length 4. 
If a dummy argument is an array, the corresponding actual argument must 
be (1) an array, or (2) an array element. In the first instance, the 
Size of the dummy array must not exceed the size of the actual array. 
In the second, the size of the dummy array must not exceed the size of. 
that portion of the actual array which follows and includes the desig- 
nated element. 
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The actual arguments can be: 
Any type of constant except hexadecimal 


Any type of subscripted or nonsubscripted variable (except one last 
defined by an ASSIGN statement) 


e An array name 
e An arithmetic or logical expression 


e The name of a FUNCTION or SUBROUTINE subprogram 


A statement number (for a SUBROUTINE subprogram only, see the sec- 
tion “RETURN Statements in a SUBROUTINE Subprogram") 


If a literal constant is passed as an argument, the actual argument 
passed is the literal as defined, without delimiting apostrophes or the 
preceding wH specification. An actual argument which is the name of a 
subprogram must be identified by an EXTERNAL statement containing that 
name. 


When the dummy argument is an array name, an appropriate DIMENSION or 
explicit specification statement must appear in the subprogram. None of 
the dummy arguments may appear in an EQUIVALENCE or COMMON statement. 


If a dummy argument is assigned a value in the subprogram, the corre- 
sponding actual argument must be a subscripted or unsubscripted variable 
name, Or an array name. A constant should not be specified as an actual 
argument unless the programmer is certain that the corresponding dummy 
argument is not assigned a value in the subprogram. 


A referenced subprogram cannot define dummy arguments such that the 
subprogram reference causes those arguments to be associated with other 
dummy arguments within the subprogram or with variables in COMMON. For 
example, if the function DERIV is defined as 


FUNCTION DERIV (X,Y, 2) 
COMMON W 


and if the following statements are included in the calling program 


COMMON B 


C = DERIV (A,B,A) 


then X, Y, Z, and W cannot be defined (e.g., cannot appear to the left 
of an equal sign in an arithmetic statement) in the function DERIV. 


Enclosing a dummy argument in slashes (e.g.e, (/X/,/Y/)) ensures that 
it will be referred to in the subprogram by location. In reference by 
location, the subprogram reserves no storage for the dummy argument. 

The subprogram uses the location of the corresponding actual argument 
for its calculations. This contrasts with reference by value in which 
the dummy argument is assigned storage in the subprogram. When the sub- 
program is entered, the value of the actual argument is brought in from 
the calling program. When the subprogram terminates, the final value is 
transmitted back to the actual argument. An argument which is an ex- 
pression (other than a simple variable) is always referred to by value 
unless its corresponding dummy argument is enclosed in slashes. 
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MULTIPLE ENTRY INTO A SUBPROGRAM 


The standard (normal) entry into a SUBROUTINE subprogram from the 
calling program is made by a CALL statement that references the subpro- 
gram name. The standard entry into a FUNCTICN subprogram is made by a 
function reference in an arithmetic expression Entry is made at the 
first executable statement following the SUBROUTINE or FUNCTION 
statement. 


It is also possible to enter a subprogram (either SUBROUTINE or 
FUNCTION) by a CALL statement or a function reference that references an 
ENTRY statement in the subprogram. Entry is made at the first execut- 
able statement following the ENTRY statement. 


EP EES EIN EES ERS AAS AES ETD GY EES NTS CEE EEE ED GEE GEE IE REDE SEY GEES CES IES CE cee aE OEE te ee ED RIED RD I CENTS enDD UaRD ce CED ED Wey CU ES ARERD GAREY Gite CUBS AOE commen Cay QED EY ERS GED UND EEE eUUEY GEES CGE GEE GP GES Gi EP EE ag ae tee ae ee 


Where: name is the name of an entry point (see the section “Naming 
Subprograms"). 


A1¢82,aAageeeean are the dummy arguments corresponding to an 
actual argument in a CALL statement or in a function 
reference. 


<a ne eee ce es es es es es es ee ee es es ee es es es ee es es es ee as es es ee ee ee es a ce ee es es es ees es ee es es es es ce et ee a ee ee es es es a ee es ee a a es es a es ee es oe 


ENTRY statements are nonexecutable and do not affect control sequenc- 
ing during execution of a subprogram. A subprogram must not reference 
itself directly or through any of its entry points. Entry cannot be 
made into the range of a DO. The appearance of an ENTRY statement does 
not alter the rule that statement functions in subprograms must precede 
the first executable statement of the subprogranm. 


The dummy arguments in the ENTRY statement need not agree in order, 
type, or number with the dummy arguments in the SUBROUTINE or FUNCTION 
statement or any cther ENTRY statement in the subprogram. However, the 
arguments for each CALL or function reference must agree in order, type, 
and number with the dummy arguments in the SUBROUTINE, FUNCTION, or 
ENTRY statement that it references. 


Entry into a sukprogram initializes the dummy arguments of the 
referenced ENTRY statement. Thus, all appearances of these arguments in 
the whole subprogram are initialized. Arguments that were referenced by 
value at some previous use of the subprogram need not appear in the 
argument list of the ENTRY statement. In this case, the reference will 
not transmit new values for the arguments not listed. A function 
reference, and hence the corresponding ENTRY statement, must have at 
least one argument. 


If a dummy argument is listed at more than one entry point, it must 
be consistently referenced either by name or by value. A dummy argument 
must not be used in any executable statement in the subprogram unless it 
has been previously defined as a dummy argument in an ENTRY, SUBROUTINE, 
or FUNCTION statement. 
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If information for an object-time dimension array is passed in a 


reference to an ENTRY statement, 


the array name and all of its dimension 


parameters (except any that are in a common area) must appear in the 
argument list of the ENTRY statement. 


In a FUNCTION subprogram, the types of the function name and entry 
name are determined by the FUNCTION and ENTRY statements. The types of 
these variables (i.e., the function name and entry names) can be dif- 
ferent; the variables are treated as if they were equivalenced. After 
one of these variakles is assigned a value in the subprogram, the others 


become indeterminate in value. 


Upon exit from a FUNCTION subprogram, the value returned is the value 
last assigned to the function name or any entry name. It is returned as 
though it were assigned to the name in the current function reference. 
If the last value is assigned to a different entry name, and that entry 
name differs in type from the name in the current function reference, 
the value of the function is undefined. 


Example 1: 


Calling Program 


TARLE(1) = FUNC(M, X, Y,Z) 
DO 5 I=2,100 
TABLE(I) = ENT(U) 


5 CONTINUE 


Subprogram 
FUNCTION FUNC (T,A,8,C) 


ENTRY ENT(T) 


FUNC = A * B+ C ¥** T 
RETURN 


END 


Explanation: The FUNCTION subprogram is entered once at entry point 
FUNC and initial values are assigned to the dummy arguments T, A, B, and 
C. Thereafter, the FUNCTION subprogram is entered at entry point ENT, 


and only one value is transmitted. 


No new values are passed for A, B, 


or C, so their values are changed only by operations in the subprogram. 
(Note that the original reference to A, B, and C must be by value -- not 


a reference by location. ) 


Each time, the result of the FUNCTION subprogram is returned to the 
main program function reference by the variable FUNC. If FUNC and ENT 


had been of different types, 


it would have been necessary to have 


returned the result by FUNC the first time and by ENT the rest of the 


times. 


94 


Example 2: 


Calling Program Subprogram 
r SUBROUTINE SUB1 (U,V,W,X,Y,2Z) 
° RETURN 
e ENTRY SUB2 (T, *, *) 
CALL SUB1 (A,E,C,D,E,F) U = V*¥ WtT 
‘ ENTRY SUB3 (*, *) 
° X = Y**Z 
° 50 IF (W) 100, 200, 300 
CALL SUB2(G, &10, &20) 100 RETURN 1 
Y= G 200 RETURN 2 
. 300 RETURN 
CALL SUB3(&10, 620) 
Y = AtB 
10 Y = CtD 
20 YY = EtF 
Explanation: 


In this example, a call to SUB1 merely performs initialization. A 
subsequent call to SUB2 or SUB3 causes execution of a different section 
of the SUB1 subroutine. Then, depending upon the result of the arith- 
metic IF statement at statement 50, control returns to the calling pro- 
gram at statement 10, 20, or the statement following the call. 


EXTERNAL STATEMENT 


a a a a a a a 1 
| General Form | 
a a ee ae | 
| EXTERNAL a, byCyece | 
| 
] Where: a,byC,e«- are names of subprograms that are passed as argu- ]j 
{ ments to other subprograms. | 
a na is oe Se ee ge re ey J 


The EXTERNAL statement is a specification statement, and must precede 
statement function definitions and the executable statements. 


If the name of a FORTRAN supplied in-line function is used in an 
EXTERNAL statement, the function is not expanded in-line when it appears 
as a function reference. Instead, it is assumed that the function is 
supplied by the user or is part of the FORTRAN-supplied library. (The 
FORTRAN-supplied in-line and out-of-line functions are given in Appendix 
Co) = | of 


The name of any subprogram that is passed as an argument to another | 
subprogram must appear in an EXTERNAL statement in the calling program. 
For example, assume’ that SUB and MULT are subprogram names in the fol- 
lowing statements: ; | 7 
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Calling Program Subprogram 


SUBROUTINE SUB(K, Y, Z) 


: IF (K) 4,6,6 

; 4 D= Y (kK, Z**2) 
EXTERNAL MULT 
CALL SUB (J, MULT,C) 6 RETURN 

; END 


Explanation: 


In this example, the subprogram name MULT is used as an argument in 
the subprogram SUR. The subprogram name MULT is passed to the dummy 
variable Y as are the variables J and C passed to the dummy variables kK 
and Z, respectively. The subprogram MULT is called and executed only if 
the value of K is negative. 


Example 2: 
e SUBROUTINE SUB (W,X, Y, Z) 
CALL SUB (A,B,MULT (C,D), 37) — 
e RETURN 
° END 
Explanation: 


In this example, an EXTERNAL statement is not required because the 
subprogram named MULT is not an argument; it is executed first and the 
result becomes the argument. 


OBJECT-TIME DIMENSICNS 


If an array is used in a FUNCTION or SUBROUTINE subprogram and its 
name is not in a COMMON statement within the subprogram, the absolute 
dimensions of the array do not have to be explicitly declared in the 
subprogram by constants. Instead, an explicit specification statement 
or DIMENSION statement appearing in the subprogram may contain integer 
variables of length 4 to specify the size of the array. When the sub- 
program is called, these integer variables receive their values from the 
actual arguments in the calling program reference or from common. Thus, 
the dimensions of a dummy array appearing in a subprogram may change 
each time the subprogram is called. 
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The absolute dimensions of an array must be declared in the calling 
program Or in a higher level calling program, and the array name must be 
passed to the subprogram in the argument list of the calling program. 
The dimensions passed to the subprogram must be less than or equal to 
the absolute dimensions of the array declared in the calling program. 
The variable dimension size can be passed through more than one level of 
Subprogram (1.e., to a subprogram that calls another subprogram, passing 
it dimension information). 


Integer variables in the explicit specification or DIMENSION state- 
ment that provide dimension information must not be redefined within the 
subprogram; i.e., they must not appear to the left of an equal sign. 


The name of an array with object-time dimensions cannot appear ina 
COMMON statement. 


Example 1: 


° SUBROUTINE SUBR1(...R,L,M...) 


DIMENSION A(5,10) ee. : 
. REAL... R(L,M)... 


CALL SUBR1i(...A,5,10...) 


: DO 10 J=1,M 
; 10 R(I,J)=0. 
END : 
RETURN 
END 


Explanation: 


This example shows the use of object-time dimensions to supply 
dimension information to a subroutine that will perform some operation 
on an array Of any specified size. In this case, the dimensions passed 
are those specified for the array in the calling program, i.e., the full 
Size of the array. 
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Example 2: 


DIMENSION A(5,10) 


° 10 


SUBR1 (0 eA, ly Jeee) 


Explanation: 


SUBROUTINE SUBR1i( ..R,L,M...) 


REALe. R(LyM) e. 


DO 10 I=1,L 
DO 10 J=1,M 
R(I,J)=0. 


RETURN 


END 


This example shows the use of object-time dimensions to specify a 


subset of the extent of an array to a subprogram. 


The effect of this 


coding is the same as if another array, B, of dimensions (4,7) had been 
defined in the calling program and had been made equivalent to array A; 
the array B and its dimensions would then have been passed to SUBR1 as 


follows: 


DIMENSION..-A(5,10) ,B(4,7) «4. 


EQUIVALENCE (A(1,1),B(1,1)) 


CALL SUBR1 (...B,I,Je..) 


END 


98 


BLOCK DATA SUBPROGRAM 


To initialize variables in a labeled (named) common block, a separate 
subprogram must be written. This separate subprogram contains only the 
DATA, COMMON, DIMENSION, EQUIVALENCE, and Type statements associated 
with the data being defined. Data may not be initialized in unlabeled 
common. 


CS er NG TO ee OCT pe Ee Teele Fe A TOE TT ae Choe TOLD ae OT ee 1 
| General Form | 
|---------------------------------------------------------------------- | 
| BLOCK DATA | 
| | 
| | 
| 
| END | 
Oe pace ema eee ar PRAT TP RE NCEE eo en in eee eR a a ana oe A a RE ROE eee TEN J 


1. The BLOCK DATA subprogram may not contain any executable 
statements. 


2. The BLOCK DATA statement must be the first statement in the subpro- 
gram. If an IMPLICIT statement is used in a BLOCK DATA subprogram, 
it must immediately follow the BLOCK DATA statement. The COMMON 
statement must precede the data initialization statements. 


3. Any main program or subprogram using a common block must contain a 
COMMON statement defining that block. If initial values are to be 
assigned, a BLOCK DATA subprogram is necessary. 


4. All elements of a common block must be listed in the COMMON state- 
ment, even though they are not all initialized; for example, the 
variable A in the COMMON statement in the following example does 
not appear in the data initialization statement: 


BLOCK DATA 

COMMON/ELN/C, A, B/RMG/Z, Y 

REAL B(4)/1.40, Le 2, 2*1 6 3/, Z2*8 (3) /3*7. 64980825D0/ 
COMPLEX C/(2.4,3.769)/ 

END 


5. Data may be entered into more than one common block in a single 
BLOCK DATA subprogram. 


6. Only one BLOCK DATA subprogram may be used to enter data into a 
particular common block. 
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APPENDIX A: SOURCE PROGRAM CHARACTERS 


Peet ge pe ne ee, eee ger eae Me er ee re re 
] Alphabetic Characters | Numeric Characters | 
}---------------------------------- }----------------------------------- { 
A | 0 | 
| B | 1 | 
Cc | 2 | 
| D | 3 | 
| E | 4 | 
| F | 5 | 
G | 6 | 
| H | 7 | 
| I | 8 | 
| J | 3 | 
| K | | 
| L | | 
| M }----------------------------------- { 
| N | Special Characters | 
| O a a ga { 
| P | (blank) | 
| Q | ‘ | 
| R | = | 
| S | / | 
| T | = | 
U | . | 
| V | ) | 
W | * | 
| X | ‘ | 
| ¥ | ( | 
i Z | © (apostrophe) | 
| $ | & | 
ESR Se Pera bon Ee Pe ng RMR Vea ree EEE So Te Ne ate eat era ee es eee So eee ee 4 


The 49 characters listed above constitute the set of characters 
acceptable by FORTRAN, except in literal data where any valid card code 
is acceptable. 


100 


APPENDIX B: OTHER FORTRAN STATEMENTS ACCEPTED BY FORTRAN IV 


This appendix discusses those features of previously implemented 
FORTRAN IV languages that are incorporated into the System/360 FORTRAN 
IV language. The inclusion of these additional language facilities 
allows existing FORTRAN programs to be recompiled for use on the IBM 
System/360 with little or no reprogramming. 


READ STATEMENT 


Where: b is the statement number or array name of the FORMAT state- 
ment describing the data. 


list is a series of variable or array names, separated by 
commas, which may be indexed and incremented. They specify 
the number of items to be read and the locations in storage 
into which the data is placed. 


AO SE EE OE aS a SY EY US RE EE ES ee ES ee ee ee ce ce ces ee ee ce ee ce ee ce ce ee ee ce es ee se ee ee es ee ee te ee ce ee ce ee ee ee ee ce es ee ee ee ee ee ee ee ee ee ee ee ee 


This statement has the effect of a READ (n,b) list statement where b 
and list are defined as above, and the value of n is installation 
dependent. 


PUNCH STATEMENT 


i Sg A SA ca A aa aaa ka a aa 1 
] General Form | 
Sah a ec a i a a a ee 4 
] PUNCH b, list | 
| | 
| 
| | 
] Where: b is the statement number or array name of the FORMAT state- | 
| ment describing the data. | 
| | 
| 
a | 
| list 1s a series of variable or array names, separated by | 
| commas, which may be indexed and incremented. They specify | 
] the number of items to be written and the locations in 
| storage from which the data is taken. | 
UE Reeder eae eet Pee ee a eee trl PE See REIT EO CY A ON ERO RSA SORES PRO SON Ne eae a SRO Tee ee J 

This statement has the effect of a WRITE (n,b) list statement where b 


and list are defined as above, and the value of n is installation 
dependent. 
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PRINT STATEMENT 


By a ee ee ay oe eae Ge ee ey GE a ye et eae ge Se ey rg 1 
| General Form | 
}~~-~~---~-~~~~------~--------~--~--~-----~-------~-------+---+-------- { 
] PRINT b, list | 
l | 
| | 

| | 
| Where: b is the statement number or array name of the FORMAT state- |] 
] ment describing the data. | 
| | 
aes | , | 
| list is a series of variable or array names, separated by ] 
] commas which may be indexed and incremented. They specify | 
| the number of items to be written and the locations in | 
! storage from which the data is taken. 
a a es ee ee J 


This statement has the effect of a WRITE (n,b) list staterent where b 
and list are defined as above, and the value of n is installation 
dependent. 
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APPENDIX C:  FORTRAN-SUPPLIED SUBPROGRAMS 


The FORTRAN-supplied subprograms are of two types: mathematical sub- 
programs and service subprograms. The mathematical subprograms corres- 
pond to a FUNCTION subprogram; the service subprograms correspond to a 
SUBROUTINE subprogram. An in-line subprogram is inserted by the FORTRAN 
compiler at any point in the program where the function is referenced. 
An out-of-line subprogram is located ina library and the compiler 
generates an external reference to it. A detailed description of out- 
of-line mathematical subprograms and service subprograms is given in the 
publication IBM System/360: FORTRAN IV Library Subprograms, Form 
C28-6596. Table 4 shows mathematical function subprograms, and Table 5 
shows out-of-line service subprograms. An asterisk in the in-line/ 
out-of-line column of Table 4 identifies the function as one defined as 
an intrinsic function in USA FORTRAN IV. 
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eee eae ee EUHE eRe ERD cD ER ca FORE ERD ER a AEP I ent EY EY SEY ERIS SOOT AAS SET RED ED CES SE SAE SE ET HY RY ONE SID ND SEED SARIS END UES cn ONE etmeS eS ceE “PRIO Ome SE cine SEER itt EERO Ge GEE eee ERED ACE SRN IES SEES ee ERY ee SE <I A ee GY cee a eS ee SY eet eine OEY oD Go ime EY cee Si ce Gn eee Oe 


104 


Table 4. Mathematical Function Subprograms (Part 1 of 3) 
(nas er ie ee ee eee te toe : a ai See te ee eo 1 
| | . | | |Type of | 
| [Entry | In-Line (I) |No. of| Type of | Function | 
] Function |Name_ | Definition jOut-of-Line (0O)] Arg. | Arguments [Value | 
}--------------4------ {--~---------------- $---------------4------}~----------}----------- { 
| Exponential | EXP fears | O | 1 |Real *4 }Real *4 | 
] {|DEXP | ea@rg | O | 1 {Real *8 {Real ¥*8 | 
| JCEXP | e@rg | O ] 1  |Complex *8 |Complex *8 | 
| {|CDEXP | e@rdJ { O i 1 |Complex *16|Complex *16| 
}~-------------}------}------------------ }--------------- $~---~--}-----------}----------- { 
{Natural JALOG |ln (Arg) ] Oo | 1 |Real *4 JReal ¥*4 | 
|] Logarithm }|DLOG fin (Arg) | O { 1 {fReal *8 [Real *8 | 
| |CLOG [ln (Arg) | O i 1 |Complex *8 |Complex *8 | 
| {CDLOG {ln (Arg) | Oo | 1 |Complex *16|Complex *16| 
}—------------ +------ }------------------ }--------------- $------}-----------}----------- { 
| Common ] ALOG10 | 10940 (Arg) | O | 1 [Real *4 |Real *4 | 
] Logarithm |] DLOG10] 10g, 0 (Arg) | Oo i 1 {Real *8 [Real *8 | 
~-------------}------}--------=--------- $--------------- 4------4-----------}----------- { 
|Arcsine JARSIN Jarcsin (Arg) ] Oo } 1 |Real *4 [Real *4 | 
| {DARSIN{arcsin (Arg) | O | 1 |Real *8 {Real *8 | 
}~-------------}------}------------------ {--------------- 4------}-----------4----------- { 
|] Arccosine JARCOS fJarccos (Arg) | @) | 1 {Real *4 {Real *4 | 
{ {DARCOS|arccos (Arg) i Oo | 1 |Real *8 {Real *8 
}~---~--------- $------}------~----------- --------------- ------+-----------}----------- { 
JArctangent? JATAN |farctan (Arg) { O | 1 |Real *4 {Real *4 | 
] {ATAN2 jfarctan (Arg,/Arga) | @) | 2 |Real ¥*4 {Real *4 | 
| |DATAN [arctan (Arg) | O | 1 |Real *8 [Real *8 | 
| {|DATAN2]}arctan (Arg,/Arga) | Oo | 2 |Real *8 |Real *8 | 
}--------------+~------ $------------------ --------------- 4------}-----------}----------- { 
|Trigonometric |SIN {sin (Arg) | O | 1 |Real *4 [Real *4 | 
| Sine |]DSIN |sin (Arg) | O { 1 jReal *8 [Real *8 | 
| (Argument in |CSIN [sin (Arg) | O | 1 [Complex *8 [Complex *8 | 
{| radians) {CDSIN Jsin (Arg) | O i 1  |Complex *16|Complex *16| 
~-------------}------}-------=---------- }~--------------4------4-----------4-----------} 
|Trigonometric [COS jcos (Arg) | O | 1 |Real ¥*4 jReal *4 | 
{Cosine }DCOS Jcos (Arg) | Oo | 1 jJReal *8 [Real ¥*8 | 
{(Argument in |CCOS |cos (Arg) ] O | 1 |Complex *8 |Complex *8 | 
| radians) {CDCOS [cos (Arg) | O | 1 |Complex *16|Complex *16| 
~--~-~~---~~--}------}------------------ 4---~-----------4~-----}-----------}-----------] 
|Trigonometric |TAN jtan (Arg) | O | 1 |Real *4 |Real *4 | 
[Tangent | | | | | 
} (Argument in ]DTAN |tan (Arg) | O ] 1 |Real *8 [Real ¥*8 | 
{] radians) | ] | | | 
}~----~-------- $------ {------------------ }--------------- 4------}-----------}----------- { 
|Trigonometric ]COTAN |cotan (Arg) ] Oo | 1 |Real ¥*4 [Real *4 | 
| Cotangent | | | | | 
| (Argument in |DCOTAN|cotan (Arg) | O | 1 |Real *8 [Real *8 | 
| radians) | | | | | | 
—~--~---------}------ femme Fp nn nnn enn nnn nnn penn nfm nnn panna ncnnnnd 
|Square Root |] SORT [<Arg) | O | 1 |Real ¥*4 [Real *4 | 
| |DSORT | (Arg) { O | 1 jReal *8 {Real *8 | 
| | CSQRT | (arg)? | O | 1 |Complex *8 |Complex *8 | 
{cpsort| (arg)”” O 1 1 [complex *16 | Complex #16 | 
{eee ee el ee fee ee aerate oer ae ee i Meera a Ma a eae | 
| |4Two ae must be supplied for ATAN2 and DATAN2. | 
J 


Table 4. Mathematical Function Subprograms (Part 2 of 3) 
(oe ei ee er Yee er Top ee 5 aan mae y Oana areas atm aha To ee 1 
| | | | | | [Type of —| 
| JEntry | | In-Line (I) JNo. of] Type of {Function | 
| Function |Name_ | Definition ]Out-of-Line (0O)] Arg. | Arguments |Value | 
--------------}------ }------------------4---------------4------}---=-------4-----------4 
{Hyperbolic {TANH |tanh (Arg) } O | 1 ]Real *4 jReal *4 | 
|} Tangent {DTANH |tanh (Arg) | O | 1 |Real *8 |Rea™ *8 | 
~---+----------4----~-}------------------ 4---------------4------}-----------4-----------4 
JHyperbolic {SINH |sinh (Arg) | O | 1 jReal *4 {Real ¥*4 | 
|Sine ]DSINH |sinh (Arg) | O | 1 |Real ¥*8 {Real *8 | 
~-------------4------}------------------ 4------------~--4------4-----------4-----------4 
|Hyperbolic |COSH Jcosh (Arg) \ O | 1 |Real *4 jReal ¥*4 | 
| Cosine |DCOSH |cosh (Arg) | O j 1 |Real *8 |Real *8 | 
Soa a SL ESS OR aaa aaa aa a aaa fe eS ee eee eae ees { 
J}Error Function|ERF | 2x -—u? | O | 1 jJReal ¥*4 |Real *4 | 
| | DERF l ir ( e du | O | 1 |Real *8 |Real *8 | 
| i 0 | | | | 
~-~-----------}------ }------------------}---------------4------4-----------4-----------J 
|Complemented |ERFC |i-erf (x) | O | 1 -|Real *4 }Real ¥*4 | 
{Error Function|DERFC |} | O | 1 jReal *8 |Real ¥*8 | 
}~-------------+------ }------------------ 4--------------- }------ -----------}----------- { 
| Gamma |GAMMA | ¢%® x-1 -u | O | 1 {Real *4 jReal *4 | 
| | DGAMMA | \ u e du ] O | 1 jReal *8 |Real *8 | 
| | | v0 | | | | | 
}-------------- }------ }------------------ 4--------------- 4~-----}----------- 4----------- { 
|] Log-gamma | ALGAMA | | O ] 1 |Real ¥*4 {Real *4 | 
| | DLGAMA| 1°de6 ee O 1 |Real *8 |Real ¥*8 | 
}-------------- +------ $------------------ }--------------- 4------ $-----------}----------- { 
|Modular {MOD JArg4 (mod Argsz)= $ tL * | 2 jInteger *4 |Integer *4 | 
| Arithmetic? JAMOD {Arg4,-([x]*Argz2 | ft | 2 |Real ¥*4 |Real*4 | 
| |DMOD |Where: [x] is the | I * | 2 |Real *8 {Real *8 | 
| Jlargest integer | | | | | 
| | jwhose magnitude | | | | | 
i }does not exceed | | | | | 
| | Jthe magnitude of | | | | | 
| | |Argi/Arga-e The | | | | 
| i {sign of the inte- | | | | | 
| | |ger is the same | | | | | 
i | Jas the sign of | | | | | 
| |Args/Arga. | | | | | 
}--------------}------}------------------ }--------------- }------ 4----------- 4---~------- { 
|Absolute value]IABS- |{Arg] | I * | 1 fInteger *4 |Integer *4 | 
| J|ABS } | I * | 1 |Real ¥*4 {Real ¥*4 | 
| {DABS | | I * | 1 |Real *8 {Real *8 | 
}-=---------------- 4--------------- 4------}-----------4----------- { 
| |CABS | ((a2+b2) for atbi | O | 1 |Complex *8 |Real *4 | 
| {CDABS | } O | 1 |Complex *16|Real *8 | 
}--------------4------ }------------------4--------------- }------}-----------}----------- { 
] Truncation | INT ]}Sign of Arg times | | a | 1 |Real *4 {Integer *4 | 
] | jlargest integer | | | | | 
| | [<]Arg| | | | | | 
] JAINT | | L- -* | 1 |Real *4 {Real *4 | 
| J IDINT | | i. st | 1 {Real *8 JInteger *4 | 
}---~----------}------ }------------------ 4--------------- 4------ 4----------- }----------- { 
jLargest value?]AMAX0O |[Max(Arg,,Arga,...-)] O * 1 22 |Integer *4 |Real *4 | 
| ]AMAX1 ] | oO * | >2 {Real *4 {Real ¥*4 | 
| |MAXO | | O * | 22 jInteger *4 |Integer *4 | 
] |MAX1 | | oO * | 22 |Real *4 [Integer *4 | 
| | DMAX1 | | O * | 22 {Real *8 {Real ¥*8 | 
}----~--~-~—-—— BS aaa Qiao SS iS Deter ees ae Dee crit eee Beha Sa Ot 4 

| |}2mMoD and AMOD are not defined when the value of the second argument is zero. 
g 

\?For the FORTRAN IV (H) compiler, these functions are in-line. | 
| }*USA FORTRAN IV intrinsic function | 
D2 a a a ts ce i aS ae ee oe he ee J 
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Table 4. Mathematical Function Subprograms (Part 3 of 3) 
aa eee Sages as aa Re a ac mas ara 5 ania | nea anlar yaa Ca lca 1 
pe | | | | [Type of ~— | 
JEntry | J} In-Line (I) JNo. of] Type of |Function | 
] Function |Name |] Definition jOut-of-Line (0)] Arg. | Arguments |Value | 
~-------------}------ }------------------}---------------}------}----------- $-----------4 
]Smallest 1AMINO |Min(Argy,Arga,«--) | O * | 22 |Integer *4 |Real *4 | 
| values JAMIN1 | | oO * | 22 |Real *4 jReal ¥*4 | 
| |MINO | | oO * | 22 {Integer *4 [Integer *4 | 
| {MINI | | Oo * | >2 JReal *4 [Integer *4 | 
{ |DMINI | | Oo * {| 22 {Real *8 [Real *8 | 
}~-------------}------}------------------ }--------------- }------ }----------- }----------- { 
| Float ] FLOAT |Convert from | I * | 1 |Integer *4 |Real *4 | 
j |] DFLOAT] integer to real | I | 1 [Integer *4 |Real *8 | 
}~------------- -~----}---------~------—- }--------------- ------ }----------- }---~------- { 
| Fix JIFIX ]Convert from ] I * | 1 {Real *4 J|Integer *4 | 
l {HFIX |real to integer } I | 1 jReal *4 {Integer ¥*2 | 
}-------~------}~-----}------------------ {~-------------- 4------ }----------- 4----------- { 
|Transfer of |SIGN J|Sign of Arga times] I * | 2 |Real *4 jReal ¥*4 | 
|sign4 | | JArga | | | | | 
j J ISIGN | [ I * ] 2 |Integer *4 |Integer *4 | 
{|DSIGN } { I * 1 2 |jReal *8 |Real *8 | 
}--------------+------ }------------------ }--------------- ------ }----------- ----------- { 
|] Positive {DIM | Arga-Min(Arg,, | I * | 2 |Real *4 JReal ¥*4 | 
| difference JIDIM {Arg2) | I * | 2 {Integer *4 |Integer *4 | 
~-~-----------}------ }-—-~-------~-~--~-}---------~----- 4 ------f--~-=------f--=--------4 
{Obtaining most|SNGL  ] ] I * ] 1 |Real *8 [Real *4 | 
|significant | | | | | | | 
|part of a | | | | | | 
|Real *8 | | | | | 
|argument | | | | | | 
~=----~-~----- }--~---}------------------4----—---------- 4 ------4----- === f--- === 
JObtain real ]REAL | | I * | 1 |Complex *8 |Real *4 | 
[part of | | | | | | | 
| complex | | | | | | | 
[argument | | | | | | | 
~----~-------- }~~-~--}-~----------------}---~-----------4------ 4 ----------- 4 --- == === 4 
{Obtain }AIMAG | { I * l 1 |Complex *8 |Real *4 | 
| imaginary | | | | 
[part of | | | 
| complex | | | | 
Jargument | i 1 [ | | | 
SSeS et eee a 1 ier cc ar ek Se aaa arias Sek aarees Gina ainsi ae ran Gb anae De saa 
|Express a Real|]DBLE ] ] I * | 1 |Real ¥*4 {Real *8 | 
{*4 argument in| | 1 | | | | 
}Real *8 form | | | 1 [ | | 
}-------------- }------ }------------------ }--------------- 4------ }----------- }----------- { 
{Express two {|CMPLX |] C=Arg,,+iArga | I * | 2 |Real *4 {Complex ¥*8 | 
{real arguments | DCMPLX] | I | 2 jReal *8 {Complex *16| 
[in complex | | | | 
| form | | | | | | 
~~--~---~-----}------}------------------ 4-----~-------~-4------4-----------4----------- 
] Obtain {CONJG {C=X-iY { I * i 1 |Complex *8 |Complex *8 | 
|] conjugate | DCONJG| For Arg=X+iY | I | 1 |Complex #16|Complex *16| 
Jof a complex | | | | | | | 
| argument | | | | | | | 
ease wwe Leelee eee ee ee a ee ae eed eee cnet ee tae eee 
}3For the FORTRAN IV (H) compiler, these functions are in-line. | 
|*SIGN, ISIGN, and DSIGN are not defined when the value of the second argument is zero. | 
{*USA FORTRAN IV intrinsic function | 
Pe a reer SNC Piro enc ye O POR  er SPIN yD TOE ee ESOT OPER CTS EP ae eS a eC on nay CARN SS gO ena een aP OE eee pee SI Ca J 
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Table 5. Out-of-Line Service Subprograms 
(oo ee ee ae ta aaa a cea a Oe a cate ate ee cre ete nee ae ee ce ce EE ee ee ce ee nee A Se RY SOT ee SS Se cee SE nnn OR aD cele ES eee 1 
j Function { CALL Statement | Argument Information | 


}---------------------- $----------------------- fo-2--------------------=---------------- { 


|Alter status of sense |CALL SLITE(i) i is an integer expression. 
jlights 


| 
| 
JIf i 
| turned off. 

{If i = 1, 2, 3, or 4, the corresponding 
| sense light is turned on. 
4 


| 
| 
If i= 0, the four sense lights are | 
| 
| 
| 


ne et ne cy eee cate eee ee ee a ce i eee SE a a we ae ee ees ee ee RO ee ee ce ae Se ce ee ee ne ee ee ee ee Se Co eee eee eee ce cee 


| 
| 
| 
+ 

|Test and record status|CALL SLITET(i, j) iis an integer expression that has a| 
| value of 1, 2, 3, or 4 and indicates | 
| which sense light to test. | 
| j is an integer variable that is set to | 
| | 
| | 
+ 


Jof sense lights | 
| 
} 1 if the sense light was on, or to 2 
| 
4 


if the sense light was off. 


a ae ae eee eee cena a ee ee SY Oe cee SS ee ee ee ee ee ce ee ae ee ee cee ete A cee te ee ee tO ee cr I cee OD ee eS EY SE ee Oe a em OE te ee epee ee ee ee ee 


— 
{ 
! 
{ 
{ 
1 
{ 
! 
! 
i 
{ 
{ 
| 
{ 
{ 
| 
| 
{ 
{ 
t 
t 
\ 
t 


a and b are variables that indicate the] 
limits of storage to be dumped. | 
(Either a or b may be the upper or | 
lower limits of storage, but both must| 
be in the same program Or subprogram | 
or in common.) | 

£ indicates the dump format and may be | 
one of the following: | 

0 - hexadecimal | 

~- logical*1 | 
~ logical*4 | 
| 

| 

| 

| 

| 

| 

| 


]Dump storage on the 
Joutput data set and 
|terminate execution 


CALL DUMP (a,,b1,f,, 
eo0 5ane Dn, En) 


- integer*2 
- integer*4 
real*4 

- real*8 

- complex*8 
~ complex*16 


WOANKDNSF WHE 
! 


- literal 


me se eet et ne es et ee es oe ee a ce es ee 


Se ee eee ee cere cee CY cee ee ee eee ee ee TP ee ee cee oe ee ee ee ee ee ee ee ee eee ee ee ee ee ee ee ee es ee ee 


a, b, and £ are as defined above for| 


ee ee eee ee ee ce em ee ee ey ee re ee ee ee ee ee oe ee es 


CALL PDUMP (a,,b4,f1, 


| Dump storage on the 
joutput data set and 
jcontinue execution 


| 

| | 

| | 

| | 

| | 

| | 

| | 

| | 

| | 

| | 

| | 

| | 

| 

| | 

| | 

| | 

| | 

| | 

+ 4 

| | 

| Be hepa 

}------------------ +---}----------------------- 4---------------------------------------- { 

]Test for divide check }| {j is an integer variable that is set to] 

| exception | | 
| | 
| | 
| | 
+ 4 
| 
| | 
| | 
| 
l | 
| | 
| 
4 + 
| 
L 1 


1 if the divide-check indicator was | 
on, or to 2 if the indicator was off. | 
After testing, the divide-check indi- | 
cator is turned off. | 


ST NS SE NE EP ON TY EP ED EO OEE EE ES OY OE ED eR SY EE a Oe ee epee ene ge ee EE Qe UE ED ee Oe ee ee EE ee ee a eee ee ee ee nee ee eee eee ee tee ene ce SG eee ee EE cee ce ee ee ee eee ee a ee me cy ee a ee ee ee es ee ee eee 


j is an integer variable that is set to] 
1 if an exponent overflow condition | 
was the last to occur, to 2 if no | 
overflow condition exists, or to 3 if | 
an exponent underflow condition was | 
the last to occur. After testing, the| 
overflow indicator is turned off. 


PP A AE ED SS EE EE A ES EY AOE LS ED OS NY ED GDS ED ED ee OD Oe Re Ane ate ED cee Ee ae cee EE ee ee ee ee ee eee ee ee ee ee ae 


|Terminate execution 


{ 
{ 
| 
| 
| 
t 
{ 
| 
{ 
t 
t 
| 
! 
1 
| 
{ 
{ 
{ 
| 
| 
{ 
{ 
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APPENDIX D: SAMPLE PROGRAMS 


SAMPLE PROGRAM 1 


The sample program (Figure 1) is designed to find all of the prime 
numbers between 1 and 1000. A prime number is an integer that cannot be 
evenly divided by any integer except itself and 1. Thus, 7 is a prime 
number. The number 9 is not prime since it can be divided evenly by 3. 


BM FORTRAN Coding Form sic nen 


vos SAMPLE PROGRAM woe ome Te a 
ee om 3764 pe (oe ee a ee a 


STATEMENT IDENTIFICATION 
NUMBER 3 | FORTRAN STATEMENT SEQUENCE 
9 
] 2.3 4 ~=5 16 | 7 8 9 W 112 +13 14 15 16 7 1 19 20 ll 22 a 24 25 26 27 28 29 30 33 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 A? 48 49 50 51 52 53 54 55 56 57 58 59 $0 61 62 63 64 65 66 67 68 oo 70 73 72473 74 75 76 77 78 79 80 


PRIIne| WonBelk leeriee| : 
Ng aan TREE POT T 
























































Bae foued We? leer) 8 


SRAM aCPReCACGHE Ee eee ies it ae SSE TR TE hells bbe 
fps Z| | | tit tte HACER 





















































5 “Eta » hes ay A 2d 


















































ae 





























Sea — TEETER 


[6 [7 8 9 10 13 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 39 3) in 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72]73 74 75 76 77 78 79 80 
*A 2 3 card form, IBM electro 888157, is available for punching statements from this form 





Figure 1 Sample Program 1 
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SAMPLE PROGRAM 2 
The n points (xj, yj) are to be used to fit an m degree polynomial by 
the least-squares method. 
Y = aot a4X + aax? + 226 + AmxM 


In order to obtain the coefficientS ag, a4,-.., am, it iS necessary to 
solve the normal equations: 


(1) Woao + Wyay + «2. + Wmam = Zo 
(2) Wiaqg + Weaq + 22- + Wyyiam = 21 
(m+1) Wmao + Win 4144 + eee + Womam = Zm 
where: 
n 
Wo oe @ | Zo = Zz Yi 
1=1 
n n 
Wy = > X= Ly = 2 Y.xX. 
: 1 ; 
1=1 1=1 ales 
n n 
W2 = 2 x.? Za =~ 2 y:x.? 
1 a 
1=1 =1 = 
Z n 
e 2m > > Vata 
. es 
: 1=1 
n 
Wom = 2 x,2m 
j=1 + 


After the W's and Z's have been computed, the normal equations are 
solved by the method of elimination which is illustrated by the 
following solution of the normal equations for a second degree polynomi- 
al (m = 2). 


(1) Woan + Waa. + Waao = Zo 
(2) W4a9 + Waay + Wa3an = Zy 
(3) Waao + Waa, + W,a2 = Zo 


The forward solution is as follows: 
1. Divide equation (1) by Wo. 


2. Multiply the equation resulting from step 1 by W, and subtract from 
equation (2). 


3. Multiply the equation resulting from step 1 by Ws and subtract from 
equation (3). 
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The resulting equations are: 


(4) Ao + Dyna, + 043an = Day 
C5) ba2a, + bazae = Day 
(6) D328, + b3a3am = Dau 
where: 
Dia = Wi/Wo, Das = W2/Wo, Dag = Zo/Wo 


Doo = Wa-biaWi », bas = Wa-bDisW, , Daw = Za~DauWa 
b32 = W3-bi2Wo , ba3 = Wy-bigW2 , Da, = Za-DiayWo 
Steps 1 and 2 are repeated. using equations (5) and (6), with baz and b3a 
instead of Wo and W,. The resulting equations are: 
(7) Aa, + Ca3am = Cay 
(8) C3342 = C34 
where: 
Ca3 = Das/Daa 1 Cay = Dau/ba2 


C33 = b337-Ca3D32 , Cay = Dau-CauDaa 


The backward solution is as follows: 


(9) Aas = Cay/Ca35 from equation (8) 
(10). “ay. = Gs eCse5, from equation (7) 
(11) ag = byy-by2a,7-~b,3a5 from equation (4) 


Figure 2 is a possible FORTRAN program for carrying out the 
calculations for the case: n= 100, ms 10. Wo, Wa, Wa, .-+, Wam are 
stored in W(1), W(2), W(3), ..., W(2M+1), respectively. Zor “4244 Lar 
«ee, Sm are Stored in Z (1), Z2 (2), Z (3), .«..., Z2(Mt1), respectively. 
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FORTRAN Coding Form 


M 
nocmm SAMPLE PROGRAM 2 Ea CN (a (| 
poe EE FF ff fied —— 


a oa 
pes ioe z IDENTIFICATION 
ae REAL) NIC4g 1% 17:18:19 20 21 23 24 25 26 27 28 29 30 31 32 35 ) 11) s/BI( 6 47 48 49 50 51 52 53 54 55 56 57 70 173 74 75 76 77 78 79 90} 


ETAL] XiC40101)]9 YICIt]00))|i((214)| MN TTT TTT TT 
FORMAT (2913 /(HFAML 7) SHAGAVONGGHSHABAEAAUEGESUOBATGENUAR 
Aa (SIE 51/6 arene ag Se TORT A ISTP Tel 

MIN xc ) BAC Pe eee 

HERES HH 
RS TT ea 
HERERE a ppp ep 
HOaHHHEAHEGALOGLESLLGGIS Ee PETRA E RCRA eae 
OTE SILI ie Leecin e e G GCSCcis ae iio e ls i 
EeePC STE ase al Toh ee ES a aire 
The TEE Se 0 ee ee ee eae riceGe his 

map pee PARR ERME RARER EER RRREE 
EET PEMD ASE Co RE RRe EEE 


ee 
Hesesetasete ase eet apeitaat etaaeee teen aaae 
oth RTE TTT 
adsrerinaeseeee tat ere Hann ee OEE 


[617 8 9 10 U2 131415 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52_53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72473 74 75 76 77 78 79 80 
IBM ‘ailable for stot this for 


SSAA 
v1 eal ax 
Spe 





Az 
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hed 





tery 
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=Eg 
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b= 
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at 
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Figure 2. Sample Program 2 (Part 1 of 3) 
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IBM FORTRAN Coding Form 
recs SAMPLE. PROGRAM 2 | 
' Z INSTRUCTIONS CARO ELECTRO NUMBER™ 
9 10 11 12 15 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 bl 62 63 64 65 66 67 68 69 70 71 72)73 74 75 76 77 78 79 80 
Tiel MOL WGI TPT TTT B a 
pb ia Hike STRRER PEER AERA 
Ho Ka EE we ee eee a ae | 
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Figure 2. Sample Program 2 (Part 2 of 3) 
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FORTRAN Coding Form 




































































fe Printed in U.S.A. 
o 
; crarnc TT [nce 3 or 3 | 
nom SAMPLE PROGRAM 2 mvcane — jctame | | TE mot 3 or 3 
T- : INSTRUCTIONS 
mo = 6/68 es a a 
5 “2 {IDENTIFICATION 
a eeey FORTRAN STATEMENT . SEQUENCE 
3) 3 
} fi 273 "4 5 Je ]7 8 9 10 1 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48-49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72/73 74 75 76 77 78 79 80 
+e) aaa Toole ct inh 7 cs oo oa 
|| | [37] ISITIGMIA) [=| [SITIGMA'+/B(T/-/4/9/J/) AIC J)D 
| TT | tl Is HEBER A Ge RRS RRE ES EEAEe Cele Sei Cele Ee IseL DI aIerer 
PTT TTT ACTD] |=] BT LIBD STIGMA | TT TT TTT TT TTT TET Et tr ET Pe 
PTT mr TEE] CUI) Guia silat) TTT TT TTT TTT TET TTT TTT TT TTT 
Li] aT IRIE! [Celsi2iy| (CAIIT io Ti=|diokzD] LIT TTT TTT TTT TTT EEE 
FEDER Bee CC CHE PEECEE LEE LEELE ECE 
cLeaae : i ! 
Lif2 3 4 s{[6f7_ 8 9 10 11 12 13 14 15 16 17 16.19 20 2 22 23 24 25 26 27 28 29 30 31 32 93 34 35 36 37 98 39--MOral 42 W344 45 46 47 48 49 50 51 52 53 54 55 56 57 5B 59 60 61 62 63 64 65 66 67 68 69 70 71 72|73 74 75 76 77 78'79 80 
5 . > 
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Figure 2. 
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Sample Program 2 (Part 3 of 3) 




























































































X28-7327-5 





























The elements of the W array, except W(1), are set equal to zero. 


W(1) is set equal to 


N. 


For each value of I, 


XI and YI are selected. 


The powers of XI are computed and accumulated in the correct W counters. 
The powers of XI are multiplied by YI, and the products are accumulated 


in the correct Z counters. 


when computing the next power of XI. 


parameters. 


In order to save machine time when the 
object program is being run, the previously computed power of XI is used 


ters have been set as follows: 


W(1) 


W(2) 


W(3) 


W(2M+1) 
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HM @ 


N 
z(1)=2 YI 
I=1 
N 
z(2) = > YIXI 
T=1 
N 
z(3) = d yYIXxXI2 
T=1 
N 
z(M+i) = 2 YIXIM 
I=1 


Note the use of variables as index 
By the time control has passed to statement 17, the coun- 


By the time control has passed to statement 23, the values Of Wo, W4,, 
ee; Wamts have been placed in the storage locations corresponding to 
columns 1 through M+1, rows 1 through Mtl, of the B array, and the 
values Of Zo, Za, «ees Zy have been stored in the locations correspond- 
ing to the column M+2 of the B array. For example, for the illustrative 
problem (M = 2), columns 1 through 4, rows 1 through 3, of the B array 
would be set to the following computed values: 


Wo Wy Wo Zo 
Wa Wo W3 4 
Wa Ws Wa Z> 


This matrix represents equations (1), (2), and (3), the normal equa- 
tions for M = 2. 


The forward solution, which results in equations (4), (7), and (8) in 
the illustrative problem, is carried out by statements 23 through 31. 
By the time control has passed to statement 33, the coefficients of the 
AIT terms in the M+1 equations which would be obtained in hand calcula- 
tions have replaced the contents of the locations corresponding to 
columns 1 through M+t1, rows 1 through Mt1, of the B array, and the con- 
Stants on the right-hand side of the equations have replaced the con- 
tents of the locations corresponding to column M+2, rows 1 through Mti, 
of the B array. For the illustrative problem, columns 1 through 4, rows 
1 through 3, of the B array would be set to the following computed 
values: 


1 bio bis Day 
0 1 C23 C24 
0 0) C33 C34 


This matrix represents equations (4), (7), and (8). 


The backward solution, which results in equations (9), (10), and (11) 
in the illustrative problem, is carried out by statements 33 through 40. 
By the time control has passed to statement 41, which prints the values 
of the AY terms, the values of the (M+1) *AI terms have been stored in 
the M+1 locations for the A array. For the illustrative problem, the A 
array would contain the following computed values for az, ai, and ag, 
respectively: 


Location Contents 

A(3) C3u/7C33 

A(2) C24u~C23a2 

A(1) bay —bD42€84-b, 385 


The resulting values of the AI terms are then printed according to 
the FORMAT specification in statement 2. 
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APPENDIX _E:__DEPUG FACILITY 


The debug facility is a programming aid that enables the user to loc- 
ate errors in a FCRTRAN source program. It is available only with the 
Operating System (G) compiler and the Disk Operating System FORTRAN com- 
piler. The debug facility provides for tracing the flow within a pro- 
gram, tracing the flow between programs, displaying the values of vari- 
ables and arrays, and checking the validity of subscripts. | 


The debug facility consists of a DEBUG specification statement, an AT 
debug packet identification statement, and three executable statements. 
These statements, alone or in combination with any FORTRAN source lan- 
guage statements, are used to state the desired debugging operations for 
a Single program unit in source language. (A program unit is a single 
main program or a subprogranm.) 


The source deck arrangement consists of the source language state- 
ments that comprise the program, followed by the DEBUG specification 
statement, followed by the debug packets, followed by the END statement. 


The statements that make up a program debugging operation must be 
grouped in one or more debug packets. A debug packet is preceded by the 
AT debug packet identification statement and consists of one or more 
executable debug facility statements, and/or FORTRAN source language 
Statements. A debug packet is terminated by either another debug packet 
identification statement or the END statement of the program unit. 


PROGRAMMING CONSIDERATIONS 


The following precautions must be taken when setting up a debug 
packet: 


1. Any DO lcops initiated within a debug packet must be wholly con- 
tained within that packet. 


2e Statement numbers within a debug packet must be unique. They must 
be different from statement numbers within other debug packets and 
within the program being debugged. 


3. An error in a program should not be corrected with a debug packet; 
when the debug packet is removed, the error remains in the program. 
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4. The following statements must not appear in a debug packet: 


SUBROUTINE 

FUNCTION 

ENTRY 

IMPLICIT 

BLOCK DATA 

Statement function definition 


5. The program being debugged must not transfer control to any state- 
ment number defined in a debug packet; however, control may be 
returned to any point in the program from a packet. In addition, a 
debug packet may contain a RETURN, STOP, or CALL EXIT statement. 


DEBUG FACILITY STATEMENTS 


The specification statement (DEBUG) sets the conditions for operation 
of the debug facility and designates debugging operations that apply to 
the entire program unit (such as subscript checking). The debug packet 
identification statement (AT) identifies the beginning of the debug pac- 
ket and the point in the program at which debugging is to begin. The 
three executable statements (TRACE ON, TRACE OFF, and DISPLAY) designate 
actions to be taken at specific points in the program. The following 
text explains each debug facility statement and contains several pro- 
gramming examples. 
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DEBUG SPECIFICATION STATEMENT 


There must be one DEBUG statement for each program or subprogram to 
be debugged, and it must immediately precede the first debug packet. 


Where: option may be any of the following: 


UNIT (a) 
where a is an integer constant that represents a data set 
reference number. All debugging output is placed in this data 
set, called the debug output data set. If this option is es 
specified, any debugging output is placed in the standard out- 
put data set. All unit definitions within an executable ae 
gram must refer to the same unit. 


SUBCHK (na, MNazgeee,Nn) 
where n iS an array name. The validity of the subscripts used 
with the named arrays is checked by comparing the subscript 
combination with the size of the array. If the subscript 
exceeds its dimension bounds, a message is placed in the debug 
output data set. Program execution continues, using the aan 
correct subscript. If the list of array names is omitted, 
arrays in the program are checked for valid subscript Re 
If the entire option is omitted, no arrays are checked for 
valid subscripts. 


{ 

| 

| 

| 

: 

| 

| 

| 

| 

| 

| 

| 

l 

| 

| 

| 

| This option must be in the DEBUG specification statement of 
| each program or subprogram for which tracing is desired. If 
| this option is omitted, there can be no display of program 
] flow by statement number within this program. Even when this 
l option is used, a TRACE ON statement must appear in the first 
] debug packet in which tracing is desired. 

| 

| 

| 

| 

] 

| 

| 

| 

| 

: 

] 

| 

| 

| 

L 


INIT (m4 Magece in) 
where m is the name of a variable or an array that is to be 
displayed in the debug output data set only when the variable 
or the array values change. If mis a variable name, the name 
and value are displayed whenever “the variable is assigned a 
new value in either an assignment, a READ, or an assigned GO 
TO statement. If mis an array name, the changed element is 
displayed. If the list of names is omitted, a display occurs 
whenever the value of a variable or an array element is 
changed. If the entire option is omitted, no display occurs 
when values change. 


| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
all| 
| 
| 
| 
: | 
TRACE | 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
SUBTRACE | 
This option specifies that the name of this subprogram is to | 
be displayed whenever it is entered. The message RETURN is to] 
be displayed whenever execution of the subprogram is | 
completed. | 
The options in a DEBUG specification statement may be given in any 
order and they must be separated by commas 
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AT DEBUG PACKET IDENTIFICATION STATEMENT 


The AT statement identifies the beginning of a debug packet and indi- 
cates the point in the program at which debugging is to begin. There 
must be one AT statement for each debug packet; there may be many debug 
packets for one program or subprogram. 


|General Form | 
}---------------------------------------------------------------------- { 
JAT statement _numker | 
| 
|Where: statement number is an executable statement number in the pro-| 
| gram or subkprogram to be debugged. | 


The debugging operations specified within the debug packet are per- 
formed prior to the execution of the statement indicated by the state- 
ment number in the AT statement. 


TRACE ON STATEMENT 


The TRACE ON statement initiates the display of program flow by state- 
ment number. Each time a statement with an external statement number 
is executed, a record of the statement number is made on the debug out- 
put data set. This statement has no effect unless the TRACE option was 
specified in the DERUG specification statement. 


a ere ee eg eg a pO Ee eT NT ee 71 
|General Form | 

~=---=------------ == === === === - $= 2-2-2222 22 22222-22222 === === 4 
[TRACE ON | 
Me care ee as it ie ee ee ae ae ee J 


For a given debug packet, the TRACE ON statement takes effect immedi- 
ately before the execution of the statement specified in the AT state- 
ment; tracing continues until a TRACE OFF statement is encountered. The 
TRACE ON stays in effect through any level of subprogram call or return. 
However, if a TRACE ON statement is in effect and control is given to a 
program in which the TRACE option was not specified, the statement num- 
bers in that program are not traced. Trace output is placed in the 
debug output data set. 


This statement may not appear as the conditional part of a logical IF 
statement. 


TRACE OFF STATEMENT 


The TRACE OFF statement may appear anywhere within a debug packet and 
stops the recording of program flow by statement number. 


This statement may not appear as the conditional part of a logical IF 
Statement. 
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DISPLAY STATEMENT 


The DISPLAY statement may appear anywhere within a debug packet and 
causes data to be displayed in NAMELIST output format 


a aay re pe ry ere IG eT eS epee Tes eR pee ee en ee ee ee eer 7 
}General Form | 


| DISPLAY list 
jJWhere: list is a series of variable or array names, separated by 
| commas. 


The DISPLAY statement eliminates the need for FORMAT or NAMELIST and 
WRITE statements to display the results of a debugging operation. The 
data is placed in the debug output data set. 


The effect of a DISPLAY list statement is the same as the following 
FORTRAN IV source language statements: 


NAMELIST /name/list 
WRITE (n, name) 


where: 
name is the same in both statements. Note that subscripted 
variables or dummy arguments may not appear in the list. 


This statement may not appear as the conditional part of a logical IF 
statement. 


DEBUG PACKET PROGRAMMING EXAMPLES 


The following examples show the use of a debug packet to test the 
Operation of a program. 


Example 1: 


INTEGER SOLON, GFAR, EWELL 


10 SOLON = GFAR * SQRT (EWELL) 
11 IF (SOLON) 40, 50, 60 


DEBUG UNIT (3) 

AT 11 

DISPLAY GFAR, SOLON, EWELL 
END 


In example 1, the values of SOLON, GFAR, and EWELL are to be examined 


as they were at the completion of the arithmetic operation in statement 
10. Therefore, the statement number entered in the AT statement is 11. 
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The debugging oreration indicated is carried out just before execu- 


tion of statement 11 
Statement, the values of SOLON, 


were before execution of statement 10. 


DIMENSION STOCK(1000),0UT(1000) 


DO 30 I = 1, 1000 
25 STOCK (I) = STCCK (I) - OUT 
30 CONTINUE 
35 A=B+C 


DEBUG UNIT (3) 
AT 35 

DISPLAY STOCK 
END 


In example 2, all of the valves of STOCK are to be displayed. 


If statement number 10 is entered in the AT 
and EWELL are displayed as they 


When 


statement 35 is encountered, the debugging operation designated in the 


debug packet is executed. 
DO loop is displayed. 


Note: 


each iteration of the DO loop. 


Example 3: 

10 Az=1.5 

12 Le=i1 

15 B=A 4#+1.5 
20 DO 22 IT = 1,5 


22 CONTINUE 
25 € = B fF 3616 
30 D= C/2 

STOP 


DEBUG UNIT (3), TRACE 
Cc DERUG PACKET NUMBER 1 
AT 10 
TRACE ON 
Cc DEBUG PACKET NUMBER 2 
AT 20 
TRACE OFF 
DO 35 I = 1,3 


35 CONTINUE 
TRACE ON 
C DEBUG PACKET NUMBER 3 
AT 30 
TRACE OFF 
END 


Appendix E: 


Debug Facility 


The value of STOCK at the completion of the 


If the AT statement indicated statement 25 as the point of execu- 
tion for the debugging operation, the value of STOCK is displayed for 
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When statement 10 is encountered, tracing begins as indicated by 
debug packet 1. When statement 20 is encountered, tracing stops as 
indicated by the TRACE OFF statement in debug packet 2 and no tracing 
occurs during the execution of the statements within this packet. Trac- 
ing resumes before leaving debug packet 2 When statement 30 is encoun- 
tered, debug packet 3 is executed, and causes tracing to stop. 


In this example, all trace information is placed in the data set 
associated with data set reference number 3. This data set contains 
trace information for the following statement numbers: 10, 12, 15, 20, 
22, 22, 22, 22, 22, 25. Note that statement numbers 35 and 30 do not 
appear. 
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APPENDIX F: FORTRAN IV FEATURES NOT IN BASIC FORTRAN IV 


The following features in FORTRAN IV are not in Basic FORTRAN IV: 


ASSIGN 

BLOCK DATA 

Lakeled COMMON 

COMPLEX 

DATA 

Debug Facility 

More than three dimensions 

Object-time dimensions 

Object-time FORMAT specifications 

Assigned GO TO 

Logical IF 

LOGICAL 

PRINT b, list 

PUNCH b, list 

READ b, list 

END and ERR parameters in a READ 

Generalized Type statement (But 
provided as an explicit type.) 

IMPLICIT 

Call by name 

Literal as argument of CALL 

ENTRY 

RETURNi (i not a blank) 

NAMELIST 

PAUSE with literal 

G, Z, and L format codes 

Complex, logical, literal, 

Generalized subscript form 


note that DOUBLE PRECISION is 


and hexadecimal constants 


The following in-line subprograms in FORTRAN IV are not in BaSic FOR- 
TRAN IV: 


REAL NT 
AIMAG AINT 
DCMPLX IDINT 
CMPLX 

CCCNIG 

CONJG 

HEF IX 

CABS 

CDAEBS 


The following out-of-line subprograms in FORTRAN IV are not in Basic 
FORTRAN IV: 


CEXP DARSIN 
CDEXP ARCOS DARCOS 
CLOG TAN DTAN 
CDLOG COTAN DCOTAN 
CSIN SINH DSINH 
CDSIN COSH DCOSH 
CCOS ERF DERF 
CDCOS 

CSORT ERFC DERFC 
CDSORT GAMMA DGAMMA 
DATAN2 ALGAMA DLGAMA 
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APPENDIX G:__FORTRAN_IV_FEATURES_NOT_IN_USA_ FORTRAN IV 


Direct Access Input/Output Statements 
Double Exponentiation | 

END and ERR parameters in READ 

ENTRY 

Generalized subscripts 

Hexadecimal constant 

IMPLICIT 

Initial data values in type statement 
Length of variables as part of type specifications 
Literal enclosed in apostrophes 

Mixed mode expressions 

More than 3 dimensions in an array 
NAMELIST 

PAUSE ‘message’ 

PRINT 

PUNCH 

T and Z format codes 

RETURN i 
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SEND statement 48 


A format code 57 
ABS function 105 
absolute value functions 105 
actual arguments 85,91 
adjustable dimensions 96 
AIMAG function 106 
AINT runction 105 
ALGAMA function 105 
ALOG function 104 
ALOG10 function 104 
AMAXO function 105 
AMAX1 function 105 
AMINO function 106 
AMIN1 function 106 
AMOD function 105 
ARCOS function 104 
arccosine functions 104 
arcSsine functions 104 
arctangent functions 104 
arguments in function or subroutine 
subprograms 91 
arithmetic assignment statements 28 
arithmetic expressions 
defined 20 
order of computation 22 
arithmetic IF 34 
arithmetic operators 21 
arrays 
arrangement of 19 
dimension information 71 
general 18 
type specification 18 
ARSIN function 104 
ASSIGN and assigned GO TO 32 
assignment statements 28 
associated variable 63 
AT debug packet identification 117 
ATAN function 104 
ATAN2 function 104 


BACKSPACE statement 62 
Basic FORTRAN: IV 121 
basic real constant a fab 
blank common 77 

blank record 51 


blanks 9 
BLOCK DATA subprogram 99 
bytes (storage locations) 72 


CABS function 105 

call by name/location 92 

CALL statement 90 

Carriage control characters al 
~Ccos function 104 

CDABS function 105 

CDCOS function 104 

CDEXP function 104 


INDEX 


CDLOG function 104 
CDSIN function 104 
CDSORT function 104 
CEXP function 104 
character set 100 
character string 13,50,57 
CLOG function 104 
CMPLX function 106 
coding form 9 
coding statements 9 
comments 9 
common logarithm 104 
COMMON statement 76 
compilers 8 
completed error function 105 
COMPLEX statement 74 
complex values 
constants 12 
in arithmetic assignment statement 28 
in FORMAT statement 54 
length specification rz 
type specification 74 
computed GO TO 32 
CONJG function 106 
constants 10 
continuation statements 9 
CONTINUE statement 39 
control statements 31 
conversion rules 
in arithmetic assignment statements 30 
in FORMAT statements 53-60 
COS function 104 
COSH function 105 
COTAN function 104 
CSIN function 104 
CSORT function 104 


D format code a3 

DABS function 105 

DARCOS function 104 

DARSIN function 104 

DATA initialization statement 69 
in BLOCK DATA subprogram 99 

data set reference number 42 

DATAN function 104 

DATAN2 function 104 

DBLE function 105 

DCMPLX function 105 

DCONJG function 105 

DCOS function 104 

DCOSH function 105 

DCOTAN function 104 

debug facility 114 

DEBUG statement 116 

DEFINE FILE statement 62 

DERF function 105 

DERFC function 105 

DEXP function 104 

DFLOAT function 105 


Index 123 


DGAMMA function 105 
DIM function 105 
DIMENSION statement 71 
object-time dimensions 96 
direct access input/output statements 62 
programming considerations 64 
DISPLAY statement 118 
DLGAMA function 105 
DLOG function 104 
DLOG10 function 104 
DMAX1 function 105 
DMIN1 function 105 
DMOD function 105 
DO statement 36 
programming considerations 38 | 
double precision number (see real numbers) 
DOUBLE PRECISION statement 75 
DSIGN function 105 
DSIN function 104 
DSINH function 105 
DSORT function 104 
DTAN function 104 
DTANH function 105 
dummy arguments 85,89,91 
DUMP subprogram 107 
DVCHK subprogram 107 


E format code 53 
elements of the language 8 
embedded blanks 9 
END FILE statement 61 
END parameter in READ 44 
END statement 

in FUNCTION subprogram 88 

in main program 41 

in NAMELIST (&END) 48 
ENTRY statement 93 
equivalence groups 80,82 
EQUIVALENCE statement = 80 
ERF function 105 
ERFC function 105 
ERR parameter in READ 
error functions 105 
executable statement, definition 8 
EXIT subprogram 107 
EXP function 104 
explicit specification 17 
explicit specification statements 74 
exponential functions 104 
exponentiation 22 
expressions 

arithmetic 20 

defined 20 

logical 24 
extended range of DO 38 
EXTERNAL statement 95 


44,65 


F format code 53 
field descriptors 50 
FIND statement 67 
fix functions 106 
FLOAT function 106 
float functions 106 
FORMAT statement 
form 50 
purpose 51 
use at object time 60 
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formatted READ statement 45 


formatted records 43 

formatted WRITE statement 46 
FORTRAN coding form 9 
FORTRAN-supplied subprograms 103 
function definition 84 

function reference 84 

FUNCTION subprogram 86 


G format code 54 
GAMMA function 105 
gamma functions 105 
GO TO statement 
assigned 32 
computed 32 
unconditional 31 
group format specification 60 


H format code 58 
hexadecimal values 

constants 14 

transmitting aye) 
HFIX function 106 
hierarchy of operations 22 
hyperbolic cosine function 105 
hyperbolic sine function 105 
hyperbolic tangent function 105 


I format code 52 
IABS function 105 
IDIM function 106 
IDINT function 105 
IF statement 
arithmetic 34 
logical 35 
IFIX function 106 
implicit specification 17 
IMPLICIT statement 72 
implied DO 42 
index 123 
INIT option of DEBUG 116 
input/output statements 42 
INT function 105 
integers 
constants 10 
I format code 52 
length specification TZ 
magnitude 10 
type specification 74 
use in arithmetic assignment 
statements 28 
INTEGER statement 74 
I/O list 
defined 42 
omitted owl 


ISIGN function 106 


L format code 57 
labeled common 77 
language elements 8 
largest value functions 105 
length specification 16 
library subprograms 103 
literals 

constants 13 

data in FORMAT statements 58 
logical assignment statements 28 
logical expressions 24 


logical IF statement 35 
logical operators 25 
LOGICAL statement 74 
logical values 
constants 12 
type specification 74 
use in arithmetic assignment 
statements 28 
use in logical expressions 24 
logical variables 16 
loop control 36 
log-gamma functions 105 


mathematical subprograms 104 
MAXO function 105 

MAX1 function 105 

MINO function 106 

MIN1i function 106 

mixed mode expressions 30 

MOD function 105 

mode (see type) 

modular arithmetic functions 105 


NAMELIST statement 4] 
natural logarithm 104 
nested DO loops 38 
numeric format codes 54 


object-time dimensions 96 
object-time format 60 
operators 
arithmetic 21 
logical 25 
order of computation 22 
relational 24 
order 
of arithmetic computation 22 
of common blocks 78 
of equivalence groups 80 
of logical expression computation 26 
of source program statements 9 
OVERFL subprogram 107 


P scale factor 56 
parentheses 
in arithmetic expressions 22 
in logical expressions 27 
in FORMAT statement 51 
PAUSE statement 40 
PDUMP subprogram 107 
positive difference functions 106 
predefined specifications 17 
primary 
arithmetic 21 
logical 24 
PRINT b, list 102 
PRINT statement 102 
printer control characters 51 
program unit, definition 8 
PUNCH b, list 101 
PUNCH statement 101 





range of DO 36, 38 
READ b, list 101 
READ statement 
direct access 65 
sequential 44,101 
REAL function 106 


real numbers 
constants 11 
in D, E, and F format codes 53 
length specification 72 
magnitude 11 
precision 11 
type specification 73 
use in arithmetic assignment 
statements 28 
REAL statement 74 
record number 62 
records 
formatted 43 
length of 63 
unformatted 43 
reference by location 92 
reference by value 92 
relational operators 24 
RETURN statement 
in FUNCTION subprogram 88 
in main program 90 
in SUBROUTINE subprogram 90 
REWIND statement 61 


scale factor 56 

sequential input/output Guy 
service subprograms 107 

SIGN function 106 

Sign transfer functions 106 
SIN function 105 

SINH function 105 

size specification, array 18 
SLITE subprogram 107 

SLITET subprogram 107 
smallest value functions 106 
SNGL function 106 

source program characters 100 
special characters 100 
specification statements 71 
SORT function 104 


square root functions 104 
statement 

categories 8 

function definitions 84 


numbers 9 
order 9 


source 8 
STOP statement Ky 
storage locations (bytes) 72 


for literals 13 
SUBCHK debug option 116 
subprograms 

arguments 91 

BLOCK DATA 99 

FUNCTION 86 

general 83-99 

multiple entry 93 

naming 83 
subprogram statements 83 
SUBROUTINE subprogram 88 
subscripts 19 
SUBTRACE debug option 116 
symbolic names 14 


T format code 59 

TAN function 104 

TANH function 105 
termination of program G1 


Index 125 


TRACE OFF statement 117 
TRACE ON statement 117 
TRACE debug option 116 
transfer of sign functions 106 
trigonometric cosine functions 104 
trigonometric cotangent functions 104 
trigonometric sine functions 104 
trigonometric tangent functions 104 
truncation functions 105 
truth values 13920 
type specification . 
of arithmetic expressions 23 
of arrays 18 
of FUNCTION subprogram 83 
of statement function definitions 83 
of variables 16 
type statement 74 
type statements 71 


unary minus 22 


unconditional GO TO 31 
unformatted READ statement 45 


126 


unformatted records 43 
unformatted WRITE statement 47 
UNIT debug option 116 

USA FORTRAN IV 7,122 


variable format statements 60 
variables 
arrangement in common 79 
arrangement in equivalence groups 82 
general 15 
length specification 16 
names 15 
type specification 16 


WRITE statement 
direct access 66 
sequential 46 

X format code 59 


Z format code 53 
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